Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,903,692	0	C	A	43.2%	62.9 / 17.4	37	intergenic (+318/+1055)	ygcE/queE	putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base C (0/21);  new base A (16/0);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 7.77e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAAC  >  NC_000913/2903544‑2903840                                                                                                                                                     |                                                                                                                                                       ttttACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAcc                                                                                                                                                       >  2:135324/1‑150 (MQ=255)             tatTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAACACCC‑AATATcc                                                                                                                                           >  1:389484/1‑149 (MQ=255)              atTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑c                                                                                                                                        <  1:264993/150‑1 (MQ=255)                     acacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACaa                                                                                                                                 <  1:366403/150‑1 (MQ=255)                      cacTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACaaa                                                                                                                                <  1:149860/150‑1 (MQ=255)                                    aGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGa                                                                                                                  <  2:110945/150‑1 (MQ=255)                                     gAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCC‑AATAAACT‑CCAAAAAAAAAAAAAAATCCAGAtt                                                                                                                >  1:432375/1‑151 (MQ=255)                                      aaTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACCC‑ATTAACCC‑CAAAAAAAAAAAAAAAATCCAAAttt                                                                                                               >  1:166241/1‑151 (MQ=255)                                       aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCC‑ATTTACCT‑CCAAACAAAAACAAAAATTCAGAttt                                                                                                               >  1:246699/1‑150 (MQ=255)                                       aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCC‑ATTTACCT‑CCAAACAAAAAAAAACTTCCACAttt                                                                                                               >  1:91124/1‑150 (MQ=255)                                       aTAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATtta                                                                                                              <  2:193868/151‑1 (MQ=255)                                           aaGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACCCCC‑ATTTGCCT‑CAAAAAAAAAAAAAAACTCCACAattatt                                                                                                            >  1:36062/1‑149 (MQ=255)                                                 aTACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAAcc                                                                                                    <  1:515956/151‑1 (MQ=255)                                                              tttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACCCAATTACCTC‑CAAAAAAAAAAAAAATTTCCA‑ATTTTTTATTTTAACAAAGaaaaaaaa                                                                                        >  1:699771/1‑150 (MQ=255)                                                               ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACaa                                                                                        <  2:255369/149‑1 (MQ=255)                                                               ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAg                                                                                       <  2:670188/150‑1 (MQ=255)                                                                 aTCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTTTTTAAAAAAACTC‑CTTTACCT‑CAAAAAAAAAAAAAAAATCCAAATTTTTTTATTTACCaaaaaaaaaaaaactt                                                                                    >  1:367265/1‑137 (MQ=255)                                                                           aCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAAACCAC‑ATTAACTT‑CAAAAAAAAAAAAAAACCCCAAAATTATTTATTTACAAATTAAAAAAAAACTTaaaaaaaa                                                                            >  2:5968/1‑147 (MQ=255)                                                                              tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACt                                                                        <  1:595785/150‑1 (MQ=255)                                                                                aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGtt                                                                     <  2:204078/151‑1 (MQ=255)                                                                                aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGt                                                                      <  2:575866/150‑1 (MQ=255)                                                                                       aaCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAATCTTAAAAAATATACTGTTATTCTaa                                                              <  2:79186/151‑1 (MQ=255)                                                                                        aCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACACCC‑CTTTACCT‑CAAAAAAAAAAAAAAATTCCCGATTTTTTTATTTAACAAAAAATAAAAAACTTTAAAAAAATACTGTttttttta                                                              >  1:497594/1‑146 (MQ=255)                                                                                         cTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATTTAAAAACAACC‑ATTATTCT‑CAAAAAAAAGACAAACCTCCCGAATTATTACTTTAACAATAAAAAAAAAAAATAAAAAAAATACTTTTTTTCTtaa                                                             >  2:251292/1‑150 (MQ=255)                                                                                          tGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTC‑ATTAGCCT‑CAAAAAAAAAAAAAAATTTCCGATTTATTATTTAAACAATGAATAACAAACTTAAAAAAAATAATGTTAtttt                                                                >  1:511347/1‑144 (MQ=255)                                                                                             atatCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaa                                                        <  2:488603/151‑1 (MQ=255)                                                                                                             gAACTATTCATCAACAGCGTATAGAGGCGGTTTTGTAAAAACACCC‑ATTATCCTCCAAAAAAAAAAAAAAATTCC‑GAATTTTTTCTTAAACAAAAAAAAAAAAAATTAAAAAAAAAAATTTTTTTTTTAAAAAAAAAATaaaaag                                              >  1:347325/1‑141 (MQ=255)                                                                                                              aaCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTC‑ATTAACCT‑CAAAAAAAAAACAAACTTCCAGATTTATTTATTAAACAAAGAATAAAAAACTTAAAAAAAATACTGTTTTTTTAAAAAAAAAAATAAAAGGCATTaa                                        >  1:145020/1‑150 (MQ=255)                                                                                                                   ttCATCAACAGCGTATAGAGGCGGTTATGTAAAAACACTC‑ATTAGCCT‑CCAAAAAAAAAAAAAAATCCCGATTTTTTTCTTTAACAAAGAAAAAAAAACTTTAAAAAAATACTGTTATTTTTAAAAAAAAAATAAAAAGCATTaa                                        >  1:687264/1‑145 (MQ=255)                                                                                                                     catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAatat                                <  1:251354/151‑1 (MQ=255)                                                                                                                     catcaACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAata                                 <  2:556557/150‑1 (MQ=255)                                                                                                                           cAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaa                          <  1:386233/151‑1 (MQ=255)                                                                                                                            aGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTaaaaaa                         <  1:348530/151‑1 (MQ=255)                                                                                                                              cGTATAGAGGCGGTTATGTAAAAACACTC‑CTTAGCCT‑CCAAAAAAAAAAAAACTTTCCGATTTATTACTTAACAAAAGAAAAAAAAAATTAAAAAAAATAATGTTATTCTTAAAAAAAAAATAAAAGGCATTAATTTAATATTaaaaaaa                        >  1:271542/1‑150 (MQ=255)                                                                                                                                      ggCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGAtt               <  2:305498/151‑1 (MQ=255)                                                                                                                                       gCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTa              <  2:283481/151‑1 (MQ=255)                                                                                                                                        cGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTaa             <  2:291093/151‑1 (MQ=255)                                                                                                                                                   aaCCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAAc  <  1:515205/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       TTTTACAAGGAATATTAACTACACTCGATGTCTTTAGAATAAAAGACAATACGAATCGCACTTTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTC‑ATTAGCCT‑CAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTAATTCAATATTAAAAAAATAATAGATTAAAATTTCTTAAC  >  NC_000913/2903544‑2903840

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