Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,487,163	0	C	T	48.9%	40.4 / 26.7	45	intergenic (+200/‑72)	hrpA/ydcF	putative ATP‑dependent helicase/DUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/22);  new base T (22/0);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 5.83e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTG  >  NC_000913/1487016‑1487306                                                                                                                                                    |                                                                                                                                                  tttttATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTtcct                                                                                                                                                 <  1:371368/150‑1 (MQ=255)       ttaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTATGAAAATAACATTTTTTCCttt                                                                                                                                            >  1:755076/1‑149 (MQ=255)         aaTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCaga                                                                                                                                        <  2:305522/151‑1 (MQ=255)         aaTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCag                                                                                                                                         <  2:7865/150‑1 (MQ=255)            tCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGt                                                                                                                                      <  1:739568/150‑1 (MQ=255)                 aaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAAGATTTTCTCCTTTCTGAGGGCTcc                                                                                                                                 >  1:194649/1‑150 (MQ=255)                      tAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAGGAAAAAACAATTTCTCCTTTTCGAGTGCGcc                                                                                                                                 >  1:197755/1‑145 (MQ=255)                        aGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCa                                                                                                                         <  2:527497/151‑1 (MQ=255)                            tAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAATAAAAATTTCTCCTTTTCGAGGGGTCCCCCtctctc                                                                                                                        >  1:418085/1‑146 (MQ=255)                               aaaTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACtattat                                                                                                                  <  1:566590/151‑1 (MQ=255)                                aaTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAGGAAAAAAGATTTTCTCCTTTTAGAGTGCTCCCCTTCTCCCtattat                                                                                                                  >  1:222905/1‑150 (MQ=255)                                  ttCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGAAAAACAATTTTCTCCTTTCTGAGGGGTCCGCCTCTCCCtattattt                                                                                                                >  1:889053/1‑148 (MQ=255)                                    cAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCa                                                                                                              <  2:172421/150‑1 (MQ=255)                                    cAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCAc                                                                                                             <  2:90067/151‑1 (MQ=255)                                       gAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCa                                                                                                          <  1:10040/151‑1 (MQ=255)                                              gAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAATAAAAATTTCTCCTTTCCGAGGGGTCCGCTTCTCCCTATTTTCTTACGCAGTA‑tttt                                                                                                  >  2:672453/1‑149 (MQ=255)                                               aaaaaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAAGAATTTCTCCCTTTTAAGTGCTCCGCCTCTCACCATTTtcttccccaggattttt                                                                                                  >  1:881655/1‑137 (MQ=255)                                                aaaaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGAAATAAGATTTTCTCCTTTCATAGGGCTCCTCTTCTCTCTATTATCTCACGCCAGA‑TTCtt                                                                                                 >  1:257831/1‑150 (MQ=255)                                                 aaaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTaa                                                                                               <  2:875856/151‑1 (MQ=255)                                                      cAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGaa                                                                                          <  1:418012/151‑1 (MQ=255)                                                       aaTTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAAc                                                                                         <  1:540851/151‑1 (MQ=255)                                                         ttCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGACATAAAAATTTTTCCTTTTCAAGAGGGCCCCCTCTCACCATTATTTCCCCCCGTA‑TTCtttaag                                                                                             >  2:239626/1‑141 (MQ=255)                                                           ctctGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAATAAGATTTTCTCCTTTTTGAGTGGTCCCCCTCTCCCTATTATTTCCCCCCGTA‑TTCtttaagg                                                                                            >  1:498857/1‑139 (MQ=255)                                                           ctctGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGACATACAATTTTCTCCCTTCTGAGTGGTCCGCTTCTCACTATTTTTTCACGCAGCA‑TTCTTAAggg                                                                                            >  2:630001/1‑144 (MQ=255)                                                            tctGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAACAATTTTTTTCTTTTTGAAGGCTCCCCTTTTTCCTCTTATCTTCCGCCATA‑TTTTTTAAGGAACGaaa                                                                                     >  2:858379/1‑148 (MQ=255)                                                               gctgcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAagg                                                                                  <  2:101205/150‑1 (MQ=255)                                                                  gcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAggagga                                                                              <  2:436699/151‑1 (MQ=255)                                                                  gcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAggagga                                                                              <  2:403563/151‑1 (MQ=255)                                                                    aaGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAAGATTTTCTCCTTTTAGAGTGCTCCGCCTCTCACTATTATCTCACCCAGGA‑TTTTTTAAGGAAAGATAAAGAGGaa                                                                             >  1:447641/1‑150 (MQ=255)                                                                    aaGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGACAAAAGAATTTCTCCTTTTTGAGTGGTCCCCCTCTCCCCATTTTTTCTCGCCGGATTTTTTAAGGGAAGAATAagagg                                                                                >  1:113452/1‑145 (MQ=255)                                                                                gTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGACATAAAAATTTCTTCTTTTTGAGGGGTCCGCCTTTTTCCATTATCTTCCACACTA‑TTCTTTAAGGGACGGTAAAGAGGAAACATGaaaa                                                                    >  2:748826/1‑145 (MQ=255)                                                                                   tATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAAAAATTTCTCCCTTCTGAGTGGTCCCCCTTTCACCATTTTTTCACCCCGTATTTTTTAAGGGAA‑GATAAAGAGGGAACATGAACATT‑Aaccc                                                            >  1:205406/1‑151 (MQ=255)                                                                                          aGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAATAAAAATTTCTCCTTTTCAAGTGCTCCCCCTTTCCCTATTTTTTTCCGCCGGATTTTTTAAGGGAAGAAAAAGGGGAAACAATAAAAATA‑ACCCCttt                                                         >  1:329469/1‑148 (MQ=255)                                                                                           gCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAACGATTTTCTCCTTTTTGAGGGGTCCGCTTCTCCCCTTTTTTTCCCCCAGGA‑TTTTTTAAGGAAAGGTAAAGAGGAAACACGAACATTAACCCCCTTTCCg                                                     >  1:381536/1‑151 (MQ=255)                                                                                                tCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCt                                                <  1:26275/150‑1 (MQ=255)                                                                                                    tAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCTTTc                                             <  2:742285/149‑1 (MQ=255)                                                                                                     aaaaGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCTTTCGcc                                          <  2:647569/151‑1 (MQ=255)                                                                                                      aaaGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAAAATTTTCTTCTTTTTTAGGGCTCCCCTTTTCCCCCTTTTTTTACGCAGGA‑TTTTTTATGGAAGGGTAAGGAGGAAACATGAAAATT‑TCCCCCttttccccccttt                                              >  2:158354/1‑136 (MQ=255)                                                                                                       aaGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCTTTCGCCgg                                        <  1:622245/151‑1 (MQ=255)                                                                                                       aaGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGAAATACAAATTTCTCCTTTTTTAGTGCTCCGCTTCTTACTTTTTTTTCACCCACTA‑TTCTTTAGAGAACGATAAAGAGGAAACATGAAAATTAACCCCCTTTcctcgcttttgc                                            >  2:611128/1‑138 (MQ=255)                                                                                                        aGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAAGATTTTCTCCCTTTAAAGGGGTCCGCCTCTCCCCATTATTTCACGCCGGATTTTTTAAGGGAA‑GATAAAGAGGAAACAAGAACAAT‑ACCCCCCTTTCCACGCCTTTGcc                                          >  1:294947/1‑148 (MQ=255)                                                                                                                   gCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATg                            <  2:737775/151‑1 (MQ=255)                                                                                                                       gCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCa                         <  1:674539/150‑1 (MQ=255)                                                                                                                       gCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCAt                        <  2:655536/151‑1 (MQ=255)                                                                                                                         ggAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATaa                      <  1:55065/151‑1 (MQ=255)                                                                                                                                             aCGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggctg  >  2:911532/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTCTTAAGGGAACGATAAGGAGGAACCATGAACATT‑ACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTG  >  NC_000913/1487016‑1487306

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