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

GCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACT  >  NC_000913/73435‑73729
                                                                                                                                                  |                                                                                                                                                    
gCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAcccc                                                                                                                                                  >  1:650031/1‑151 (MQ=255)
   gcgcCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCAt                                                                                                                                                <  2:684780/150‑1 (MQ=255)
    cgcCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATcg                                                                                                                                              >  1:343400/1‑151 (MQ=255)
     gcCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATcg                                                                                                                                              >  2:66099/1‑150 (MQ=255)
      cccGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATcgac                                                                                                                                            <  2:407063/151‑1 (MQ=255)
         gCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGAt                                                                                                                                         >  1:532279/1‑151 (MQ=255)
          cAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACCCCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGAt                                                                                                                                         <  2:456254/150‑1 (MQ=255)
           aGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATc                                                                                                                                        >  1:658072/1‑150 (MQ=255)
                       ccccAAAGCAGCATCAGGGTCAGCACTCCCCAAAATCCCCCGGCCGCCGGCGAAATACGCAACGAGGCCCACAGCGCCTGCCACAGCACGGGTTGTGCCAACCCTCCCGGCAACGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACg                                                                                                                           <  1:108176/149‑1 (MQ=255)
                            tAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGgc                                                                                                                      >  1:760237/1‑151 (MQ=255)
                                  catcatGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAc                                                                                                                >  2:360034/1‑151 (MQ=255)
                                       gggTTAGCACTCCCCAAAACACCCCTCCCCCCGGCGAAATACATAGCGAGGCCCACGGCGCCCGCCACAGCCCCGTTTGGGCAAGCACTTCGGGCAACTGGCTATTTCCCCCACCGATGATCTCCGCCAGTAACGGTGGCAGCAAcag                                                                                                              <  1:369824/147‑1 (MQ=255)
                                         gTCAGCACCCCGCAAAATCCCCCGGCCGCCAGCGCAATCCGCAACGGGGTCCACGGCCCCGGCCACAGCCCGGGTTGGGCCAGCCCTCCCGGCAACGGGCATTTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAAcagcagc                                                                                                          <  1:472110/150‑1 (MQ=255)
                                                        aaTACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAAtt                                                                                           >  2:152743/1‑150 (MQ=255)
                                                         aTACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTa                                                                                          >  1:412848/1‑150 (MQ=255)
                                                           acacCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAac                                                                                        >  2:1566/1‑150 (MQ=255)
                                                            ccccTGCCGCCGGCGAAATCCGCAAGGAGGCCCACAGCGCCTGCCACAGCCCCGTTTGTGCCAGCACTTCCGGAAACGGGCGATTTCCCCCATCGCCGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAacac                                                                                      <  1:638263/149‑1 (MQ=255)
                                                             cccTGCCGCCAGCGCAATCCGCAACGAGGCCCACAGCCCCTCCCACAGCCCCGGTTGTGCCAGCCCTTCCGGCAACGGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACAcc                                                                                     <  1:313549/150‑1 (MQ=255)
                                                               cTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCgt                                                                                   >  1:680034/1‑151 (MQ=255)
                                                                tGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCgtg                                                                                  >  1:760414/1‑151 (MQ=255)
                                                                            aaTACGAAACGAGGTCCACAGCCCCTGCCACAGCACCGTTTGTGCCAGCCCTTCCGGCAACGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCgt                                                                                   <  1:414797/138‑1 (MQ=255)
                                                                                      ggggcccACACCCCCTGCCACAGCCCCGTTTGCCCCACCCCTCCCGGCAACTGGCGATTTCCCCCATCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCt                                                                    <  1:251105/138‑1 (MQ=255)
                                                                                        gtcccACAGCGCCTGCCACACCCCCGTTTGTGCCAGCACTTCCGGCAACTGGCGTTTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGa                                                          <  1:606871/148‑1 (MQ=255)
                                                                                           ccccgccgcCTCCCCCACCCCCGGTTGGGCCAGCACTCCGGGCAACGGGCTATTCCCCCAACCAACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGAt                                                         <  2:71157/143‑1 (MQ=255)
                                                                                                 cgcCTGCCACAGCCCCGGTTGTGCCAGCTCTTCCGGCAACTGGCGTTTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCgg                                                  <  1:91958/150‑1 (MQ=255)
                                                                                                   ccggCCACCGCAGCGGGCGGGCCATCTCTTCCGACGACGGCAGATTTCCCCCCGCGCAAATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATc                                                        <  2:99909/139‑1 (MQ=255)
                                                                                                   ccTGCCACAGCCCCGGTTGTGCCAGCCCTTCCGGCAACGGGCGATTTCCCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGt                                                <  1:585863/150‑1 (MQ=255)
                                                                                                     tGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACTCCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTcgc                                             >  2:406471/1‑151 (MQ=255)
                                                                                                      gCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTcgcg                                            >  1:498368/1‑151 (MQ=255)
                                                                                                        cacaGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGcc                                          >  2:307724/1‑151 (MQ=255)
                                                                                                                     gggCCAGCACTTCCGGCAACTGGCGATTTCCCCCAGCGCCGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTgcagca                              <  2:201270/148‑1 (MQ=255)
                                                                                                                      tgCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCagcagc                            >  2:464486/1‑151 (MQ=255)
                                                                                                                       gCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCg                           >  2:587448/1‑151 (MQ=255)
                                                                                                                           gCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCgtggt                       >  1:63944/1‑151 (MQ=255)
                                                                                                                                 ccGGCAACGGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCgg                  <  1:651307/150‑1 (MQ=255)
                                                                                                                                   ggtaacgggctctttCCCCCACCACCGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgc                <  2:767307/138‑1 (MQ=255)
                                                                                                                                   ggCAACGGGCGATTCCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgc                <  2:294832/150‑1 (MQ=255)
                                                                                                                                   gacAACTGGCAATTTCCCCCATCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCg                   <  1:744136/145‑1 (MQ=255)
                                                                                                                                      aaCTGGCGATTTCCCCCACCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGgcgc              <  1:360035/149‑1 (MQ=255)
                                                                                                                                          ggCGTTTTCCCCCAGCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAAt           <  2:711486/148‑1 (MQ=255)
                                                                                                                                          ggCGATTCCCCCCCCCGACGACCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCa             <  2:271586/146‑1 (MQ=255)
                                                                                                                                              aTTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCt      >  1:641079/1‑149 (MQ=255)
                                                                                                                                                ttACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACt  >  1:592420/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACT  >  NC_000913/73435‑73729

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: ‑