Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913827,5070TG56.4% 29.4 / 50.0 39D492A (GAC→GCC) ybhFputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (3/14);  new base G (22/0);  total (25/14)
Fisher's exact test for biased strand distribution p-value = 4.51e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-05

TGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTT  >  NC_000913/827364‑827631
                                                                                                                                               |                                                                                                                             
tGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGAGGGGGTCAACgt                                                                                                                         >  2:325742/1‑149 (MQ=255)
 ggCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGcc                                                                                                                        <  2:175974/150‑1 (MQ=255)
  gCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGTGaaacgcgg                                                                                                                      >  1:741794/1‑145 (MQ=255)
       tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGt                                                                                                                  <  1:757373/150‑1 (MQ=255)
        aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCg                                                                                                                <  2:4632/151‑1 (MQ=255)
         aTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGTGGGGGGCAACGCCCGAAGGCgg                                                                                                               >  1:592303/1‑151 (MQ=255)
          tttCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAggagg                                                                                                               >  1:203764/1‑146 (MQ=255)
               cgcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCgggt                                                                                                             >  2:79871/1‑145 (MQ=255)
                gcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgtcg                                                                                                        <  2:573067/151‑1 (MQ=255)
                 cgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGTCACGCCGGGAAGCGGgtcgtc                                                                                                         >  1:509859/1‑149 (MQ=255)
                   gTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCCGAAGTCGGGtcgt                                                                                                          >  2:331309/1‑146 (MQ=255)
                      cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCCGAAGGCGGTTTGGCGAGgaa                                                                                                   >  1:709378/1‑150 (MQ=255)
                        ccAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGACGGGTTGTCGAGgaa                                                                                                   >  2:688915/1‑148 (MQ=255)
                          aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAAAAGa                                                                                               >  1:736688/1‑150 (MQ=255)
                          aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAAAGCCGGAAGGCGGGTCGTCGAGGAACa                                                                                                 >  1:727573/1‑148 (MQ=255)
                            gCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGAGCAGCCAAAATTCACGGCGGGAGAGGGGGGCAACGCCGCGAGACGGTTCGCCGCGAAACAg                                                                                                >  2:656690/1‑147 (MQ=255)
                                     gTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCa                                                                                   <  1:5527/151‑1 (MQ=255)
                                      tCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGGCGGTTCGTCGAGAAACAGAAAGTCCGGTTCAt                                                                                  >  2:343536/1‑151 (MQ=255)
                                                ccGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGAAGTCGGGTTGTCGAGAAAAAGAATGTCCGGGTCCTGCATCAg                                                                           >  1:742759/1‑148 (MQ=255)
                                                ccGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGTGAAAAAAAATGTCCGGGTCATGGCTCAGc                                                                          >  1:743821/1‑149 (MQ=255)
                                                       aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcc                                                                 <  2:942061/151‑1 (MQ=255)
                                                        tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAAGCGGTTCGTCGAGAAAAAGAAAGTCCGGTTCATGCATCAGCGAACAAgcc                                                                 >  2:369353/1‑150 (MQ=255)
                                                               aGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgccag                                                         <  1:465590/151‑1 (MQ=255)
                                                                      ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTg                                                   <  1:300448/150‑1 (MQ=255)
                                                                        tgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGttt                                               >  1:528637/1‑151 (MQ=255)
                                                                             aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAc                                          <  2:815114/151‑1 (MQ=255)
                                                                               cACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCCTCCGCGAACAGGCCCGCGCCCGCCGCCGTTTtaaaac                                         >  1:530050/1‑148 (MQ=255)
                                                                               cACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGTTCGTCGAGAAACAGAATGTCCGGTTCATTCAACAGCGAACAGGCCAACGCCCGCCGCCGTTTTAAAcc                                         >  2:142976/1‑150 (MQ=255)
                                                                                 ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAACCCTa                                       <  1:770295/150‑1 (MQ=255)
                                                                                   gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATg                                    <  1:645698/151‑1 (MQ=255)
                                                                                      aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCa                                 <  2:239342/151‑1 (MQ=255)
                                                                                      aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAATCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGAAAAAGCCAGCGCCAGCCCCCGGTTTAAACCTAATgg                                   >  2:411212/1‑149 (MQ=255)
                                                                                                  aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTg                      <  1:303215/150‑1 (MQ=255)
                                                                                                    cATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTTCGGGTCATGCATCAACCAAAAGGCCAGCGCCAGCCGCTGTTTAAAAACTAATGGCAGGTCATCGGTgg                     >  2:729220/1‑149 (MQ=255)
                                                                                                     aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCCGTTTTAAACCCAAAGGCAGTTCATCGGGGGCg                   >  1:631035/1‑150 (MQ=255)
                                                                                                                 tgCAGCCAAAATTCACGGCGGGAGAGGGGGGGAACGGCGGAAGTCGGGTCGGCGAGAAACAGAATGTCCGGGTCATGCATGAACGGACAAGCCAGGGACAACCGCTGTTTAAAAACCAATAGCAGGTCCACGGGGGCgggggggggggta       >  2:825910/1‑145 (MQ=255)
                                                                                                                   cAGCCAAAATTCACGGCGGGTGAGGGGGGCCACGCCGGGAGACGGGTCGTCGAGAAAAAGAATGTCCGGGTCCTGCATCAGCGAACAGGCCAACGCCAGCCGCCGTTTTAAACCTAATTGCAGGTTATCGGGGGGGGGGGGGGGGGTACtt    >  1:489206/1‑151 (MQ=255)
                                                                                                                     gCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATTCtttt  <  1:666501/151‑1 (MQ=255)
                                                                                                                      ccAAAATTCACGGCGGGGGAGGGGGGGAAAGCCCGAAATCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCCGTTTTAAAACCAAATGCAGGTCATCGGGGGGGGGGGGGGGGGTACtttt  >  1:935032/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                             
TGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTT  >  NC_000913/827364‑827631

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