Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA1,341,7321.A30.8% 13.8 / 14.8 13intergenic (+5/‑189)yciM/pyrFLPS regulatory protein; putative modulator of LpxC proteolysis/orotidine‑5'‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base . (4/5);  new base A (2/2);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGCCGTTTACC  >  NC_000913_3_pae_tpiA/1341584‑1341876
                                                                                                                                                     |                                                                                                                                                
gACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTa                                                                                                                                                 >  2:1573/1‑150 (MQ=255)
                                                                   acacCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATa                                                                                                                     <  3:45769/111‑1 (MQ=255)
                                                                   acacCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATa                                                                                                                     >  4:45769/1‑111 (MQ=255)
                                                                      cccTCTACCGGCATTGTCCGTCTTGTCTTGTCTTGGCATTCTTTAAGCGTATTCGCGTTCTTGATTGCCTTTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCa                                                                           <  4:54740/150‑1 (MQ=255)
                                                                              tGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTaaaaaaaa                                                                                                                                          >  1:113731/1‑79 (MQ=255)
                                                                              tGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTaaaaaaaa                                                                                                                                          <  2:113731/79‑1 (MQ=255)
                                                                                gCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCAtt                                                                                                 >  1:134872/1‑119 (MQ=255)
                                                                                gCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCAtt                                                                                                 <  2:134872/119‑1 (MQ=255)
                                                                                     gTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTc                                                                                           >  1:3384/1‑119 (MQ=255)
                                                                                     gTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTc                                                                                           <  2:3384/119‑1 (MQ=255)
                                                                                                                                  ttGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATg               >  2:45346/1‑150 (MQ=255)
                                                                                                                                         ccTGTAATTTTTAAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTc                            >  1:110326/1‑131 (MQ=255)
                                                                                                                                         ccTGTAATTTTTAAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTc                            <  2:110326/131‑1 (MQ=255)
                                                                                                                                             tAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGTAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGCCGTTTa    >  4:35893/1‑150 (MQ=255)
                                                                                                                                              aaTTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGCCGTTTAc   <  1:59029/150‑1 (MQ=255)
                                                                                                                                               aTTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGTAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGCCGTTTAcc  <  3:35893/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGCCGTTTACC  >  NC_000913_3_pae_tpiA/1341584‑1341876

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