Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,638,026 (C)6→7 coding (627/1179 nt) bamB ← outer membrane protein assembly factor BamB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,638,0201.C90.9% 101.0 / 3.6 33coding (633/1179 nt)bamBouter membrane protein assembly factor BamB
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (15/15);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CAACGCCGTTAACAACGACGGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTGACCGTTAC  >  NC_000913/2637881‑2638153
                                                                                                                                             |                                                                                                                                       
caacGCCGTTAACAACGACGGGAGTCGAGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCg                                                                                                                                          <  2:124331/140‑1 (MQ=255)
    gCCGTTAACAACGACGGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGcccc                                                                                                                                    >  2:508699/1‑140 (MQ=255)
         tAACAACGACGGGAATCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGcccccccac                                                                                                                              >  2:905769/1‑137 (MQ=255)
                 acgGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGAACAGTATCGCCCCCCCACGACGGCCg                                                                                                                       >  1:350130/1‑140 (MQ=255)
                    ggAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCAc                                                                                                                    >  2:30825/1‑140 (MQ=255)
                     gAGTCGTGTCAACATCGGTCAGACGGTCAATCTCGTTAGTACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCAcc                                                                                                                   <  2:409375/140‑1 (MQ=255)
                      aGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCa                                                                                                                  >  1:560964/1‑140 (MQ=255)
                       gTCGTGTCAACATCGCTCAGACGGTCAATTTTGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCaa                                                                                                                 <  2:360656/140‑1 (MQ=255)
                       gTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCaa                                                                                                                 <  1:228766/140‑1 (MQ=255)
                                 cATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGt                                                                                                       <  1:656717/140‑1 (MQ=255)
                                      cTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGGCCCCCCAGGACCGCAGCACCAAAAGCCGCTGTCGgcg                                                                                                  >  1:386204/1‑140 (MQ=255)
                                      cTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGgcg                                                                                                  >  2:892830/1‑140 (MQ=255)
                                          gACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCTCCACGACAGCCGCATCAAAAGCCGTTGTCGGCGCAGa                                                                                              >  1:416036/1‑140 (MQ=255)
                                               tCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGc                                                                                         <  1:58620/140‑1 (MQ=255)
                                               gcAATTTCGGTAGAACTCGTCGCCTGGGAAATATGCCGCTGCCAAATCGTCCGGCCCTGTTCCCACCGCCCTGGGCTGGCGCGGCCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGc                                                                                         <  2:448858/139‑1 (MQ=255)
                                                              ccGGTCGCCTGGGAAATACGCTCCTGCCAAATCATATGGCCCTGTTCCTTCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCg                                                                          <  2:368077/140‑1 (MQ=255)
                                                                 gTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAg                                                                       <  1:508699/140‑1 (MQ=255)
                                                                   cGCCTGGGAAATACGCTGCTGCCAAATCATCTAGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGc                                                                     >  2:162638/1‑140 (MQ=255)
                                                                        gggAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATc                                                                <  1:30825/140‑1 (MQ=255)
                                                                                            atcatcTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGa                                            <  1:892830/140‑1 (MQ=255)
                                                                                                          ttCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCt                              >  1:778933/1‑140 (MQ=255)
                                                                                                               tCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGGGGTTAACGGTCTATATGACAGGGCCGTCAGCTTCGtt                         >  1:267010/1‑140 (MQ=255)
                                                                                                               tCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGtt                         >  1:839334/1‑140 (MQ=255)
                                                                                                               tCAGCACTGCGCTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACGGTCCATTTGACAGGGCCGTCAGCTTCGtt                         >  1:337229/1‑140 (MQ=255)
                                                                                                                       gcgcTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCtt                 >  2:618101/1‑140 (MQ=255)
                                                                                                                        cgcTGACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTg                >  1:955324/1‑140 (MQ=255)
                                                                                                                        cgcTGACGCGACCATTATCGCCACCCCACGACGGCTGCATCAAAAGCCATTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTg                <  2:416036/140‑1 (MQ=255)
                                                                                                                            gACGCGACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAAc            <  2:955324/140‑1 (MQ=255)
                                                                                                                              cgcgACCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTg          <  2:560964/140‑1 (MQ=255)
                                                                                                                                  aCCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAgcgc                   >  1:49881/1‑127 (MQ=255)
                                                                                                                                  aCCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAgcgc                   <  2:49881/127‑1 (MQ=255)
                                                                                                                                  aCCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTGACCg      >  1:139456/1‑140 (MQ=255)
                                                                                                                                  aCCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTGACCg      >  1:68286/1‑140 (MQ=255)
                                                                                                                                  aCCATTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTCGTAACTGACCg      >  2:912659/1‑140 (MQ=255)
                                                                                                                                     aTTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTGACCGTTa   <  1:905769/140‑1 (MQ=255)
                                                                                                                                     aTTATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTGACCGTTa   <  1:138053/140‑1 (MQ=255)
                                                                                                                                      ttATCGCCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTGACCGTTAc  <  2:148746/140‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                       
CAACGCCGTTAACAACGACGGGAGTCGTGTCAACATCGCTCAGACGGTCAATTTCGGTAGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGCTTGTAACTGACCGTTAC  >  NC_000913/2637881‑2638153

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: