Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,801,591	0	T	C	64.3%	26.7 / 23.5	28	Y327C (TAC→TGC)	thrS	threonine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (5/5);  new base C (7/11);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCG  >  NC_000913/1801348‑1801830                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                caCATATCATAGACTAAACGGATACACCCGTTAACTTCATCTCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCTGGCGCCTCAGGCCATGCAGCTAACCTGACGGCTCGTTACGGTGGCAGCTACCAATCTCGTCCATACGTAGCGGCAGATCGCGATTAGACTTTAGCTCCTGGTTGAAAATTTGTTCTTGACCCGGTCAGTTCATCGGCTTAATGCAGCAGTCACg                                                                                                                                                                                                                                          <  2:143603/250‑1 (MQ=255)                 aaCGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGAtgtg                                                                                                                                                                                                                          <  1:155126/251‑1 (MQ=255)                    ggATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATTCACGCGCGCGCATCAGGCCATGCAGCGAACCTGACTGCTCGTTTCGGTTGCAGCTACCAAACTCGGCCATACGTAGCGGCAGATCGCGATAAGACTTCAGCTCCTGGTGGAAAATTTGTACTTGACCCGGGCTGTTCATCGGCTGAATGCAGCATGCACGTTTCTCAGATGATgtggtg                                                                                                                                                                                                                       <  2:289340/251‑1 (MQ=255)                             cGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGc                                                                                                                                                                                                               >  2:152586/1‑250 (MQ=255)                             cGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCagaag                                                                                                                                                                                                                               <  2:65988/234‑1 (MQ=255)                             cGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCagaag                                                                                                                                                                                                                               >  1:65988/1‑234 (MQ=255)                                      cATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAg                                                                                                                                                                                                     <  1:152586/251‑1 (MQ=255)                                           cgAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCTAACCTGACGGCTCGTTACGGTGGCAGCTTCCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCATCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTGCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTc                                                                                                                                                                                                <  1:204606/251‑1 (MQ=255)                                                             cAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTc                                                                                                                                                                              >  1:78150/1‑251 (MQ=255)                                                               gaagaTATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTccc                                                                                                                                                                            >  2:110159/1‑251 (MQ=255)                                                                                        aaTCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACgg                                                                                                                                                   <  2:95304/251‑1 (MQ=255)                                                                                              cgcacgcGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAAAATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCttt                                                                                                                                             >  2:244012/1‑251 (MQ=255)                                                                                              agaatgcgaaTCAGGCCATGCAGCGAACCTGGCGGCTCGTTACGTGGGCAGCTTCTAATCTCGGCCATACGCATCGTGAGATCGCGATTAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTTTCAGAATATGGGGTGAACATTGCATCGTTGTAGTTGTCCCAGTGACCGGTTTTTTTCCTCAGGACACGGTCAATCATGAACGGTCCttt                                                                                                                                             <  2:137668/242‑1 (MQ=255)                                                                                                    cgcATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGTTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTc                                                                                                                                       <  1:100589/251‑1 (MQ=255)                                                                                                                      gAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCcatcat                                                                                                                                                         >  1:129277/1‑215 (MQ=255)                                                                                                                      gAACCTGACGGCGCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGAGAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTGCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGGCcatcat                                                                                                                                                         <  2:129277/215‑1 (MQ=255)                                                                                                                         ccTGACTGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAg                                                                                                                  <  2:303982/251‑1 (MQ=255)                                                                                                                               ggCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAga                                                                                                            <  2:254740/251‑1 (MQ=255)                                                                                                                                                gCTACCACACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTcc                                                                                            >  2:77307/1‑250 (MQ=255)                                                                                                                                                   aCCAAACGCGGCCATACGCAGCTGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCATAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACTGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAg                                                                                          <  2:174518/249‑1 (MQ=255)                                                                                                                                                                       gCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCcg                                                                    <  1:110159/251‑1 (MQ=255)                                                                                                                                                                                 gcgATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGcc                                                          >  2:205174/1‑251 (MQ=255)                                                                                                                                                                                     tAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCGGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGa                                                       >  1:300298/1‑250 (MQ=255)                                                                                                                                                                                                            aaTTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGCATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATCCCCGGCGCTTCTTCCt                               >  2:148798/1‑251 (MQ=255)                                                                                                                                                                                                                tGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACg                                                                                                                                                                  <  1:56117/116‑1 (MQ=255)                                                                                                                                                                                                                tGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACg                                                                                                                                                                  >  2:56117/1‑116 (MQ=255)                                                                                                                                                                                                                       gACCCGGGCAGTTCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCGTCCTGCATATGGTAc                    <  1:197150/251‑1 (MQ=255)                                                                                                                                                                                                                                          cTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCg  >  1:78165/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                TACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGCAGCGGCAGATCGCGATAAGACTTCAGCCCCTGGTTGAAAATTTGTACGTGACCCGGGCAGTTCATCGGCTTAATGCAGTATTCACGGTTCTCAGAAGATGTGGTGAACATTGCATCTTTGTAGTTGTCCCAGTGACCGGTTTTTTCCCACAGGACACGGTCCATCATGAACGGACCTTTAACTTCCTGATACTGGTACTCTTTCAGTTTAGAACGAACAAACACTTCCAGTTCACGGAAGATGGTCCAGCCGTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCG  >  NC_000913/1801348‑1801830

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