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

GTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATAC  >  NC_000913/4429846‑4430115
                                                                                                                                          |                                                                                                                                   
gTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgtt                                                                                                                                    <  2:147213/140‑1 (MQ=255)
 tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGGCATTGGTACGGGggtg                                                                                                                                   >  1:98497/1‑140 (MQ=255)
 tACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgttg                                                                                                                                   >  2:466764/1‑140 (MQ=255)
     ggAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTa                                                                                                                               <  1:129182/140‑1 (MQ=255)
     ggAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTa                                                                                                                               <  1:171501/140‑1 (MQ=255)
     ggAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTa                                                                                                                               <  1:748729/140‑1 (MQ=255)
      gaagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGTACGGGGGTGttttt                                                                                                                              >  2:605967/1‑138 (MQ=255)
      gaagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTAt                                                                                                                              >  1:411406/1‑140 (MQ=255)
       aagaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATg                                                                                                                             <  2:225589/140‑1 (MQ=255)
          aaaaaaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTtggg                                                                                                                          >  2:639246/1‑140 (MQ=255)
          aaaaaaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATgggg                                                                                                                          >  1:672078/1‑140 (MQ=255)
                aCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGc                                                                                                                    >  2:6264/1‑140 (MQ=255)
                  aTGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCaa                                                                                                                  <  2:68409/140‑1 (MQ=255)
                        gATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGAt                                                                                                            <  1:222808/140‑1 (MQ=255)
                                           ccGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCg                                                                                         <  2:116025/140‑1 (MQ=255)
                                                gatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGGACGGGGGTGTTTTTGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGGTCAt                                                                                    >  2:224636/1‑140 (MQ=255)
                                                gatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAt                                                                                    <  1:165765/140‑1 (MQ=255)
                                                 atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGACAAAACGATTAGCCTTGCCGGGCCGTCGATCAtt                                                                                   <  1:653867/140‑1 (MQ=255)
                                                   cAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAtttt                                                                                 <  1:575648/140‑1 (MQ=255)
                                                     ggCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCg                                                                               <  2:280819/140‑1 (MQ=255)
                                                     ggCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCg                                                                               <  1:629248/140‑1 (MQ=255)
                                                       cTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGtt                                                                             <  1:217555/140‑1 (MQ=255)
                                                         ccGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTta                                                                           <  2:272149/140‑1 (MQ=255)
                                                          cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGGACGGGGGTGTTTTTGGGGGCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTtat                                                                          >  2:296132/1‑140 (MQ=255)
                                                          cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTtat                                                                          <  1:316709/140‑1 (MQ=255)
                                                           ggCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGTACGGGGGTGTTTTTGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTTGTTtata                                                                         >  1:437002/1‑140 (MQ=255)
                                                            gCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTatat                                                                        <  2:645987/140‑1 (MQ=255)
                                                            gCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTatat                                                                        <  1:466764/140‑1 (MQ=255)
                                                             cTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATg                                                                       <  2:357297/140‑1 (MQ=255)
                                                                       gCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTTGGACGGGGGTGTTTTTTGGGGCTGGCCAAAAGATTAGCCTTGCCGGGGCGTCGGTCCTTTTTGTTTTTATGGTCAtttggt                                                             >  2:250595/1‑136 (MQ=255)
                                                                       gCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTTTGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTTGTTTATATGATCATTGGtt                                                             >  1:122041/1‑140 (MQ=255)
                                                                              cgcAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGct                                                      <  2:42745/140‑1 (MQ=255)
                                                                               gcAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTTGTACGGGGGTGTTTTTGGGGGCTGGCAAAAAGATTAGCCTTGCCGGGGCGTCGATCATTTTTGTTTTTATGATCATTTGTTTTTTGctc                                                     >  2:671591/1‑140 (MQ=255)
                                                                                             cGACATATTCAGCTTATTGCCATTGGCGGTGCCCTTGGTACGGGGGTGTTTTTTGGGCCTTGCAAAACGATTAGCCTTGCCGCGGCGTCGATCATTTTTTTTTTTATTATCATTTGTGTTTTGCTCCTTTTCCTGATgtg                                       >  2:549113/1‑138 (MQ=255)
                                                                                             cGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCg                                       <  1:6264/140‑1 (MQ=255)
                                                                                             cGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCg                                       >  1:266588/1‑140 (MQ=255)
                                                                                               aCATATTCAGCTTATTGCCATTGGCGGTGGCATTGGTACGGGGGTGTTTTTGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTTTGCTCTTTTTCGTGATGCggg                                     >  2:46613/1‑140 (MQ=255)
                                                                                               aCATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCggg                                     <  1:662502/140‑1 (MQ=255)
                                                                                                         cTTATTGCCATTGGCGGTGCCCTTTGTACGGGGTTGTTTTTGGGGGCTGGCCAAACGGTTAGGCTTGCCGGGGCCGCGATCATTTTTGTTTATATTATTATTGGGTTTTTTCTCTTTTTTGTGTTGCGGGGAATtggggg                           >  2:720329/1‑139 (MQ=255)
                                                                                                           tATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTCGCAAAACGAGTAGCCTTGCCGGGGCGTCGGTCATTTTCGGTTATATGAATATTGGTTTTATTCTCTTTTTCGTGATGCGGGCGACggggggaa                         >  2:302744/1‑137 (MQ=255)
                                                                                                            aTTGCCATTGGCGGTGCCCTTGGTACGGGGGTGTTTTTGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGGTCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAAtt                        >  1:325261/1‑140 (MQ=255)
                                                                                                              tGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAGCCTTGCCGGGGCGGCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATtgc                      >  1:356084/1‑140 (MQ=255)
                                                                                                               gCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATtgct                     <  2:301482/140‑1 (MQ=255)
                                                                                                                   ttGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGGCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATACTCTTTTTCGTGATGCGGGCAATGGGGGGATTGCTGCtt                 >  2:244514/1‑140 (MQ=255)
                                                                                                                      gCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCg              <  2:672078/140‑1 (MQ=255)
                                                                                                                         gTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAAt           <  1:302744/140‑1 (MQ=255)
                                                                                                                             ccttGGGACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAGCCTTGCCGGGGCGTCGATCATTTTCGTTTATATGATCATTGGTTTTTTGCTCTTTTTTGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGCATCTgg       >  2:17163/3‑140 (MQ=255)
                                                                                                                              aTTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGa      <  2:193205/140‑1 (MQ=255)
                                                                                                                                  gTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATAc  <  2:24434/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                   
GTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATAC  >  NC_000913/4429846‑4430115

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