Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,429,9840TG29.7% 62.9 / 10.2 37L41V (TTG→GTG) cycAserine/alanine/glycine/cycloserine:H(+)symporter
Reads supporting (aligned to +/- strand):  ref base T (4/22);  new base G (11/0);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 1.60e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAA  >  NC_000913/4429845‑4430117
                                                                                                                                           |                                                                                                                                      
ggTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt                                                                                                                                        <  2:28535/140‑1 (MQ=255)
ggTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt                                                                                                                                        <  2:159892/140‑1 (MQ=255)
ggTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt                                                                                                                                        <  1:251196/140‑1 (MQ=255)
  tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgttg                                                                                                                                      >  1:684027/1‑140 (MQ=255)
  tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgttg                                                                                                                                      >  2:680392/1‑140 (MQ=255)
     aGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGttt                                                                                                                                   >  1:294664/1‑140 (MQ=255)
      ggAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTa                                                                                                                                  <  1:680392/140‑1 (MQ=255)
       gaagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTAt                                                                                                                                 <  1:683569/140‑1 (MQ=255)
       gaagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTAt                                                                                                                                 >  1:66462/1‑140 (MQ=255)
       gaagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTAt                                                                                                                                 >  1:119901/1‑140 (MQ=255)
        aagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATg                                                                                                                                <  2:660619/140‑1 (MQ=255)
         agaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATgg                                                                                                                               >  2:157754/1‑140 (MQ=255)
          gaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATggg                                                                                                                              <  2:684027/140‑1 (MQ=255)
           aaaaaaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATgggg                                                                                                                             >  1:176795/1‑140 (MQ=255)
           aaaaaaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGCGCCACTTGTACGGGGGTGTTTCTgtgg                                                                                                                             >  2:395406/1‑140 (MQ=255)
                   aTGGTAGATCAGGTAAAAGTCGTTGCCGATGTTCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCaa                                                                                                                     <  1:55449/140‑1 (MQ=255)
                     ggTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCaaaa                                                                                                                   >  2:461791/1‑140 (MQ=255)
                        aGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGa                                                                                                                <  1:270320/140‑1 (MQ=255)
                             aGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTTGGTTGGGGGTGTTTTTGGGGGCTGGCCAAAAGaattac                                                                                                           >  2:185108/1‑135 (MQ=255)
                               gTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGTACGGGGGTGTTTATTGGGGCTGGCAAAACGATTAACCt                                                                                                        >  2:541128/1‑140 (MQ=255)
                                      tCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCggg                                                                                                 <  1:514944/140‑1 (MQ=255)
                                         ttGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCg                                                                                              <  1:734188/140‑1 (MQ=255)
                                                tgatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTTGTACGGGGGTGTTTTTGGGGGCTGGCAAAAAGATTTGCCTTGCCCGGGCGGCGGTcc                                                                                       >  2:399065/1‑139 (MQ=255)
                                                 gatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAt                                                                                      <  2:338456/140‑1 (MQ=255)
                                                  atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAtt                                                                                     <  2:176795/140‑1 (MQ=255)
                                                   tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAttt                                                                                    <  1:541128/140‑1 (MQ=255)
                                                      ggCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCg                                                                                 <  2:66462/140‑1 (MQ=255)
                                                      ggCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCg                                                                                 <  2:663711/140‑1 (MQ=255)
                                                       gCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGt                                                                                <  1:157754/140‑1 (MQ=255)
                                                           cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTtat                                                                            <  1:495360/140‑1 (MQ=255)
                                                           cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTtat                                                                            >  2:767771/1‑140 (MQ=255)
                                                             gCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTatat                                                                          <  1:780385/140‑1 (MQ=255)
                                                             gCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTatat                                                                          >  1:700236/1‑140 (MQ=255)
                                                             gCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTatat                                                                          >  2:42302/1‑140 (MQ=255)
                                                                 aaCAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATc                                                                      >  1:220985/1‑140 (MQ=255)
                                                                     gTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGGCCATGGGAAGGGGGTGTTTTTTGGGGCTTGCAAAAAGAGTAGCCCTTGCGGGGCGGTGATCCTTTTTGTTTTTATTATCAttt                                                                  >  2:279390/1‑139 (MQ=255)
                                                                          tACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGtttt                                                             <  1:395406/140‑1 (MQ=255)
                                                                             ggCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATg                                                          <  2:700236/140‑1 (MQ=255)
                                                                                    tctcACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGTACGGGGGTGTTTTTGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCtttt                                                   >  2:687513/1‑140 (MQ=255)
                                                                                     ctcACAAACCGACATATTCAGCTTATTGCCATTGGCGGGGCCCTTGGGACGGGGGTGTTTATGGGGTCTGTCAAAAAGAGTAGCCCTGCCGGGCCGTCGAACACTTTTGTTTATATGATCATTGGTTTTTTGCTCttttt                                                  >  2:477984/1‑140 (MQ=255)
                                                                                                 cATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGc                                      <  1:461791/140‑1 (MQ=255)
                                                                                                    atTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAAt                                   <  2:119901/140‑1 (MQ=255)
                                                                                                        aGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTTTTTTTGGGGGCTGGCAAAAAGATTAGCCTTGCCGGGCCGGCGGTCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGGTGCGGGCAATgggg                               >  2:280991/1‑140 (MQ=255)
                                                                                                        aGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTTTGCTCTTTTTCGTGATGCGGGCAATgggg                               >  2:712142/1‑140 (MQ=255)
                                                                                                             aTTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTGGGGTCTGGCAAAACGATTAGCCTTGCCGGGGCGTCGATCATTTTTGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGGAtt                          >  2:538498/1‑140 (MQ=255)
                                                                                                              ttGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATtg                         <  1:423401/140‑1 (MQ=255)
                                                                                                              ttGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATtg                         <  1:295256/140‑1 (MQ=255)
                                                                                                                  cATTGGCGGTGCCATTGGTACGGGGGTGTTTATTGGGGCTGGCAAAAAGATTAGCCTTTCCGGGCCGTCGGTCATTTTTGTTTATATTATCATTGGTTTTTTGCTCTTTTTTGTTATGCGGGCAATGGGGGGATttctgc                     >  2:295256/1‑140 (MQ=255)
                                                                                                                  cATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGGCGGCGATCATTTTCGTTTATATGATCATTTGGTTTATGCTCTTTTTTGTGATGCGGGCAATGGGGGGATtgctgt                     >  2:423401/1‑139 (MQ=255)
                                                                                                                                     aCGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACaa  >  2:267000/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                      
GGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAA  >  NC_000913/4429845‑4430117

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