Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 1,812,261 C→A 8.4% T150T (ACC→ACA)  rpoS → RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,812,2610CA 8.4% 139.2 / 2.0 46T150T (ACC→ACA) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base A (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

GCGCCTGGTTGTCAAAATTGCCCGTCGATACGTGAACCGTGGCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGAACTGAACGTCTACCTGCGCGCCGCGCGCGAGCTGAC  >  ACB143/1812126‑1812386
                                                                                                                                       |                                                                                                                             
gcgcCTGGTTGTCAAAATTGCCCGTCGATACGTGAACCGTGGCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCGGGGGCTGAACCGGGCGTTCGAGAAGTTCGACCAGGAGCTTGGTTTACGGTTTTCGAcc                                                                                                                               >  8:141152/1‑136 (MQ=255)
 cgcCTGTTTTTCAAAATTGCCCGTCGATACGTGAACCGTGGCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGAcc                                                                                                                               <  8:157570/135‑1 (MQ=255)
   ccTGGTTGTCAAAATTGCCCGTCGATACGTGAACCGTGGCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTAt                                                                                                                            <  1:34473/136‑1 (MQ=255)
         tGTCAAAATTGCCCGTCGATACGTGAACCGTGGCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGAcc                                                                                                                      >  2:35549/1‑136 (MQ=255)
            cAAAATTGCCCGTCGATACGTGAACCGTTGCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCtgg                                                                                                                   <  3:83263/136‑1 (MQ=255)
                       gTCGATACGTGAACCTTGGCCTGTCGTTGCTAGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCTGGCGGTCGAGAAGTTCGACCCGGAGCGTTGTTTCCGTTTTTCGACCTATGCTACCTGGTGGATTCGCCa                                                                                                        <  6:303504/136‑1 (MQ=255)
                                   aCCGTTGCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACg                                                                                            <  8:108224/136‑1 (MQ=255)
                                   aCCGTGGCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACg                                                                                            <  1:378131/136‑1 (MQ=255)
                                      gTGGCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGc                                                                                         >  3:55925/1‑136 (MQ=255)
                                         gCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGa                                                                                       <  3:355208/135‑1 (MQ=255)
                                         gCCTGTCGTTGCTCGACATGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCTACCCGGAGCGTGGTTTCAGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGCACGGGCGAt                                                                                      >  7:272148/1‑136 (MQ=255)
                                          ccTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGAt                                                                                      >  7:152857/1‑135 (MQ=255)
                                          ccTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGAt                                                                                      >  2:332807/1‑135 (MQ=255)
                                            tGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGACGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGCCCTATGCGACCTGGTGGATTCGCCAGACCCTCGAACGGGCGAGCAt                                                                                   >  7:107512/1‑136 (MQ=255)
                                                 ttGCTGGACCTTATCGAAGAGGGCAACCTGGGGCTGATCCGCGCGGTAGAGAAGTTTGACCCGGAACGTGGTTTCCGCTTCTCAACATACGCAACCTGGTGGATTCGCCAGACGATTGAACGGGCGATTATGAAc                                                                               >  6:169965/1‑135 (MQ=255)
                                                    ctggACCTTATCGAAGAGGGCAACCTGGGGCTGATCCGCGCGGTAGAGAAGTTTGACCCGGAACGTGGTTTCCGCTTCTCAACATACGCAACCTGGTGGATTCGCCAGACGATTGAACGGGCGATTATGAACCaa                                                                            <  5:169965/132‑2 (MQ=255)
                                                    cTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGTTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAg                                                                            <  5:88232/135‑1 (MQ=255)
                                                       gACCTGATCGAGGAGTGCAACCTGGGGCTGATCCTGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGAc                                                                          <  8:96865/134‑1 (MQ=255)
                                                         ccTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCg                                                                       >  5:310970/1‑135 (MQ=255)
                                                            gATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGGTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATAATTAACCAGACCCGCAc                                                                    >  6:31041/1‑135 (MQ=255)
                                                                      ggCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGc                                                         <  7:252070/136‑1 (MQ=255)
                                                                      ggCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGc                                                         >  1:239904/1‑136 (MQ=255)
                                                                         aaCCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGAGCATGAACCAGACCCGCACCATTCGCCTGcc                                                        >  4:73584/1‑134 (MQ=255)
                                                                           ccTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATc                                                    <  2:271338/136‑1 (MQ=255)
                                                                           ccTGGGGCTGATCCGGGAGGTCTTGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATc                                                    <  2:223031/136‑1 (MQ=255)
                                                                               gggCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCa                                                  <  2:16686/134‑1 (MQ=255)
                                                                                   tGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCCCGTGGt                                            >  8:96675/1‑136 (MQ=255)
                                                                                        cgcgCGGTAGAGAAGTTTGACCCGGAACGTGGTTTCCGCTTCTCAACATACGCAACCTGGTGGATTCGCCAGACGATTGAACGGGCGATTATGAACCAAACCCGTACTATTCGTTTGCCGATTCACATCGTAAAg                                        >  5:278024/4‑135 (MQ=255)
                                                                                        cGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTGCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGCCCATCGAACGGGCGAACATGAACCAGACCCGCCCCATTCGCCTGCCGCTCCACGTGGTCAAgg                                       >  6:223500/1‑136 (MQ=255)
                                                                                        cGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAg                                        >  8:144693/1‑135 (MQ=255)
                                                                                        cGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAg                                        >  3:83324/1‑135 (MQ=255)
                                                                                        cGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAg                                        >  2:330322/1‑135 (MQ=255)
                                                                                        cGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAg                                        >  3:253612/1‑135 (MQ=255)
                                                                                        cGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAg                                        >  2:264314/1‑135 (MQ=255)
                                                                                         gggCGGTCTAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGa                                      <  2:340843/136‑1 (MQ=255)
                                                                                         gggCGGTCTAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACATATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGa                                      <  7:132573/136‑1 (MQ=255)
                                                                                         gggCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGa                                      <  1:332807/136‑1 (MQ=255)
                                                                                           gCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGAAc                                    <  6:171858/136‑1 (MQ=255)
                                                                                           gCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGAAc                                    <  7:180590/136‑1 (MQ=255)
                                                                                           gCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTACGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGa                                      <  3:262083/134‑1 (MQ=255)
                                                                                           gCGGTCGAGAAGTTCGACACGTAGCGTGGTTTCCGTTTTTCGAACTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATTAAACAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGAAc                                    <  5:366765/136‑1 (MQ=255)
                                                                                             ggTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGCTCCACGTGGTCAAGGAACt                                   >  2:362009/1‑135 (MQ=255)
                                                                                                  agaAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCCAACGGGCGATCATGAACCAGCCCCGCACCATACGCCTGCCGATCCACGTGGTCAAGGAACTGCACGt                             >  1:307534/1‑136 (MQ=255)
                                                                                                  agaAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACAATCGAACGGGCGATCATGAACCAGACACGCACCATTCGCCTGCCGATCCACGTGGTCAAGGAACTGAACGt                             >  2:261974/1‑136 (MQ=255)
                                                                                                            cccGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGAACTGAACGTCTACCTgcgc                   >  5:308735/1‑136 (MQ=255)
                                                                                                                             gTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGAACTGAACGTCTACCTGCGCGCCGCGCGCGAGCTGAc  <  5:223500/136‑1 (MQ=255)
                                                                                                                             gTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGAACTGAACGTCTACCTGCGCGCCGCGCGCGAGCTGAc  <  5:37145/136‑1 (MQ=255)
                                                                                                                             gTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGAACTGAACGTCTACCTGCGCGCCGCGCGCGAGCTGAc  <  1:264314/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                             
GCGCCTGGTTGTCAAAATTGCCCGTCGATACGTGAACCGTGGCCTGTCGTTGCTCGACCTGATCGAGGAGGGCAACCTGGGGCTGATCCGGGCGGTCGAGAAGTTCGACCCGGAGCGTGGTTTCCGTTTTTCGACCTATGCGACCTGGTGGATTCGCCAGACCATCGAACGGGCGATCATGAACCAGACCCGCACCATTCGCCTGCCGATCCACGTGGTCAAGGAACTGAACGTCTACCTGCGCGCCGCGCGCGAGCTGAC  >  ACB143/1812126‑1812386

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑