Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,903,715	0	G	A	32.0%	36.6 / 11.2	25	intergenic (+341/+1032)	ygcE/queE	putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (1/16);  new base A (8/0);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 8.32e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGATTTAAGAATCATACAAAT  >  NC_000913/2903569‑2903860                                                                                                                                                    |                                                                                                                                                   cGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGAc                                                                                                                                                  <  2:172533/149‑1 (MQ=255)        tttAGAATAAAAGACAATACGAATCGCTCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGGTTTAT‑ATAAAAACTCAGTATCCTCaaaaaaaaaaaaaa                                                                                                                                               >  1:381870/1‑145 (MQ=255)             aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGAttt                                                                                                                                   <  2:14081/151‑1 (MQ=255)              aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAAACACTCATTAACCTCAAAAAAAAAAAAAAA‑ATCCAGAttt                                                                                                                                   >  2:11335/1‑150 (MQ=255)                aaaaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACCAC‑CAATAACCCCAAAACAAAAAAAAACATTCCA‑ATttatt                                                                                                                                >  1:18082/1‑151 (MQ=255)                  aaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑ATCCAGATttagta                                                                                                                               <  1:69796/150‑1 (MQ=255)                       aaTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATAT‑AAAAAAACTCAATATACTCACAAAAAAAACAAAA‑ATTCCGATttatt                                                                                                                                >  1:21346/1‑144 (MQ=255)                                 gCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTT‑AAAAAAACACCTTAACCTCAAAAAAAAAAAAAAA‑ATCCCGAATTTTTTTTTAACCAAataata                                                                                                                >  2:326643/1‑150 (MQ=255)                                   aCTTTATCAATTTTCGCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACa                                                                                                             <  2:130820/151‑1 (MQ=255)                                   aCTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGT‑TAAAACACACAATTTTCCCAAAACAAAAAAAAAA‑AACTACAATTTTTATTTAACCAAACAATaaa                                                                                                              >  1:225332/1‑149 (MQ=255)                                       tATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGT‑AAAAACACACAATTTTCCCAAAAAAAAAAAAAAC‑ATTCAGATTTATTAATTTTACAATGaaataa                                                                                                              >  1:135892/1‑142 (MQ=255)                                        aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCtt                                                                                                        <  1:342158/151‑1 (MQ=255)                                        aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCtt                                                                                                        <  2:261332/151‑1 (MQ=255)                                         tCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTT‑AAAAACACTCATTAACCTTAAAAAAAAAACAAAA‑ATTCACATTTTTTATTTAAACAAAGAAAAAAAAACtt                                                                                                        >  2:323988/1‑150 (MQ=255)                                              tttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaaaaa                                                                                                   <  2:27634/150‑1 (MQ=255)                                                 cACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAatat                                                                                               <  1:225257/151‑1 (MQ=255)                                                         gTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTAt                                                                                       <  2:281104/151‑1 (MQ=255)                                                                                  aaGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCAt                                                               <  2:191199/150‑1 (MQ=255)                                                                                         aTTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCCCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCa                                                        <  2:93696/150‑1 (MQ=255)                                                                                          ttCATCAACAGCGTATAGAGGCGGTTTTTT‑AAAAACACCCCTTTACCCCCAAACAAAAACAAAC‑ATCCACATTTATTTTTTTAACAAAAAATAAAAAACTTAAAAAAAATAATTTTTTTTTTAAAAAAAAAATAAAAAGCATTAATTTaa                                                       >  2:233170/1‑150 (MQ=255)                                                                                                aaCAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaa                                                <  2:229754/151‑1 (MQ=255)                                                                                                   aGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                                             <  2:21349/151‑1 (MQ=255)                                                                                                     cGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAaat                                           <  2:78397/151‑1 (MQ=255)                                                                                                            agGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAt                                    <  1:42090/151‑1 (MQ=255)                                                                                                                    aTGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATGaaaagtg                            <  1:117584/151‑4 (MQ=255)                                                                                                                                               aCAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGATTTAAGAATCATACAAAt  >  2:354687/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGT‑AAAACCACTCATTAGCCTCAAAACAAAGACAAAC‑TTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAACGATTTAAGAATCATACAAAT  >  NC_000913/2903569‑2903860

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