Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913898,1920CT30.0% 42.4 / 10.1 31V68V (GTC→GTT) ybjOputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (15/6);  new base T (0/9);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 7.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGATGCTGCAGAAGCTGCCGGATATGCTCATCCTGATGC  >  NC_000913/897946‑898430
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
cAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCtttt                                                                                                                                                                                                                                              <  1:175510/251‑1 (MQ=255)
  aCCTGGCTACATAATGCCGCGCATGACGCGGCATTGTTTTAATGAAAGCCGAAACGTTGGAATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGAGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGCTTTTCAATACGCTGGGCGTGCGCGGGATTGGCGAGTTCCATCATCGCAGCCTACAAACCTGGACATTAACGCTGGCCttttt                                                                                                                                                                                                                                             >  1:202232/1‑250 (MQ=255)
             tAATGCCGCGCATGTCGCGGCATTGTTTTCATGGCAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGGTCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGGGTTTATTCATCGCAGAGTACAAACCTGGAGTTTAACGCCTGTCTTTTTAAGCAggtgtgg                                                                                                                                                                                                                                 >  2:305091/1‑245 (MQ=255)
                                                                   tttAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGGGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACACACCTGGCGGTTAACGCTGGGCTTTTTAAGCAGTCTGGGGCTTGTTTTCATTGCGATCTG‑GTG‑TGCGTTTTCACTGGTGGTAGGGCGTCGCTg                                                                                                                                                                         >  2:1191/1‑251 (MQ=255)
                                                                      aaGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTggg                                                                                                                                                                       <  1:175499/250‑1 (MQ=255)
                                                                                      ctcATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACCAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGACTGGTGCTGGTTTTCATTGTGATCTG‑GTG‑TGCGTTTTCACTGGTGACAGGGCGTCGCTGGGCGCGCTGGATATATctg                                                                                                                                                      >  2:312385/1‑251 (MQ=255)
                                                                                        cATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGCGTTCATTCATCGCAGCGTACACACCTGGGGTTTAACGCTGGTCTTTTTAAGCCGTCTGGTGCTGGGTTTCATTGAGATCTG‑GTG‑TGCGTTTTCCCTGGGTAAAGGGCGTCGCTGGGCGCGCTGGCTACATccgcc                                                                                                                                                    >  2:328899/1‑247 (MQ=255)
                                                                                              cGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGacac                                                                                                                                               <  1:149563/250‑1 (MQ=255)
                                                                                                    aaTGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGGCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCAc                                                                                                                                        >  1:185604/1‑251 (MQ=255)
                                                                                                       gtgCCTGCGCTGGGGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCAccg                                                                                                                                      <  2:254953/250‑1 (MQ=255)
                                                                                                                cTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGtt                                                                                                                             >  1:70963/1‑250 (MQ=255)
                                                                                                                     gtaGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGATTTCATTCATCGCAGCGTACAACCCTGGAGTTTAACGCTGGTTTTTTTAAGCAGTCGGGTGCTGGTTTTCATTGAGATCTG‑GTT‑TGCTTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGTTACTtgt                                                                                                                       <  1:1191/249‑1 (MQ=255)
                                                                                                                              ggcgCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGGGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGGGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACCAATCACCGACGCCAGTTACGTGTGGGCGGCtt                                                                                                              >  2:38397/1‑251 (MQ=255)
                                                                                                                               gcgcTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTc                                                                                                             >  1:119569/1‑251 (MQ=255)
                                                                                                                                                                gTGCTGATGCTTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCAGAGCGTACAAACCTGGAGTTTAACGCTTGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTTAAAGGGCGTCACTGGGCGCGTTGGCCTTATCTGCTGACCCAAATCACCGCCGCAGGTTACTTGTGGGCGGCTTCGGTGGGGGTCGGTTATCCGGAGCTGTTCAGcct                                                                            >  2:267858/1‑249 (MQ=255)
                                                                                                                                                                        gatTTTAAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCCTTCATCGCAGCGTACAAACCTGGGGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCCTTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGATGACACCAATCACCTCCGCAAGTTACTTGTGGGCGGCTTCGCTGTGGTACGGTTATCCGGAGCTGTTCCGCCTTCCCGGTg                                                                    >  2:282168/1‑251 (MQ=255)
                                                                                                                                                                              cAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCa                                                               >  1:213084/1‑250 (MQ=255)
                                                                                                                                                                                aTACGCTGGGGGTGCGCTGTATTGGCGAGTTCATTCATCGAGGCGTACAAGCCTGGAGTTTAACGCTGGTTTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAGGGGCGTCGCTGGGCGCGGTGGCTATACCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCaaaa                                                            <  2:213084/251‑1 (MQ=255)
                                                                                                                                                                                            ggcgtggTATGGGCGAGTTCTTTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTTTTTTTAAGAAGTCGGGTGGGGGTTGTCATTGAGATCTG‑GTG‑TGCGTTTTCATTGGGGAAAGGGCGTCGCGGGGCGCGTGGGCTATATCTGATGACTAAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTc                                                <  1:312385/246‑1 (MQ=255)
                                                                                                                                                                                               gcgttATTGGCGAGTTCATTCACCGCAGCGTACAAACCGGGAGTTTAACTCTGGTTTTTTTAAGCATTCTGGTCCTGTTTTTCATTGAGATCTG‑GTG‑TTCTTTTTCCCTGTTAAAAGGGCGTCGCGGGGCGCGCTGGCTATATCTGCTTCCACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCAt                                             <  1:100736/247‑1 (MQ=255)
                                                                                                                                                                                                     ttGGCGAGTTTTTTCGTCGCAGCTTAAAAACCGGGTGTTTAACGCTGTTTTTTTGAAGAATTTTGGTGTGGGTTTTC‑TTGGAAACTGTGTG‑GGCTTTTTCACTGGGGAAGGGGCGTTGCGGGGCCCGTGGGTTATATCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTGATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCt                                        <  1:321318/250‑1 (MQ=255)
                                                                                                                                                                                                       ggCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGAt                                     >  1:281548/1‑251 (MQ=255)
                                                                                                                                                                                                        gCGAGTTCTTTCAGCGCAGCGTATAAACCGGGAGTTTAACTCTGTTTTTTTTAAGCAGTCGTGTGCTGGTTTTCATTGAATTCTG‑GTG‑TGCTTTTTCACTGGTGAAAGGGCGGCGCTGGGCGCGTTGGCTATATTTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGAtg                                    <  1:282777/251‑1 (MQ=255)
                                                                                                                                                                                                           aGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGAtgctg                                 >  2:215807/1‑251 (MQ=255)
                                                                                                                                                                                                                       gcatcGTAAAAACCTGGAGTTTACCGCTTTTGTTTTTAAGTAGTTTGGGGCTGTTTTTCATTGAGATCTG‑GTT‑TGCTTTTTCATTGTTGAAAGGGCGTCGCTGGGCGCGTTGGCTATATCTGCTGACCCAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGTGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGATGCTGCAGAAGCTGCCg                     <  1:256534/247‑1 (MQ=255)
                                                                                                                                                                                                                       gcagcttccAAACCTGGATTTTAACTCTGTTTTTTTTAAGCAGTCTGGTGCTGTTTTTCATTGAGTTCTG‑TTT‑TGTGTTTTTACGGTTGAAAGGGCGTCGCTGGTCGCGTTGGCTATTTCTGCGGACCCAAACCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGTTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGATGCTGCAGACGCTTCCg                     <  1:256539/243‑1 (MQ=255)
                                                                                                                                                                                                                            gtaaaaaCCTGGAGTTAAAGGCGGGTCTTTTTAAGCAGTCGTGTGCTGGTTTTTATTGAGATCTG‑GTG‑TGCGTTTTCACTGGGAAAAGGGCGTCGCTGGGCGCGCTGGTTATATCTGCTGACACAAATCGCCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGATGCTGCAGAAGCTGCCGGatat                <  1:43452/247‑1 (MQ=255)
                                                                                                                                                                                                                            gTACAACCCTGGAGTTTAACGCTGGTTTTTTTAAGCAGTCGGGTGCTGGTTTTCATTGAGATCGG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGATGCTGCAGAAGCTGCCGGatat                <  1:215807/251‑1 (MQ=255)
                                                                                                                                                                                                                            gTACAAACCTGGAGTTTAACGTGGGGTTTTTTAAGAAGTCTGGTGCTGGTTTTCATGGAGATCTG‑GGG‑TGCGTTTTCTTTGGTGAAAGGGCGGCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGATGCTGCAGAAGCTGCCGGatat                <  1:38397/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGATGCTGCAGAAGCTGCCGGATATGCTCATcc        >  1:267736/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ttaacgtgggTCTTTTTAAGCAGTTGGGTGCTGGTTTTCTTGGAGATC‑G‑GTGGTGCGTTTCCACTGGTAAAAGGGCGCCGCGGGGCGCGCGGGCAATATCTGCTGACAAAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGATGCTGCAGAAGCTGCCGGATATGCTCATCCTGAtgc  <  2:70963/242‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
CAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTTTCAATACGCTGGGCGTGCGCGGTATTGGCGAGTTCATTCATCGCAGCGTACAAACCTGGAGTTTAACGCTGGTCTTTTTAAGCAGTCTGGTGCTGGTTTTCATTGAGATCTG‑GTG‑TGCGTTTTCACTGGTGAAAGGGCGTCGCTGGGCGCGCTGGCTATATCTGCTGACACAAATCACCGCCGCAAGTTACTTGTGGGCGGCTTCGCTGGGGTACGGTTATCCGGAGCTGTTCAGCATTCCCGGTGAATCAAAACGTGAAATCTTCCATAGCCTGATGCTGCAGAAGCTGCCGGATATGCTCATCCTGATGC  >  NC_000913/897946‑898430

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