Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913827,5070TG45.1% 61.7 / 11.2 51D492A (GAC→GCC) ybhFputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (9/19);  new base G (23/0);  total (32/19)
Fisher's exact test for biased strand distribution p-value = 1.62e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCG  >  NC_000913/827369‑827655
                                                                                                                                          |                                                                                                                                                     
aTTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAg                                                                                                                                           <  1:1012848/151‑1 (MQ=255)
  tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTc                                                                                                                                         <  1:907618/151‑1 (MQ=255)
              gTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGCCCGGCGTCGGTttgtcccga                                                                                                                             >  2:317726/1‑147 (MQ=255)
                acacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGGGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCCACGCCGGGAGTggggttttccggaaa                                                                                                                           >  1:169091/1‑139 (MQ=255)
                 cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGGCGAGaaaa                                                                                                                          >  1:78974/1‑150 (MQ=255)
                 cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGCGaaaa                                                                                                                          >  1:357343/1‑150 (MQ=255)
                                    cAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCa                                                                                                       <  2:1148356/151‑1 (MQ=255)
                                      atatTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTGAACGGCGGGAGGCGGGTTGTCGAGACAAAGAAAATCCGGtgcatgcac                                                                                                      >  2:740697/1‑149 (MQ=255)
                                      atatTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCCCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCAt                                                                                                      >  1:1099745/1‑150 (MQ=255)
                                        tttCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGTGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                                                                                                   <  2:1120046/150‑1 (MQ=255)
                                                   tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCCGAAGTCGGTTCGGCGGGAAAAAGAATGGCCGGTTCATGCATCAGCGAAAAgggc                                                                                        >  1:242523/1‑148 (MQ=255)
                                                    ccATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAACAGAAAGTCCGGGTTATGCATCAGCGGAAAGgcccg                                                                                      >  1:310708/1‑149 (MQ=255)
                                                       tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgc                                                                                   <  2:824295/151‑1 (MQ=255)
                                                       tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgc                                                                                   <  1:343253/151‑1 (MQ=255)
                                                       tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATGCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgc                                                                                   <  2:564184/151‑1 (MQ=255)
                                                       tAAAGAGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTAAACGCCCGGAGTCCGTTTGTCGAGGAAAAGAATGTTCGGGTCTTGGATTAGCGAAAAAgcccgcg                                                                                    >  1:575037/1‑150 (MQ=255)
                                                                   tgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGGGCCAGCCGCTGttt                                                                       <  1:932720/151‑1 (MQ=255)
                                                                    gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTCCGGGTCCTGCATCAGCGAACAGGCCAGCGCCCGCCCCCGTTTa                                                                      >  1:244453/1‑151 (MQ=255)
                                                                    gACCATCACCGTGACGCCTTCCTCTACCATGCTGTTGATTTGCAGCCAAAATTCACGGCGGGCGGGGGGGTTCATGCCCGGACTCGCGTCGTTTGGGAACACCAAGTGCGGTTCAATCATCAGCGCACAAGCCAGCGCCAGcccaggtttt                                                                      >  2:687777/1‑143 (MQ=255)
                                                                        aacACCGTGACGCCTGTCTCTACCATGCGGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATACATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAc                                                                  <  2:61697/149‑1 (MQ=255)
                                                                           aCCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTGCGGGTCATGCCTCAGCGAAAAGGCCAGCGCCCGCCCCCGTTTTAAACCTa                                                               >  2:926446/1‑151 (MQ=255)
                                                                            ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGTGGGGGGCCACGCCCGGAGGCGGTTCGTCGGGAAACAGAAAGTCCGGTTCATGGATCAGCGAACAAGCCCGCGCCAACCGCCGTCTTTAACCtta                                                              >  2:253589/1‑149 (MQ=255)
                                                                                 aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                          <  1:764275/150‑1 (MQ=255)
                                                                                 aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                          <  1:838906/150‑1 (MQ=255)
                                                                                    ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGTGGAGTCGTCCAGAAACAGAAAGTTCGGGTTATGCATCAACGAACAGGCCAACGGCCGACGCTGTTTAAAACCTAATGTCCGtt                                                      >  2:818616/1‑151 (MQ=255)
                                                                                    ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGTCGGGTCGTCGGGAAAAAGAATGTCCGGGTCATGCCTCAGCGAACAGGCCCGCGCCCGCCGCCGTTTTAAAACTAATGGGAggt                                                      >  1:371503/1‑149 (MQ=255)
                                                                                    ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGGAGTCGGTTCGTCGAGAAACAGAAAGTCCGGGTCATTCATCAACGAAAAAGCCAGCGCCCGGCCCTGTTTTAAAACCAATTGCAggt                                                      >  1:121796/1‑149 (MQ=255)
                                                                                      tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCTGGAGTCGGTTCGTCGAGGAAAAGAATGTCCCGGTCCTGGCTCCGCGAAAAGGCCCGCGCCCGCCGCCGTTTAAAACCCAAACGCAggtta                                                    >  2:75150/1‑147 (MQ=255)
                                                                                      tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTTGGGTTTTCGAGAAAAAGAATGTCCGGGTCATGTCTCCGCGAACAAGCCAGCGCCAGCCGCTGTTTTAAAACTAATGGccggt                                                      >  1:540062/1‑145 (MQ=255)
                                                                                      tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGGGGGGGGCGACGCCCGAAGCCCGTTCTTCCGGAAAAAGAGTGGCGGGTGCCTCCATCTGCGAACAGGCCCGCGCCCCCCGCCGTTCAAAACCTAATGCCAGtgtt                                                    >  2:82671/1‑148 (MQ=255)
                                                                                        tctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGTAACGCCCGAAGGGGGTTTCGCTAGAAAAAGAGTGTTCGGGTCATTTGACAGCGGAAAGGACCGCGCGAGCCCGCGTTTTTAACACAATGGCAgcttatg                                                  >  2:766095/1‑145 (MQ=255)
                                                                                        tctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGGCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAAAAAGCCAGCGCCAGCCGCCGTTTAAAAACTAATGGCAggttat                                                   >  1:198089/1‑145 (MQ=255)
                                                                                           ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGTCAACGCCCGAAGTCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCAACAGCGAAAAGGCCCGCGCCCGCCCCTGGTTAAAAACTAAAGGGAGGTCATCggg                                               >  1:1066461/1‑150 (MQ=255)
                                                                                           ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAAGGGGGCAACTCGGGCCCACGGTTCCGTGGGAAACAGAAAGTTCGGTTAAAGCACCAGCCAACACGCCAGGCCCCACCGCTTTTTTAAACCCAATCGGAGTTTAtaagt                                               >  1:822920/1‑147 (MQ=255)
                                                                                             aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTTGTCGGGAAAAAGAATGTCCGGTTCATGCCTCAGCGAACAGGCCCGCGCCCCCCCCCGTTTTAAACCCCAAGGCCGTTCATCCGTgg                                             >  1:572381/1‑151 (MQ=255)
                                                                                                aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGGCGGGTCGTCGAGAAAAAAAAAGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAACCAATGGGCGGTCATCGGTGGCgg                                          >  2:523805/1‑150 (MQ=255)
                                                                                                 tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGGCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCCGTTTAAAACCCAAAGGCAggttctcggggggggg                                         >  2:483212/1‑136 (MQ=255)
                                                                                                   cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAACCCGGGAATCGGGGTGGGGGAAAAAAGAAAGGCCGGTTCATGGATTAGCAAACACGCCCAGGGCCACCGCCGTTTTAAAACTAAAGGCAGTTTAACGgggggggggg                                       >  1:586747/1‑151 (MQ=255)
                                                                                                   cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTCAACGCCGGAAGGCGGTTTGTCGAGAAAAAGAATGTCCGGGTCATGCATCAACGAACAGGCCCGCGCCCGCCGCTGGTTTAAACCTAAAGGGAGGTTATTGGGGGCGgggg                                       >  1:1141383/1‑151 (MQ=255)
                                                                                                     gTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGTAACGCCGGGAGGCGGGTTGTCGAGAAAAAAAAAGTCCGGGTCATGCAACAGCGAAAAGGCCCGCGCCCGCCGCTGTTTTAAACCCAAAGGCCGGTCAAcgggggcggggggg                                     >  1:146492/1‑149 (MQ=255)
                                                                                                          tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGAt                                <  1:842453/151‑1 (MQ=255)
                                                                                                            tgCAGCCAAAATTCACGGCGGGTGAGGGTGTCAACGCCGGAAGTCGGTTCGTCGATAACCAGAATGTCCGGTTCATGCATCAGCGACCAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTTGCGTGGGAGGCGATAc                              <  2:801238/151‑1 (MQ=255)
                                                                                                            tgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATAc                              <  1:199515/151‑1 (MQ=255)
                                                                                                            tgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATAc                              <  2:318577/151‑1 (MQ=255)
                                                                                                                 ccAAAATTCACGGCGGGTGAGGGGGTCACCGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTc                         <  2:614911/151‑1 (MQ=255)
                                                                                                                 ccAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTc                         <  2:284696/151‑1 (MQ=255)
                                                                                                                                     gggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGCTACTTTTCAGGCCGAACGCCTCGCTCAt     >  2:1014452/1‑151 (MQ=255)
                                                                                                                                      ggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGACACTTTTCAGGCCGAACGCCTCGCTCATg    >  2:353977/1‑151 (MQ=255)
                                                                                                                                       gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATg    <  1:1100380/150‑1 (MQ=255)
                                                                                                                                       gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATg    <  1:715007/150‑1 (MQ=255)
                                                                                                                                        ggTCAACCCCTGATGTCGCTTCGACGAGAAACCGAATGTCCGGTTCATGCATCTGCGAACAGGCCAGCCCCAGCCGCTGTTTAAAACCTAATGGCAGTTTATTGGTGACATCGGCGGCGATACTTTCGAGGGCCAACGACTCGCTCATGCg  >  2:327776/1‑151 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
ATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCG  >  NC_000913/827369‑827655

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