Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,536,197	0	A	T	45.3%	59.3 / 43.4	53	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 (24/0);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 4.17e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGT  >  NC_000913/4536047‑4536341                                                                                                                                                        |                                                                                                                                                   cATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTa                                                                                                                                                    <  1:760912/151‑1 (MQ=255)       aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AAt                                                                                                                                              <  2:638814/150‑1 (MQ=255)       aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATg                                                                                                                                             <  2:473054/151‑1 (MQ=255)        cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATga                                                                                                                                            <  1:595400/151‑1 (MQ=255)              tagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTGG‑AGTgaaaaaa                                                                                                                                       >  2:449304/1‑150 (MQ=255)               agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATgaaagaaa                                                                                                                                      <  2:353200/150‑1 (MQ=255)               agtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGaaaa                                                                                                                                     <  1:508948/151‑1 (MQ=255)                 taAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTCACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGCTTTTTTG‑AAATAAAGaaaa                                                                                                                                     >  1:629548/1‑149 (MQ=255)                   aaGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATgca                                                                                                                                 <  1:713619/151‑1 (MQ=255)                    aGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATgcaa                                                                                                                                <  1:838562/151‑1 (MQ=255)                      ggCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTGG‑GAATAAAGAAAAtaaaa                                                                                                                                >  1:758379/1‑145 (MQ=255)                       gCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATgcaagca                                                                                                                             <  2:161052/151‑1 (MQ=255)                          tatTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGCTTTTTTG‑AATGAAAAAAAAAacaaacaa                                                                                                                            >  2:800736/1‑149 (MQ=255)                           atTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTGG‑AATGAAAAAAAAAgcaagcaaaa                                                                                                                          >  1:60360/1‑148 (MQ=255)                                aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAataatataa                                                                                                                    <  2:257839/151‑1 (MQ=255)                                  cAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTTG‑AATGAAAAAAAAAGCAAACaataaaa                                                                                                                       >  1:239327/1‑144 (MQ=255)                                  cAATGTTAAGGAGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAAtaatataata                                                                                                                  <  2:261892/151‑1 (MQ=255)                                   aaTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATAtaataa                                                                                                                 <  2:814743/151‑1 (MQ=255)                                   aaTGTTAAGGCGACTTTCCGTTTTATTACGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATAtaataa                                                                                                                 <  2:62300/151‑1 (MQ=255)                                        tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGG‑AAT‑AAAAAAAAAACAAACAAAAAaataataattgg                                                                                                             >  2:613942/1‑147 (MQ=255)                                                  ttCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCt                                                                                                    <  1:692418/149‑1 (MQ=255)                                                   tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCttt                                                                                                  <  2:349257/150‑1 (MQ=255)                                                      gTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGTAGCTTTTTTG‑TATGAAATAAAAAGCAAGGAAATATATaataatagtgttattttttt                                                                                                  >  2:52987/1‑133 (MQ=255)                                                              gCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTa                                                                                      <  2:168837/150‑1 (MQ=255)                                                                    cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCata                                                                                <  2:385863/150‑1 (MQ=255)                                                                     gCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATatc                                                                              <  2:132147/151‑1 (MQ=255)                                                                         ttCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTGG‑AATGAAAAAAAAAGTAAGAAATAAAATAATAAAAGGCTCAATTTTTTGAAAA‑ACAATAAACa                                                                                  >  1:773277/1‑143 (MQ=255)                                                                            cTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTG‑GATGAAAAAAAAAGCAAACAATAATATAATAAAAGGCTAATTTTTTTTTACA‑TCATTTAACATATAATCATCg                                                                       >  1:574162/1‑151 (MQ=255)                                                                             tGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGTTTTTTTG‑AATGAAAAAAAAAGCAAGCAAAAAAATAATAAAGGCCTAATTTTTTTTAACA‑ACATTAAACAAATAAACATCGt                                                                      >  2:499338/1‑151 (MQ=255)                                                                                      aTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTTG‑GATGAAAGAAAAAGCAAGCAATAATATAATAATAGCCTAATTTTTTTTAACA‑ACATTAAGCATATCATCATCGTCTATTTTcc                                                             >  1:725649/1‑151 (MQ=255)                                                                                       tCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGCTTTTTTG‑AATTAAAAAAAATGCAAGCAAAAATATAATAAAGGGCTAATTTTTTTTAACA‑TCATTAAGCATATAATCATCGTCTATTTTcc                                                             >  2:388768/1‑150 (MQ=255)                                                                                        cTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTTG‑AATTAAAAAAAAAGCAAACAATAAAAAAATAAAAGCCTAATTTTTTTTTAAA‑TCATTAAACATAATAACAACGTTTAtttt                                                               >  2:64424/1‑146 (MQ=255)                                                                                          tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTGG‑GATTAAAGAAAATGCAAACAATAATATAATAATGGCCTAATTTTTTTTAACA‑TCATTAAACATATCATCATCGTCTTTTTTCCAGa                                                          >  1:176538/1‑150 (MQ=255)                                                                                          tATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGcc                                                             <  2:791528/147‑1 (MQ=255)                                                                                                 tACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGGTTTTTGG‑GATGAAAGAAAAAGCAAGCAATAATATAATAAAAGCCTAATTTTTTTTAAAA‑TCATTAAACATATCATTATCGTCTATTTTCCAGaa                                                         >  1:422169/1‑144 (MQ=255)                                                                                                  aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGa                                                 <  1:682157/151‑1 (MQ=255)                                                                                                  aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTG‑GATGAAAGAAAAAGGAAACAAAAATAAAATAAATGCCTAATTTTTTTTAAAA‑ACATTAAACATATCATTATTCTCTATTTTCCAGAATccc                                                     >  1:461219/1‑146 (MQ=255)                                                                                                   ccTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGAc                                                <  2:250675/151‑1 (MQ=255)                                                                                                     tAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGCTTTTTTG‑AATGGAAAAAAAAACAAGAAATAAAATAAAAATTGGCTCAATTTTTTTAACA‑TCATTAAAAATAAAATCATCATCTATTTTTCAGAATTCGTCGACtt                                              >  2:835218/1‑151 (MQ=255)                                                                                                      aaTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTta                                             <  2:678946/151‑1 (MQ=255)                                                                                                         cATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTtata                                           <  1:79780/150‑1 (MQ=255)                                                                                                          aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAAt                                         <  1:366576/151‑1 (MQ=255)                                                                                                          aTTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGCTTTTTGGAAAGGAAAAAAAAAACAA‑CAATAATATAATAAAAAGCTATATTTTTTGAACATTCCTTAA‑CATAAAATCATCCTCGTTTTTCCACGATCCCTCCGCTtat                                            >  1:769117/1‑148 (MQ=255)                                                                                                            tttCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGGGGTTTTTTG‑GATTGAAGAAAAAACAAGCCAAAAAATAATAAAAGGCTAATTTTTTTTTAAA‑AAAATAAGAAAATCATTATCGTTTATTTTTCAGAATTCGTCGACTTATaa                                          >  2:290490/1‑148 (MQ=255)                                                                                                                    gCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAAc                                <  1:344889/150‑1 (MQ=255)                                                                                                                      aTAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTTGG‑AATGAAAGAAAAAGCAAGCAAAAATATAATAAATGCCTAATTTTTTTTAACA‑TCATTAAACATATTATCATCGTCTATTTTCCAGAATCCGTCGACTTATAAACGCCAGAACgg                              >  1:692844/1‑150 (MQ=255)                                                                                                                        aaTGCGCCTGGGATTAATGGC‑GGGGGCTTTTTGG‑AATGTAAGAAAATTCAAGAAAAAAAATAATAATTGCCCAATTTCTTTGTAAA‑CTATTAAACATATCTTCATCGATTATTTTCCCTAATCCCTCGTGTTATAATCTCCCGa                                  >  1:410502/1‑144 (MQ=255)                                                                                                                               cTGGGATTAATGGC‑GGGGGGTTTTTGG‑GATGAAAAAAAAAGCAAGCAATAATATAATAATAGCCTAATTTTTTTTTACA‑TCATTAAACATAATATCATCGTCTATTTTCCAGAATCCGTCGACTTTTAATTACCAAAAAGGCCTATCAgg                     >  2:70197/1‑150 (MQ=255)                                                                                                                                   gATTAATGGC‑GGGTGCTTTTTGG‑GATGAAAAAAAAAGGAAACAAAAAAATAATAATGGGCTAATTTCTTTTTACA‑ACATTAAGCATATCATAATAGTCTTTTTTCCAGAATCCGTCGGCTTATAATCTCCAGAAAGGCATATCAGGACgg                 >  2:232906/1‑150 (MQ=255)                                                                                                                                          ggC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg          >  1:819971/1‑150 (MQ=255)                                                                                                                                          ggC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCg          <  2:538251/150‑1 (MQ=255)                                                                                                                                              gggTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAAc       <  2:158634/150‑1 (MQ=255)                                                                                                                                                  gCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtgt  >  1:248207/1‑151 (MQ=255)                                                                                                                                                  gCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtg   <  1:155795/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   CATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGC‑GGGTGCTTTTAGG‑AATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACA‑TCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGT  >  NC_000913/4536047‑4536341

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