Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,621,1180AC23.7% 72.1 / 10.4 38V7G (GTG→GGG) eamAcysteine/O‑acetylserine exporter EamA
Reads supporting (aligned to +/- strand):  ref base A (22/7);  new base C (0/9);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 7.02e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTA  >  NC_000913/1620982‑1621255
                                                                                                                                        |                                                                                                                                         
tGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAccc                                                                                                                                        >  1:770054/1‑140 (MQ=255)
tGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAccc                                                                                                                                        >  2:286259/1‑140 (MQ=255)
          aaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcg                                                                                                                              <  1:699450/140‑1 (MQ=255)
          aaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcg                                                                                                                              <  1:376373/140‑1 (MQ=255)
            gATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgcg                                                                                                                            <  2:64003/140‑1 (MQ=255)
                      aaaGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTAcc                                                                                                                  >  2:435624/1‑140 (MQ=255)
                        aGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGc                                                                                                                <  1:316480/140‑1 (MQ=255)
                                   tAAAGCGCAAATCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTc                                                                                                     >  2:676024/1‑140 (MQ=255)
                                           aaaCCGGCCAGCATCAGCGGTGGCATGTTATGACGCCCCACTTTGAGGACCACAAAATTTAGACCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCttttt                                                                                             >  2:592589/1‑139 (MQ=255)
                                            acccGGCCAGCATCGGCGGTGGCATGTTATGAACCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                            <  2:802191/138‑1 (MQ=255)
                                            aaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGCCCCCTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                            <  1:1169109/140‑1 (MQ=255)
                                             aCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                           >  2:739578/1‑140 (MQ=255)
                                                 gCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                       <  1:850790/140‑1 (MQ=255)
                                                 gCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                       >  2:531706/1‑140 (MQ=255)
                                                      cATCAGCGGTGGCATTTTATGAAGCCCCATTTTGATGACCAAAAAATTTAGCCCCCATACGCCCACTACCAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCaa                                                                                  <  1:286259/140‑1 (MQ=255)
                                                           gCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAg                                                                             <  1:676024/140‑1 (MQ=255)
                                                             ggTGGCATGTTATGAACCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACCACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                           <  2:465113/140‑1 (MQ=255)
                                                                     tttATGAAGCCCCATTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                   <  2:496515/139‑1 (MQ=255)
                                                                     gTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                   <  1:755820/140‑1 (MQ=255)
                                                                            aGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCtt                                                            >  2:869957/1‑140 (MQ=255)
                                                                              ccccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTTAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTa                                                          >  1:477673/1‑140 (MQ=255)
                                                                                ccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                        >  1:1001967/1‑140 (MQ=255)
                                                                                      tgatgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGa                                                  <  1:739578/140‑1 (MQ=255)
                                                                                           cccACAAAATTTACCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGt                                             <  2:586295/139‑1 (MQ=255)
                                                                                                  aaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGc                                      >  1:1033208/1‑140 (MQ=255)
                                                                                                       aGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                 >  2:755848/1‑140 (MQ=255)
                                                                                                           cccATGCAACCACTACCAGTAGCGCCACCCCCCCATTTTTTCGCAAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacga                             <  2:552765/140‑1 (MQ=255)
                                                                                                                       cTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                 >  1:664495/1‑140 (MQ=255)
                                                                                                                       cTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                 >  2:267004/1‑140 (MQ=255)
                                                                                                                        tACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcg                >  2:953198/1‑140 (MQ=255)
                                                                                                                        tACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcg                >  1:574019/1‑140 (MQ=255)
                                                                                                                         aCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcgc               >  1:1124745/1‑140 (MQ=255)
                                                                                                                            aGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaaa            >  2:1118522/1‑140 (MQ=255)
                                                                                                                              tAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATa          <  1:755848/140‑1 (MQ=255)
                                                                                                                              tAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATa          <  2:664495/140‑1 (MQ=255)
                                                                                                                              tAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATa          >  1:306611/1‑140 (MQ=255)
                                                                                                                                gcgcCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAg        >  2:32069/1‑140 (MQ=255)
                                                                                                                                gcgcCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAg        >  1:603915/1‑140 (MQ=255)
                                                                                                                                  gcCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc      >  2:407342/1‑140 (MQ=255)
                                                                                                                                   ccAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg     >  1:658837/1‑140 (MQ=255)
                                                                                                                                      acacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCta  >  2:224405/1‑140 (MQ=255)
                                                                                                                                      acacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt   >  1:694357/1‑139 (MQ=255)
                                                                                                                                      acacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt   <  2:694357/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                         
TGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTA  >  NC_000913/1620982‑1621255

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