Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA2,637,4440AG26.1% 31.2 / 11.3 25intergenic (‑88/+30)der/bamBGTPase; multicopy suppressor of ftsJ/BamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (8/9);  new base G (2/4);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTT  >  NC_000913_3_pae_tpiA/2637297‑2637574
                                                                                                                                                   |                                                                                                                                  
gTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGc                                                                                                                                  >  1:53309/1‑150 (MQ=255)
     aTAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCTCCTGGACAGGGGCCGttt                                                                                                                             <  2:57173/150‑1 (MQ=255)
          gTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCTCCTGGACAGGGGCCGTTTTcctg                                                                                                                        <  1:22085/150‑4 (MQ=255)
                                gcCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                        <  1:158704/128‑1 (MQ=255)
                                gcCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCaaa                                                                                                                        >  2:158704/1‑128 (MQ=255)
                                  ccGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgt                                                                                                <  3:57998/150‑1 (MQ=255)
                                   cGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgtg                                                                                               <  1:23127/150‑1 (MQ=255)
                                    gACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACgtgt                                                                                              <  2:53309/150‑1 (MQ=255)
                                               ccACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCTCCTGGACAGGGGCCGTTTTcct                                                                                                                         >  1:93138/1‑110 (MQ=38)
                                               ccACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCTCCTGGACAGGGGCCGTTTTcct                                                                                                                         <  2:93138/112‑3 (MQ=38)
                                               ccACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGAc                                                                                                         >  1:49845/1‑128 (MQ=255)
                                               ccACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGAc                                                                                                         <  2:49845/128‑1 (MQ=255)
                                                  caGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGAGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCTGGTGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACgg                                                                                <  3:23025/150‑1 (MQ=255)
                                                               tAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGc                                                                   <  3:110028/150‑1 (MQ=255)
                                                                      tCATTTTTATACATCATCAGCCAACGCACGCTATATTCACGTCGATGTTAAAAACAGGAAAACGGCCCCTGTCCAGATGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGCACACGGTTCCGTCTTTTGCCTGGATc                                                            >  2:165844/1‑150 (MQ=255)
                                                                      tCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATc                                                            >  1:159068/1‑150 (MQ=255)
                                                                                 aaTCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGcc                                                 <  2:60827/150‑1 (MQ=255)
                                                                                                    cTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCt                                                                 <  2:81757/115‑1 (MQ=255)
                                                                                                    cTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCt                                                                 >  1:81757/1‑115 (MQ=255)
                                                                                                    cTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTg                                                   >  3:68430/1‑129 (MQ=255)
                                                                                                    cTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTg                                                   <  4:68430/129‑1 (MQ=255)
                                                                                                                     ttAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACt             >  4:31924/1‑150 (MQ=255)
                                                                                                                                aaaaCGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACttttt  >  2:46895/1‑150 (MQ=255)
                                                                                                                                  aaCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCt                          <  2:175736/124‑1 (MQ=255)
                                                                                                                                  aaCGGCTCCTGGACAGGGGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCt                          >  1:175736/1‑124 (MQ=255)
                                                                                                                                                   |                                                                                                                                  
GTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAATCATCAGACAACGCACGCTATATTCGCGTCGTTGTTAAAAACAGGAAAACGGCCCCTGTCCAGGAGCCGTTTTCAAAGTGAACGACAGAGACGATTAACGTGTAATAGAGTACACGGTTCCGTCTTTTGCCTGGATCAGCAGTTTGCCGTCAGCGGCAACCGGTTCAGTCTGGAAACCGGAACTATCAACTTTTT  >  NC_000913_3_pae_tpiA/2637297‑2637574

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