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

ATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCG  >  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)
                                                   tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCCGAAGTCGGTTCGGCGGGAAAAAGAATGGCCGGTTCATGCATCAGCGAAAA‑gggc                                                                                        >  1:242523/1‑148 (MQ=255)
                                                    ccATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAACAGAAAGTCCGGGTTATGCATCAGCGGAAA‑Ggcccg                                                                                      >  1:310708/1‑149 (MQ=255)
                                                       tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑Ggccagcgc                                                                                   <  2:824295/151‑1 (MQ=255)
                                                       tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑Ggccagcgc                                                                                   <  1:343253/151‑1 (MQ=255)
                                                       tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATGCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑Ggccagcgc                                                                                   <  2:564184/151‑1 (MQ=255)
                                                       tAAAGAGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTAAACGCCCGGAGTCCGTTTGTCGAGGAAAAGAATGTTCGGGTCTTGGATTAGCGAAAA‑Agcccgcg                                                                                    >  1:575037/1‑150 (MQ=255)
                                                                   tgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGGGCCAGCCGCTGttt                                                                       <  1:932720/151‑1 (MQ=255)
                                                                    gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTCCGGGTCCTGCATCAGCGAACA‑GGCCAGCGCCCGCCCCCGTTTa                                                                      >  1:244453/1‑151 (MQ=255)
                                                                    gACCATCACCGTGACGCCTTCCTCTACCATGCTGTTGATTTGCAGCCAAAATTCACGGCGGGCGGGGGGGTTCATGCCCGGACTCGCGTCGTTTGGGAACACCAAGTGCGGTTCAATCATCAGCGCACA‑AGCCAGCGCCAGcccaggtttt                                                                      >  2:687777/1‑143 (MQ=255)
                                                                        aacACCGTGACGCCTGTCTCTACCATGCGGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATACATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAAc                                                                  <  2:61697/149‑1 (MQ=255)
                                                                           aCCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTGCGGGTCATGCCTCAGCGAAAA‑GGCCAGCGCCCGCCCCCGTTTTAAACCTa                                                               >  2:926446/1‑151 (MQ=255)
                                                                            ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGTGGGGGGCCACGCCCGGAGGCGGTTCGTCGGGAAACAGAAAGTCCGGTTCATGGATCAGCGAACA‑AGCCCGCGCCAACCGCCGTCTTTAACCtta                                                              >  2:253589/1‑149 (MQ=255)
                                                                                 aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                          <  1:764275/150‑1 (MQ=255)
                                                                                 aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                          <  1:838906/150‑1 (MQ=255)
                                                                                    ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGTGGAGTCGTCCAGAAACAGAAAGTTCGGGTTATGCATCAACGAACA‑GGCCAACGGCCGACGCTGTTTAAAACCTAATGTCCGtt                                                      >  2:818616/1‑151 (MQ=255)
                                                                                    ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGTCGGGTCGTCGGGAAAAAGAATGTCCGGGTCATGCCTCAGCGAACA‑GGCCCGCGCCCGCCGCCGTTTTAAAACTAATGGGAggt                                                      >  1:371503/1‑149 (MQ=255)
                                                                                    ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGGAGTCGGTTCGTCGAGAAACAGAAAGTCCGGGTCATTCATCAACGAAAA‑AGCCAGCGCCCGGCCCTGTTTTAAAACCAATTGCAggt                                                      >  1:121796/1‑149 (MQ=255)
                                                                                      tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCTGGAGTCGGTTCGTCGAGGAAAAGAATGTCCCGGTCCTGGCTCCGCGAAAA‑GGCCCGCGCCCGCCGCCGTTTAAAACCCAAACGCAggtta                                                    >  2:75150/1‑147 (MQ=255)
                                                                                      tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTTGGGTTTTCGAGAAAAAGAATGTCCGGGTCATGTCTCCGCGAACA‑AGCCAGCGCCAGCCGCTGTTTTAAAACTAATGGccggt                                                      >  1:540062/1‑145 (MQ=255)
                                                                                      tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGGGGGGGGCGACGCCCGAAGCCCGTTCTTCCGGAAAAAGAGTGGCGGGTGCCTCCATCTGCGAACA‑GGCCCGCGCCCCCCGCCGTTCAAAACCTAATGCCAGtgtt                                                    >  2:82671/1‑148 (MQ=255)
                                                                                        tctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGTAACGCCCGAAGGGGGTTTCGCTAGAAAAAGAGTGTTCGGGTCATTTGACAGCGGAAA‑GGACCGCGCGAGCCCGCGTTTTTAACACAATGGCAgcttatg                                                  >  2:766095/1‑145 (MQ=255)
                                                                                        tctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGGCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAAAA‑AGCCAGCGCCAGCCGCCGTTTAAAAACTAATGGCAggttat                                                   >  1:198089/1‑145 (MQ=255)
                                                                                           ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGTCAACGCCCGAAGTCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCAACAGCGAAAA‑GGCCCGCGCCCGCCCCTGGTTAAAAACTAAAGGGAGGTCATCggg                                               >  1:1066461/1‑150 (MQ=255)
                                                                                           ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAAGGGGGCAACTCGGGCCCACGGTTCCGTGGGAAACAGAAAGTTCGGTTAAAGCACCAGCCAACA‑CGCCAGGCCCCACCGCTTTTTTAAACCCAATCGGAGTTTAtaagt                                               >  1:822920/1‑147 (MQ=255)
                                                                                             aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTTGTCGGGAAAAAGAATGTCCGGTTCATGCCTCAGCGAACA‑GGCCCGCGCCCCCCCCCGTTTTAAACCCCAAGGCCGTTCATCCGTgg                                             >  1:572381/1‑151 (MQ=255)
                                                                                                aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGGCGGGTCGTCGAGAAAAAAAAAGTCCGGGTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAAACCAATGGGCGGTCATCGGTGGCgg                                          >  2:523805/1‑150 (MQ=255)
                                                                                                 tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGGCCGGGTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCCGTTTAAAACCCAAAGGCAggttctcggggggggg                                         >  2:483212/1‑136 (MQ=255)
                                                                                                   cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAACCCGGGAATCGGGGTGGGGGAAAAAAGAAAGGCCGGTTCATGGATTAGCAAACACGCCCAGGGCCA‑CCGCCGTTTTAAAACTAAAGGCAGTTTAACGgggggggggg                                       >  1:586747/1‑151 (MQ=255)
                                                                                                   cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTCAACGCCGGAAGGCGGTTTGTCGAGAAAAAGAATGTCCGGGTCATGCATCAACGAACA‑GGCCCGCGCCCGCCGCTGGTTTAAACCTAAAGGGAGGTTATTGGGGGCGgggg                                       >  1:1141383/1‑151 (MQ=255)
                                                                                                     gTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGTAACGCCGGGAGGCGGGTTGTCGAGAAAAAAAAAGTCCGGGTCATGCAACAGCGAAAA‑GGCCCGCGCCCGCCGCTGTTTTAAACCCAAAGGCCGGTCAAcgggggcggggggg                                     >  1:146492/1‑149 (MQ=255)
                                                                                                          tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGAt                                <  1:842453/151‑1 (MQ=255)
                                                                                                            tgCAGCCAAAATTCACGGCGGGTGAGGGTGTCAACGCCGGAAGTCGGTTCGTCGATAACCAGAATGTCCGGTTCATGCATCAGCGACCA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTTGCGTGGGAGGCGATAc                              <  2:801238/151‑1 (MQ=255)
                                                                                                            tgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATAc                              <  1:199515/151‑1 (MQ=255)
                                                                                                            tgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATAc                              <  2:318577/151‑1 (MQ=255)
                                                                                                                 ccAAAATTCACGGCGGGTGAGGGGGTCACCGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTc                         <  2:614911/151‑1 (MQ=255)
                                                                                                                 ccAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTc                         <  2:284696/151‑1 (MQ=255)
                                                                                                                                     gggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGCTACTTTTCAGGCCGAACGCCTCGCTCAt     >  2:1014452/1‑151 (MQ=255)
                                                                                                                                      ggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGACACTTTTCAGGCCGAACGCCTCGCTCATg    >  2:353977/1‑151 (MQ=255)
                                                                                                                                       gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATg    <  1:1100380/150‑1 (MQ=255)
                                                                                                                                       gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATg    <  1:715007/150‑1 (MQ=255)
                                                                                                                                        ggTCAACCCCTGATGTCGCTTCGACGAGAAACCGAATGTCCGGTTCATGCATCTGCGAACA‑GGCCAGCCCCAGCCGCTGTTTAAAACCTAATGGCAGTTTATTGGTGACATCGGCGGCGATACTTTCGAGGGCCAACGACTCGCTCATGCg  >  2:327776/1‑151 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
ATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACA‑GGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGCG  >  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: ‑