Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,637,439 C→A 28.6% intergenic (‑83/+35) der ← / ← bamB GTPase; multicopy suppressor of ftsJ/BamABCDE complex OM biogenesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,637,4390CA28.6% 37.2 / 15.1 21intergenic (‑83/+35)der/bamBGTPase; multicopy suppressor of ftsJ/BamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (8/7);  new base A (3/3);  total (11/10)
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.63e-01

AGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGC  >  NC_000913/2637363‑2637581
                                                                            |                                                                                                                                              
aGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAgaga                                                                                                                 >  4:773641/1‑108 (MQ=255)
aGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAgaga                                                                                                                 <  3:773641/108‑1 (MQ=255)
                                cgcTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATag                                                                                                <  3:408050/93‑1 (MQ=255)
                                cgcTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATag                                                                                                >  4:408050/1‑93 (MQ=255)
                                    atatTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                               >  1:1013108/1‑58 (MQ=255)
                                    atatTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                               <  2:1013108/58‑1 (MQ=255)
                                    atatTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGa                                                                     >  4:848760/1‑116 (MQ=255)
                                    atatTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGa                                                                     <  3:848760/116‑1 (MQ=255)
                                                                aaCGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATc                                                                   >  4:807791/1‑90 (MQ=255)
                                                                aaCGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATc                                                                   <  3:807791/90‑1 (MQ=255)
                                                                 aCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCttt                                                                             >  4:914221/1‑79 (MQ=37)
                                                                 aCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCttt                                                                             <  3:914221/79‑1 (MQ=37)
                                                                 aCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGa                       >  3:876562/1‑133 (MQ=255)
                                                                 aCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGa                       <  4:876562/133‑1 (MQ=255)
                                                                 aCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACtttt          <  3:44783/146‑1 (MQ=255)
                                                                 aCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACtttt          >  4:44783/1‑146 (MQ=255)
                                                                    gCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATagag                                                                                              >  4:440179/1‑59 (MQ=255)
                                                                    gCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATagag                                                                                              <  3:440179/59‑1 (MQ=255)
                                                                       ccTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTg                                                          <  2:206563/92‑1 (MQ=255)
                                                                       ccTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTg                                                          >  1:206563/1‑92 (MQ=255)
                                                                       ccTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGc  >  1:596532/1‑148 (MQ=255)
                                                                            |                                                                                                                                              
AGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTTGCTGGGC  >  NC_000913/2637363‑2637581

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