Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913827,5070TG50.0% 54.9 / 13.9 46D492A (GAC→GCC) ybhFputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (3/20);  new base G (23/0);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 6.32e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGG  >  NC_000913/827362‑827657
                                                                                                                                                 |                                                                                                                                                        
gCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACg                                                                                                                                                     >  1:51337/1‑151 (MQ=255)
   ggCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGcc                                                                                                                                                   <  1:192010/150‑1 (MQ=255)
   ggCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCg                                                                                                                                                  <  1:95364/151‑1 (MQ=255)
      gATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGcccg                                                                                                                                                 >  2:526265/1‑147 (MQ=255)
       atcaaTTTCCCGCGGTTCACCAGGCCCATGAGGTCGCCATATTCCGCTTCATCCATAAAGCGGGTGGTGACTATCACCGTGACGCGTTTCTCTTCCATGCTGTTGATGTACAGCCAAAATTCGCGGCGCGTGAGGGGGtgat                                                                                                                                                       <  2:605912/139‑4 (MQ=255)
        ttAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGTCAACGcc                                                                                                                                                   >  2:504345/1‑145 (MQ=255)
          aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTc                                                                                                                                            <  1:891036/150‑1 (MQ=255)
             ttCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCgg                                                                                                                                          >  2:412120/1‑149 (MQ=255)
                ccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGcccgggatcgg                                                                                                                                          >  2:419067/1‑137 (MQ=255)
                 cgcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTtcg                                                                                                                                      >  1:313749/1‑149 (MQ=255)
                       acacCGGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACTCCGGAAGTCGGTTCGGCGAGgaa                                                                                                                              <  2:34331/151‑1 (MQ=255)
                       acacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGGGGGTTTGCCGAGaa                                                                                                                               >  1:723593/1‑150 (MQ=255)
                                cGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                                                                                    <  1:450099/151‑1 (MQ=255)
                                cGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGGCGAGAAAAAGAATGt                                                                                                                     >  1:350661/1‑150 (MQ=255)
                                         cGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATg                                                                                                           <  1:239281/151‑1 (MQ=255)
                                          gCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTCCgggt                                                                                                               >  1:360590/1‑144 (MQ=255)
                                             atatTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGGCGGTTTCTCGAGAAAAAGAATGGTCGGGTCATGCa                                                                                                         >  2:736340/1‑149 (MQ=255)
                                                            cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagc                                                                                        <  1:137608/151‑1 (MQ=255)
                                                            cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagc                                                                                        <  1:620747/151‑1 (MQ=255)
                                                            cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagc                                                                                        <  1:520893/151‑1 (MQ=255)
                                                              tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTTGTCGAGAAACAAAAAGTCCGGGTCCTGCATCAGCGAAAAAgcc                                                                                           >  1:625993/1‑146 (MQ=255)
                                                               aaaGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAAAAAgcccgcgcc                                                                                     >  1:70561/1‑151 (MQ=255)
                                                                aaGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgcca                                                                                    <  1:490569/151‑1 (MQ=255)
                                                                  gtgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAAACGGGTCCTCGAGAAAAAAAAAGTCCGGTTCATGCATCAGCGAACAAgcccatgcc                                                                                     >  2:733919/1‑142 (MQ=255)
                                                                      gtggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGATGGTTTGTCGAGAAAAAGAATGTCCGGGTCAAGCATAAGCGAACAGgccagcgcccacc                                                                                 >  1:816567/1‑144 (MQ=255)
                                                                              cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGCGAAACAAAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGGTTtaa                                                                      >  2:683048/1‑150 (MQ=255)
                                                                              cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGGAAAGGCGGGAGTCGGTTCTTCGCGAAAAAGAATATCCGGGTTATGGATCAGAGAACAGGGCAGCGCCAGCCGCTGtttta                                                                       >  2:334997/1‑147 (MQ=255)
                                                                                 cACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCt                                                                  <  2:838993/151‑1 (MQ=255)
                                                                                  aCCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAcc                                                                   <  2:405478/149‑1 (MQ=255)
                                                                                      tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCCAGAAACAAAAAGGCCGGGTCATGCATCAGCGGACAAGCCCGCGCCAGCCGCTGTTTTAAACCTAAttg                                                             >  2:502744/1‑149 (MQ=255)
                                                                                       gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                            <  2:563406/151‑1 (MQ=255)
                                                                                          gCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGt                                                         <  2:413433/151‑1 (MQ=255)
                                                                                                  ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATc                                                    <  1:946311/148‑1 (MQ=255)
                                                                                                  ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCCGCCGCTGTTTAAAACCTAATGGCAGGTCATCgg                                                  >  1:569741/1‑150 (MQ=255)
                                                                                                   tACCATGCTGTTGATGTGCAGCCAAAATTCACGGGGGATGGGGGGGGCAAAACCGGGAGGCGGTTCGTCGAGTAAAAGAAATGGGCGGTCCTGCAACCGCGAGCAAGGCAGCGGCAACACCCGGTTTTAAACCTATGGCAGTTTATCgg                                                  >  2:914192/1‑149 (MQ=255)
                                                                                                   tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGGAGGCGGGTTGTCGAGAAAAAGAAAGTCCGGGTCATGGATCCGCGGACAGGCCCGCGCCAGCCGCTGTTTTAAAACTAATTGCAGGTCATCGGTg                                                >  2:138668/1‑151 (MQ=255)
                                                                                                   tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAAAAGAATGTACGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTTAAACCTAATGGCAGTTCATCGGTg                                                >  2:719960/1‑151 (MQ=255)
                                                                                                        tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGCGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGt                                            <  2:413356/150‑1 (MQ=255)
                                                                                                         gCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTg                                           <  2:842773/150‑1 (MQ=255)
                                                                                                                 tgtgCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCCGTTTTAAACCTAATGGCAGTTCATCCGTGGGGGGGGAgggg                                    >  1:502693/1‑147 (MQ=255)
                                                                                                                     cAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACt                               <  2:568493/150‑1 (MQ=255)
                                                                                                                       gCCAAAATTCACGGCGGGTGGGGGGGGCAACCCCGGGAGTCGGGTCGTCGTGAAACAGAATAGCCGGGTCATGCTTCCGCGAACAAGCCCGCGCCCGCCGCTGTTTTAAAAATAATGGGGGTTCTTCGGGGGGGGGGGGGGGGGGACtttt                            >  2:770648/1‑151 (MQ=255)
                                                                                                                            aaTTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCCGGTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAAAAGGCCAGCGCCAGCCGCTGTTTTAAAACCAATGGCAGTTCATCGGGGGCGTGGGGGGGGGTACTTTTTAgg                        >  2:192286/1‑150 (MQ=255)
                                                                                                                                 aCGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAAc                  <  1:531047/151‑1 (MQ=255)
                                                                                                                                     cggGTTAGGGGGTTGACGCCGGCTGGCGGTTCGTAGAGAAAGATTAAGTCCGGTTCAGAGATCAGCGAACTGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCATTTCATCGGGCGCGTGGGAGGCGGTTCTTTTCAGGCCGAATGc                <  2:791297/149‑1 (MQ=255)
                                                                                                                                            gggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCAt       >  2:263274/1‑151 (MQ=255)
                                                                                                                                            gggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCAt       >  1:908373/1‑151 (MQ=255)
                                                                                                                                                 tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCggg  >  1:257082/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                        
GCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCGGG  >  NC_000913/827362‑827657

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: