Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,536,197	0	A	T	43.1%	68.9 / 33.6	51	intergenic (+166/+417)	yjhR/nanS	pseudogene, helicase family; putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base A (2/27);  new base T (22/0);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.77e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCA  >  NC_000913/4536050‑4536347                                                                                                                                                    |                                                                                                                                                        gATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGa                                                                                                                                                      <  1:794959/151‑1 (MQ=255)    aCAGGAGTAGTAAAGGGCATGTTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATg                                                                                                                                                   >  1:352055/1‑151 (MQ=255)      aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGAATg                                                                                                                                                   >  2:9617/1‑149 (MQ=255)       ggAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaa                                                                                                                                                <  2:45943/151‑1 (MQ=255)        gAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATgaaa                                                                                                                                                >  2:263176/1‑150 (MQ=255)          gtagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaa                                                                                                                                             <  2:320739/151‑1 (MQ=255)          gtagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATgaaagaa                                                                                                                                             <  1:559820/151‑1 (MQ=255)               aaaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATg                                                                                                                                         <  1:627522/150‑1 (MQ=255)                   ggCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagc                                                                                                                                    <  1:929212/151‑1 (MQ=255)                     cATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAGAAAAAgcaaaca                                                                                                                                   >  1:397924/1‑150 (MQ=255)                      atatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAGAAAAAacaagcaa                                                                                                                                  >  1:930210/1‑150 (MQ=255)                        atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATTaaaaaaaaaaaaaacaata                                                                                                                                >  2:213251/1‑139 (MQ=255)                        atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaata                                                                                                                                <  1:342408/150‑1 (MQ=255)                        atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaata                                                                                                                                <  2:1008838/150‑1 (MQ=255)                          tttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaataat                                                                                                                              <  1:112050/150‑1 (MQ=255)                           ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATGAAAGAAAAAGCAAGCAaaaata                                                                                                                             >  1:781551/1‑150 (MQ=255)                            tACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatat                                                                                                                            <  2:374607/150‑1 (MQ=255)                             aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatata                                                                                                                           <  1:646150/150‑1 (MQ=255)                              ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAGCAAGCAaaaatataa                                                                                                                          >  1:154531/1‑150 (MQ=255)                              ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGAATGAAAGAAAAAGAAAAAAaaaatataa                                                                                                                          >  1:915115/1‑150 (MQ=255)                              ccAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGAAAAAAAAAACAAGAAaaaaaataa                                                                                                                          >  2:1007666/1‑150 (MQ=255)                               cAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAAtaatataata                                                                                                                        <  1:305367/151‑1 (MQ=255)                                aaTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAtaataa                                                                                                                       <  1:938965/151‑1 (MQ=255)                                aaTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGAGCTTTTAGGAATGAAAGAAGATGCAAGCAATAATAtaataa                                                                                                                       <  2:315148/151‑1 (MQ=255)                                  tGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGGTTTTGGGATGAAAAAAAAAACAAGaaaaaaaaaaaaaa                                                                                                                       >  2:768085/1‑141 (MQ=255)                                    ttAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATgg                                                                                                                    <  2:990797/150‑1 (MQ=255)                                          cGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATGAAAAAAAAAGAAAGCAAAAATATAATAATGGGCTAAtt                                                                                                             >  2:584959/1‑151 (MQ=255)                                          cGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAAtt                                                                                                             <  1:156600/151‑1 (MQ=255)                                           gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAAttt                                                                                                            <  1:252410/151‑1 (MQ=255)                                                          tGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAAAGAAAAAAAAAGCAAGCAAAAAAATAAAAATAGCCTAATTTTTTTTTAAATCAtt                                                                                              >  1:694154/1‑150 (MQ=255)                                                            cGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAg                                                                                           <  1:756864/151‑1 (MQ=255)                                                             gAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGAAGAAAAAAAAAGCAAGAAATAATATAATAAAAGCCCAATTTTTTTTAAAATCAtt                                                                                              >  1:741142/1‑147 (MQ=255)                                                               gTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCat                                                                                        <  2:833893/151‑1 (MQ=255)                                                                       tCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTGGAATGAAAAAAAAAGCAAGCAAAAAAATAAAAAAGGGCTAATTTTTTTTTAAATCATTAAAAatataaaca                                                                                 >  2:369952/1‑145 (MQ=255)                                                                                 aaaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTg                                                                     <  2:741418/151‑1 (MQ=255)                                                                                   aTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAGAAAAAGCAAGCAATAAAATAATAATAGCCTAATTTTTTTTAACATCATTAAACATATCATCATCCTC‑TAtttttc                                                                   >  1:848296/1‑148 (MQ=255)                                                                                           gTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTTATGAAAAAAAAAGCAAGCAAAAAAATAATAATGGGCTAATTTTTTTGAAAATCATTAAACAAATTATAATTCTC‑TAtttttcaaaa                                                               >  2:116091/1‑140 (MQ=255)                                                                                             ttACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATGAAAAAAAAAGGAAGCAATAAAATAATAAAAGCCTAATTTTTTTTAAAATCATTAAACAAATTATCATCTTC‑TATTTTTCAGAATCCGt                                                          >  2:188387/1‑150 (MQ=255)                                                                                              tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTc                                                         <  1:497058/150‑1 (MQ=255)                                                                                              tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTc                                                         <  2:672631/150‑1 (MQ=255)                                                                                                   aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTCGACTta                                                   <  2:122043/151‑1 (MQ=255)                                                                                                          ttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCCTTTTGGGATGAAAAAAAAAAAAAACAATAATATAATAATGTGCTAATTTTTTTGTACATCAATAAGAATATTATAATCCTT‑TATTTTTCAGAAAACGTCGTCTTATaa                                                >  2:36004/1‑147 (MQ=255)                                                                                                                ttcATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAAAAAAAGCAAGCAATAATATAATAAAAGCCCAATTTTTTTTTAAATCATTAAGCATATTATTATCCTCTTTTTTTCCAAATTCCG‑CGACTTTTAAttaccacaa                                       >  1:774801/3‑142 (MQ=255)                                                                                                                tGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATGAAAAAAAAAGCAAGCAAAAATATAATAAAGGGCTAATTTTTTTTTAAATAATTAAACAAATAATCATCATC‑TATTTTCCAGAAACCGTCGACTTATAATCACCAGaa                                       >  2:86766/1‑150 (MQ=255)                                                                                                                 gCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACg                                     <  2:781477/151‑1 (MQ=255)                                                                                                                     aaTGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCat                                 <  1:263192/151‑1 (MQ=255)                                                                                                                        gcgcCTGGGATTAATGGCGGGGGCTTTTTGGGATGAAAGAAAAAACAAACAAAAAAATAAAAAATGCCTAATTTTTTTTAACATCATTTAACATATTATCATCGTC‑TTTTTTCCAGAATCCGTCGACTTATAATTGCCAGAAAgg                                    >  1:722384/1‑145 (MQ=255)                                                                                                                         cgcCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCa                             <  1:790969/151‑1 (MQ=255)                                                                                                                                       ggCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg                <  1:213267/150‑1 (MQ=255)                                                                                                                                              gCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgt        >  1:953018/1‑151 (MQ=255)                                                                                                                                                  ttAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGcc    >  2:327966/1‑151 (MQ=255)                                                                                                                                                    aGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCa  <  1:751856/151‑1 (MQ=255)                                                                                                                                                    aGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATAATAATAGTC‑TATTTGCCAGAATCAGTCGACTTATAATAGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGccc   >  2:93155/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                        GATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTC‑TATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCA  >  NC_000913/4536050‑4536347

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