Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,5810AC41.7% 87.8 / 11.8 48V314G (GTA→GGA) thiPthiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (25/3);  new base C (0/20);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.06e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTAC  >  NC_000913/73432‑73728
                                                                                                                                                      |                                                                                                                                                   
tCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTa                                                                                                                                                     <  2:429925/150‑1 (MQ=255)
tCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTa                                                                                                                                                     <  2:676445/150‑1 (MQ=255)
     cgcgCCGGCAGTTCGCGACTGCTGCATAGCAGCATGATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCAt                                                                                                                                               <  2:913941/151‑1 (MQ=255)
     cgcgCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCAt                                                                                                                                               >  1:340409/1‑151 (MQ=255)
       cgcCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATc                                                                                                                                              <  1:206197/150‑1 (MQ=255)
            gCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGAt                                                                                                                                        >  1:603943/1‑151 (MQ=255)
              aGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATc                                                                                                                                       >  1:446155/1‑150 (MQ=255)
                     gACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACATCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGcc                                                                                                                                 >  1:794151/1‑149 (MQ=255)
                            ccATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGCCGATCACCGCCAGTAACgg                                                                                                                         >  2:872349/1‑150 (MQ=255)
                             cATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTg                                                                                                                       >  1:701023/1‑151 (MQ=255)
                                       tcatGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAc                                                                                                               >  2:268051/1‑149 (MQ=255)
                                                gCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAgcgc                                                                                                    >  2:584311/1‑151 (MQ=255)
                                                 cACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGcc                                                                                                   >  1:654337/1‑151 (MQ=255)
                                                  aCTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGcc                                                                                                   >  1:627070/1‑150 (MQ=255)
                                                      cGCACAATCCCCCTGCCGCCAGCGAAATCCGCAACGAGGCCCACAGCGCTTGCCACAGCCCCGGTTGTGCCAGCCCTTCCGGCAA‑CTGGCGTTTTCCCCCAGCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAg                                                                                                 <  1:66866/148‑1 (MQ=255)
                                                      cGCACAATACCCCTGCCGCCAGCGCAATACGCAAGGAGGTCCACAGCGCCTGCCCCAGCCCCGGTTGGGCCAGCACTTCCGGCAA‑CGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAcagcagca                                                                                                        <  1:864280/141‑1 (MQ=255)
                                                        cacaATCCCCCGCCCCCCAGCGAAATACACAACGGGTTCCTCAGCGCTTGCCACAGCACCGGTAGGGCCAGCCCTTCCGCTAC‑CGGGGGATTTCCCCCACCGCCGACCGCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACaa                                                                                            <  1:632338/151‑1 (MQ=255)
                                                        caaaaTACCCCTCCCGCCAGCGCAATCCGCAACGAGGCCCACAGCGCCTGCCACAGCCCCGGTTGTGCCAGCCCTTCCGGCAA‑CGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAgcg                                                                                                     <  1:332587/139‑1 (MQ=255)
                                                           aatcccccTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCCCCTGCCACAGCCCCGGTTGCGCCAGCACTTCCGGAAA‑CGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGca                                                                                               <  1:239041/139‑1 (MQ=255)
                                                                 ccTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCg                                                                                   >  1:582970/1‑151 (MQ=255)
                                                                 ccTCCCGCCAGCGCAATACGCACCGAGGCCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCg                                                                                   <  1:200727/151‑1 (MQ=255)
                                                                     ccccccgcacactACGAAACGAGGCCCCCAGCCCCGGCCCCCGCACCGGTTGCGCCTCCTCTCCGGGCAA‑GTGGCGTTTTCCCCCATCCGCGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCgtgt                                                                                <  2:380396/138‑1 (MQ=255)
                                                                             gcAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGc                                                                       >  2:100378/1‑151 (MQ=255)
                                                                              aaatccgaagcGAGGCCCACGGCCCCTCCCACAGCCCCGTTCGAGCAAGCCCTCCCGGCAA‑CTGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAAt                                                                         <  1:234044/138‑1 (MQ=255)
                                                                               accccacAACGCGGCCCAAGGCGCCTCCCACACCCCCGGTTCCGCCTGCTCTCCCGACAA‑CTGGCATTTCCCCCCCTCGCCATTCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGc                                                                       <  1:899738/143‑1 (MQ=255)
                                                                                       aCGAGTCCCACGCCGCGTGCCACAGCCCCGGTTGTGCCACCACTCCCGGCAA‑CTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGc                                                               <  1:907602/149‑1 (MQ=255)
                                                                                       aCGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGc                                                               >  2:413922/1‑149 (MQ=255)
                                                                                         gAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGAc                                                           >  2:861301/1‑151 (MQ=255)
                                                                                         gAGGCCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CCGGCTATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGCGTCGCAAATGCGGCTATGCAGAc                                                           >  2:382634/1‑151 (MQ=255)
                                                                                                   gccccTGCCACAGCACCGGTGGTGCCAGCCCTCCCGGAAA‑CGGGCGTTTTCCCCCAGCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGAt                                                        <  2:503403/141‑1 (MQ=255)
                                                                                                      ccTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTc                                              >  1:814529/1‑151 (MQ=255)
                                                                                                               gCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGcc                                     >  1:382178/1‑151 (MQ=255)
                                                                                                                      ttGAGCCAGCAC‑TCCGGCAACCGGGCGATTGACCCCATCGCCAATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgcag                               <  2:326341/150‑1 (MQ=255)
                                                                                                                       tgtgCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgcagc                              >  2:473579/1‑150 (MQ=255)
                                                                                                                         ggCCACCACTCCCGGCAA‑CGGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgcagc                              <  1:905128/147‑1 (MQ=255)
                                                                                                                           ccAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCg                          >  2:131791/1‑150 (MQ=255)
                                                                                                                                cctcccGGCAG‑CGGGCGTTTTCCCCCCCCGGCACCCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCagcagc                           <  2:597934/140‑1 (MQ=255)
                                                                                                                                aCTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGc                    >  2:319286/1‑151 (MQ=255)
                                                                                                                                aCTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGc                    >  1:673418/1‑151 (MQ=255)
                                                                                                                                    ccGGCAA‑CGGGCGATTTCCCCCAGCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCg                  <  1:666229/149‑1 (MQ=255)
                                                                                                                                      ggCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgcg              >  1:523797/1‑151 (MQ=255)
                                                                                                                                       gccaccgggCGTTTTCCCCCGCCGCCGACCCCCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgc               <  2:290630/143‑1 (MQ=255)
                                                                                                                                        cAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgcgc             >  1:152760/1‑150 (MQ=255)
                                                                                                                                         aa‑CTGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCa            <  1:136405/150‑1 (MQ=255)
                                                                                                                                          accgggCGATTTCCCCCATCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCaa           <  2:753808/147‑1 (MQ=255)
                                                                                                                                          a‑CTGGCGTTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAAt          <  1:416768/151‑1 (MQ=255)
                                                                                                                                             tGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATgg        >  1:109664/1‑151 (MQ=255)
                                                                                                                                              ggcattttCCCCCATCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGc       <  2:899746/146‑1 (MQ=255)
                                                                                                                                                   tttCCCCCACCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCtt    <  1:41816/149‑1 (MQ=255)
                                                                                                                                                   tttACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTAc  >  1:390670/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAA‑CTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTAC  >  NC_000913/73432‑73728

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