Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,265,6170CA43.8% 51.7 / 13.6 34F56L (TTC→TTAyeiPelongation factor P‑like protein
Reads supporting (aligned to +/- strand):  ref base C (17/1);  new base A (0/14);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 3.18e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTATGGTACTGAATTACAACGGCAAACTGCTGCTGGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAAAGAAGACTATACCCCGTATACCTTCACCAAAGATCAGATTGAAGAAGAGTTGCTGT  >  NC_000913/2265476‑2265765
                                                                                                                                              |                                                                                                                                                     
aGGTATGGTACTGAATTACAACGGCAAACTGCTGCTGGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGtgat                                                                                                                                              >  2:39098/1‑151 (MQ=255)
    aTGGTACTGAATTACAACGGCAAACTGCTGCTGGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTgatgat                                                                                                                                          >  2:58793/1‑150 (MQ=255)
       gTACTGAATTACAACGGCAAACTGCTGCTGGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCg                                                                                                                                      >  2:554615/1‑151 (MQ=255)
            gAATTACAACGGCAAACTGCTGCTGGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGac                                                                                                                                 >  2:307044/1‑151 (MQ=255)
             aaTTACAACGGCAAACTGCTGCTGGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGac                                                                                                                                 >  2:331209/1‑150 (MQ=255)
             aaTTACAACGGCAAACTGCTGCTGGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGac                                                                                                                                 >  2:130244/1‑150 (MQ=255)
                    aCGGCAAACTGCTGCTGGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGAc                                                                                                                         <  2:477361/151‑1 (MQ=255)
                      ggCAAACTGCTGCTGGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGc                                                                                                                       >  1:204271/1‑151 (MQ=255)
                             tgctgctgGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACccg                                                                                                                >  2:232117/1‑151 (MQ=255)
                                   tgGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGCGCCGATACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACccgccg                                                                                                             >  1:619578/1‑148 (MQ=255)
                                                  atattCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACCCGCCGTTACGTTGACTTCTCCt                                                                                            >  2:64523/1‑150 (MQ=255)
                                                  atattCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACCCGCCGTTACGTTGACTTCTCCt                                                                                            >  2:106749/1‑150 (MQ=255)
                                                           ccccccCTGCCCGCGGTGCAGCTACGCTGAATAAAGGGCTTTTTTGTAATGCCGGTCCCGGGATAAAAGAAGAAGGGTGTTAAAAGGGTGATAATATTGTTGACCCGGGGCCGCTCCCCCGCCGTTACGTTGACTTCTCCTATGTCg                                                                                      <  2:429078/145‑1 (MQ=255)
                                                           cGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATgg                                                                                  >  1:52980/1‑151 (MQ=255)
                                                           cGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATgg                                                                                  >  1:86373/1‑151 (MQ=255)
                                                            gCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGc                                                                                 >  2:30546/1‑151 (MQ=255)
                                                             ccccCTGCCGGGGGTGCCGCTACGTTGAAAAAATTGTTTTTTTCTGAGGCCCGCGCCGGGCAGAAAGTAGAGGTGTGTTAAAAGGGTGATGTTATCGTTGCCACCGTGACGCTGCCCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCa                                                                                <  2:435023/151‑1 (MQ=255)
                                                             cccGCTGCCCGCGGGCCCGCTCCGTTATAAAAAATGCGTTTTTTGGTGGCCGGTGCGGGGTTAAAAGAGGAGGTGTTTTAAAAGGGTGATGATGGTGTTGACGCGGGGCCGCTGAGCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCa                                                                                <  1:24514/151‑1 (MQ=255)
                                                              cccctcccgggggTGCCGCTAGGTTAAAAAAATTGTTTTTTTCTGATGTCCGTCCGGGGCAAAAAGTAAAAGGGTTTTAAAAAGGTGATGATATCTTTGACCCCGTGACGCTGCCCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCaa                                                                               <  1:673585/140‑1 (MQ=255)
                                                               cccgccccggggccccGCTACGGTAAAAAAATTGTTTTTTTTTTATGCCCGCCCCGGGATAAAAGAAGAGGTGTTTTTAAAGGGTAATTATTTCTTTCACCCCGGCCCCCTCACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAAc                                                                              <  2:304842/137‑1 (MQ=255)
                                                                       ggggTGCCGCTACGCTGTAAAAAATGCTTTTTTCTGATGCCCGTCCGGGGCTAAAAGAAGAAGAGCTTTAAAAAGGTGATGATATCGTTGACCCCGTGACGCTGCCCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATg                                                                       <  1:297156/148‑1 (MQ=255)
                                                                          gtcccGCTGCGCTGTAAAAAATGCTTTTTTCTGAGGTCCGTCCGGGGCTAAAAGAAGAAGAGCTTTTAAAGGGTATTGTTATCGTTGCCCCCGTGACGCGGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCt                                                                    <  1:470361/147‑1 (MQ=255)
                                                                             ccccTACGTTGTAAAAATTGTTTTTTTCGTATGCCCGTCCGGGGCAAAAAGTAGAAGAGGTTTTAAAGGTTTTTATTTTCTTTGCCCCCGGGACCCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTa                                                                 <  2:494573/147‑1 (MQ=255)
                                                                             ccGCTACGCTGAACAAAGTGTTTTTTTCTGATGTCCGTGGCGGGATAAAAGGAGAAGAGTTTTTAAAAGTTGATGATTTCTTTGACCCCGGGACCGTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCt                                                                    <  2:350371/147‑1 (MQ=255)
                                                                                    gCTGTAAAAAATGCTTTTTTCTGATGTCCGTCCCGGGCAAAAAGTAGAAGAGCTTTTAAAGGGTGATGATTTGGTGGACACGGGGAGGCTGCCCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGa                                                             <  1:671834/147‑1 (MQ=255)
                                                                                         aaaaagTGTGTTTTTCTGTTGTCCGTCCCGGGTAGAAAGTAGAAGGGTGTTAAAAGGGTAATATTTTCTTAGACCCGGTGACGCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAaaga                                                       <  1:296951/146‑1 (MQ=255)
                                                                                                  gTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAAAGAAGACTATAcccc                                           >  2:574046/1‑151 (MQ=255)
                                                                                                       ttggATGCCCGCCCCGGGCTAAAAGTAGAAGAGTGTCAAAAGGGTGAGTATTTCGTCGACCGGGTGACGCTGTCCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAAAGAAGACTATACCCCGtat                                       <  2:353219/147‑1 (MQ=255)
                                                                                                        cTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAAAGAAGACTATACCCCGTATAc                                     >  1:14070/1‑151 (MQ=255)
                                                                                                            tgcccGTCCGGGACTGAAAAAAGAAGAGCTTTTTAAGGTTGATGATTTTGTTAACACGGCGACGCTGCCCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAAAGAAGACTATACCCCGTATACCt                                   <  2:66677/146‑1 (MQ=255)
                                                                                                                ccggccggggataaaaGTAGAAGAGCGTTTAAAGGGTGATGTTTTCGTTGCCCCGGTGACGCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAAAGAAGACTATACCCCGTATACCt                                   <  1:225717/133‑1 (MQ=255)
                                                                                                                     cggggagaaaaGTAGAAGAGCTTTTAAAGGGTGATGTTTTGTTTACCCCCGTGCCGCTGCCCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAAAGAAGACTATACCCCGTATACCTTCAc                               <  2:669845/137‑1 (MQ=255)
                                                                                                                             aaaGAAGAAGAGCTTTAAAAAGGTGATTATATCGTTGACCCGGTGACGCTGCCCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAAAGAAGACTATACCCCGTATACCTTCACCAAAGATCAGa                    <  1:74507/148‑1 (MQ=255)
                                                                                                                                            ttCAAAGGTGATGATATCGTTGACACCGTGACGCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAAAGAAGACTATACCCCGTATACCTTCACCAAAGATCAGATTGAAGAAGAGTTGCTGt  >  2:679190/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
AGGTATGGTACTGAATTACAACGGCAAACTGCTGCTGGTTAAGGATATTGATATTCAGTCGCCCACTGCCCGCGGTGCCGCTACGCTGTACAAAATGCGTTTTTCTGATGTCCGTACCGGGCTGAAAGTAGAAGAGCGTTTCAAAGGTGATGATATCGTTGACACCGTGACGCTGACCCGCCGTTACGTTGACTTCTCCTATGTCGATGGCAACGAATATGTCTTTATGGATAAAGAAGACTATACCCCGTATACCTTCACCAAAGATCAGATTGAAGAAGAGTTGCTGT  >  NC_000913/2265476‑2265765

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