Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,745,089	0	A	C	30.3%	136.5 / 12.6	68	V84G (GTG→GGG)	trmD	tRNA m(1)G37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/12);  new base C (0/20);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 5.55e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGC  >  NC_000913/2744846‑2745320                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        gAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCCCCTTTgc                                                                                                                                                                                                                                  <  2:378552/251‑1 (MQ=255)           gTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCAccc                                                                                                                                                                                                                        <  1:722281/251‑1 (MQ=255)                 gTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCcgcg                                                                                                                                                                                                                  >  1:566297/1‑251 (MQ=255)                          gAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTtgct                                                                                                                                                                                                         >  1:290356/1‑251 (MQ=255)                           aaTAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTtgctg                                                                                                                                                                                                        >  2:277029/1‑251 (MQ=255)                              aaaCCGGGAAACGGAGCAAATCAGCGTCATTCCGGGTACCTCACCCCCACTGAGAACGTAATCGCCGATTGCCCATTCGGCGGCAATTTCGGTTTGGACCACGCGCTCATCTATGCCTGCGTAGCGACCGCCCACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACCCCCGCTGGATCAGGCTTGCGCCCCTGTGGT‑GACATATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCAt                                                                                                                                                                                                     <  2:568373/251‑1 (MQ=255)                                  cGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAAt                                                                                                                                                                                                 >  1:338133/1‑251 (MQ=255)                                   gggAAACGGATTCAATCAGCGTCATTGCGGGTAACTCCCCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCTCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGTT‑GACAGATAAAACACCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATg                                                                                                                                                                                                <  2:255742/251‑1 (MQ=255)                                        aCGGAGTCAATCAGCGTTATTGCTGGTCACTCACCACCACGGAGAACGTATTCGCCGATTGCCCATTCTTCGTCAATTTCGGTTTGGATAACGCGCTCCTCTATACCTTCGTAGCGACCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATAAAGCTTGCGTCCCTGTGTT‑GACAGATAAATCTCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTc                                                                                                                                                                                           <  2:516605/251‑1 (MQ=255)                                          ggtgcgaaTCAGCGTCATTTCTGGAACCTCTCCGCCCCGGAGACCGTAATCGCCGATTGACTATTCTTCGTCAATTTGGGTTTGTATCACGCTCTCATTTATCCCTTCGTAGCGCCCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGTTTGCGTCCCGGGGTT‑GACAGTTAAATCACCTTTCCCCCTTCTCCCGCCGCGGCTTTTGCTGCATGAATGGCGTccc                                                                                                                                                                                         <  1:187559/245‑1 (MQ=255)                                               cAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa                                                                                                                                                                                      >  2:104074/1‑249 (MQ=255)                                                    ngCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAggg                                                                                                                                                                                  >  2:341128/2‑248 (MQ=255)                                                    ggCGTCATTGCTGGTAACTCCCCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAggg                                                                                                                                                                                  <  1:341128/247‑1 (MQ=255)                                                    ggCGTCATTGCTGGAAACTCCCCCCCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTg                                                                                                                                                                               <  2:609874/250‑1 (MQ=255)                                                     gCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGc                                                                                                                                                                              >  2:123765/1‑251 (MQ=255)                                                       gTCATTGCTGTTAACTCACCACCCCTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc                                                                                                                                                                            <  2:190124/251‑1 (MQ=255)                                                         cATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAcc                                                                                                                                                                           >  2:447507/1‑250 (MQ=255)                                                         cATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAcc                                                                                                                                                                           <  2:290356/250‑1 (MQ=255)                                                          aTTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcat                                                                                                                                                                         <  1:70419/251‑1 (MQ=255)                                                          aTTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcat                                                                                                                                                                         >  2:603852/1‑251 (MQ=255)                                                            tgctggtaccgcccccccacGGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTCCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACCCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatc                                                                                                                                                                        <  2:566297/234‑1 (MQ=255)                                                                    aCTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCTCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATc                                                                                                                                                               >  2:559904/1‑251 (MQ=255)                                                                       caccaccacTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc                                                                                                                                                            >  1:493496/1‑251 (MQ=255)                                                                         ccaccacTCAGAACGTATTCGCCGCTTGCCCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCCTCTCTACCTTCGTCGCGACCGCCCACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACGCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCgg                                                                                                                                                          >  2:268484/1‑251 (MQ=255)                                                                          ccccacGGAGAACGTAACCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCCCGCGCTCATCTATCCCTTCGTAGCGACCGCCCACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGt                                                                                                                                                         <  2:232639/249‑1 (MQ=255)                                                                           cccacTGAGAACGTAATCGCCGATTGCCCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATCCCTTCGTAGCGACCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTc                                                                                                                                                        <  2:640042/250‑1 (MQ=255)                                                                               cgtagaAGGTAATCGCGGATTGACCTTTTGTCGTCAATTTCGGTTTGGATCACGCGCTCTTCTATCCCTTCGTAGCGACCGCACACAAGAATCAATTTTTGATTCGTTGCCAGTGCGCTCACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAAATAAATCCCCTTGGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccg                                                                                                                                                      <  2:607548/246‑1 (MQ=255)                                                                                      cGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTa                                                                                                                                              >  2:185284/1‑250 (MQ=255)                                                                                       gTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTa                                                                                                                                              >  1:193732/1‑249 (MQ=255)                                                                                         aaTCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAgg                                                                                                                                           >  1:733070/1‑250 (MQ=255)                                                                                                    gcccATTCTTCGTCAATTTCGGTTTGGATCGCGCGCTCATCTATCCCTTCGTAGCGCCCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCa                                                                                                                                <  1:123765/248‑1 (MQ=255)                                                                                                      ccATTCTTCGACAATTTGGGTTTGAATCACCCTCTCATCTTTCCCTTCGTAGCGACCGCACCCCAAAATCATTTTTTGATTCGTGGCCATTTCGCTGCCGCCCGCTTAATCAACCTTGCGCCCTTGTGGT‑GACAAAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacgg                                                                                                                             <  2:398259/251‑1 (MQ=255)                                                                                                        aTTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCacggt                                                                                                                            >  2:765980/1‑250 (MQ=255)                                                                                                              cGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTACGGACGATCGTCCACGGTACGGTGc                                                                                                                     >  2:523464/1‑251 (MQ=255)                                                                                                               ttCAATTCCGGTTTGGACCCCGCGCTCATCTATCCCTTCGTAGCTCCCGCACACCAGAATCAATTTTTTATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatc                                                                                                                                                                        <  2:513090/198‑1 (MQ=255)                                                                                                               gTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatc                                                                                                                                                                        >  1:513090/1‑199 (MQ=255)                                                                                                                       ggttttGGACCCCGCGCCCATCTATACTTCCGTAGCGACCCCCCACCAGAATCAATTTTTGTGTCGTTGCGTTTTCGCTGCCCCCCGCTTGTTAAAGTTTGCGTCCTTGTGGT‑GACAGAAAAATCACCTGTGGCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGc                                                                                                            <  2:411678/248‑1 (MQ=255)                                                                                                                               aTCACGCGCTCATCTATCCCTTCGTAGCGACCGCACCCCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTc                                                                                                    <  1:104074/251‑1 (MQ=255)                                                                                                                                            ttatcccTTCGTAGCGCCCGCACACCAGAATCATTTTTTGATTCTTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGTT‑GACAGATAAATCACCTTTGCGCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCa                                                                                        <  1:447507/245‑1 (MQ=255)                                                                                                                                              ataCCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGct                                                                                     >  1:627211/1‑251 (MQ=255)                                                                                                                                                  cgtcggggCGACCGCACCCAAAAATCAATTTTGGATTCGTGGCCGGTCCGCTGACGCCCGCTGGATCAAGCTGGCGTCCCTGTGTG‑GACAGATAAATCACCTTGGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTgg                                                                                 <  2:658455/244‑1 (MQ=255)                                                                                                                                                    tCGTAGCGCCCGCCCACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGa                                                                                <  1:277029/250‑1 (MQ=255)                                                                                                                                                      gTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGc                                                                             >  1:212721/1‑251 (MQ=255)                                                                                                                                                       gggCGACCGCACACCAAAATCAATTGTTGTTTCTTTCCCAGTCCGCTGACCCCCGCTTGATCAGGCTTGCGTCCCTGTGGT‑GACAGATAAATCCCCTTGGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCt                                                                            <  2:509940/249‑1 (MQ=255)                                                                                                                                                        aGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGCGGACTCCAGCTCTGGATGCTc                                                                           >  2:526580/1‑251 (MQ=255)                                                                                                                                                         gCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTc                                                                           >  1:46360/1‑250 (MQ=255)                                                                                                                                                          cGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTc                                                                           >  1:725333/1‑249 (MQ=255)                                                                                                                                                          cGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCCGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTggtttctc                                                                           >  2:416747/1‑243 (MQ=255)                                                                                                                                                                        aaTCAATTTTTGTTTCGTTGCCATTTCGCTGCCCCCCGTTTGATCAAGCTTGCTTCCGTGTGGT‑GACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTaa                                                           <  1:471354/251‑1 (MQ=255)                                                                                                                                                                        aaTCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTa                                                            >  1:494898/1‑250 (MQ=255)                                                                                                                                                                               ttttGATTCGTTGCCAGTTCGCTGACGCCCGTTTGATCAAGCTTGCGCCCCTGTGGT‑GACAGAAAAATCACCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGccc                                                     <  1:406757/250‑1 (MQ=255)                                                                                                                                                                                tttGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCCTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCgg                                                   >  1:340612/1‑251 (MQ=255)                                                                                                                                                                                ttggtttctttgccattccccTGACCCCCGCTTGACCAGGTTTGCGTCCCT‑TGGTTGAAAGAAAAACCCCCTTGGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCg                                                    <  2:627211/230‑1 (MQ=255)                                                                                                                                                                                          tGCCAGTTCGCTAACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCACGGCCAGTTAccc                                         >  1:560362/1‑251 (MQ=255)                                                                                                                                                                                              agtccgcggACGCCCGCTGGATCAAGTTTGCGTCCCTGTGGT‑GACAGATAAATCCCTTTTGCGCCTTCACCCGCCGCGGCTTTGGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCAtttt                                                              <  1:456951/218‑1 (MQ=255)                                                                                                                                                                                              aGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCGCAGCAGGCCAtttt                                                              >  2:456951/1‑226 (MQ=255)                                                                                                                                                                                                 tCGCTGCCCCCCGCTTGATCAACCTTGCTCCCCTGTGTT‑GACAAATAAACCCCCTTGCCCCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAAt                                   <  2:668763/250‑1 (MQ=255)                                                                                                                                                                                                  cGCTGACGCCCGCTTGATCAAGCTGGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCCCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATc                                  <  1:487575/250‑1 (MQ=255)                                                                                                                                                                                                   gCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCgg                                >  1:223666/1‑251 (MQ=255)                                                                                                                                                                                                     ngACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGt                               >  2:343225/2‑250 (MQ=255)                                                                                                                                                                                                       ccGCCCGTTGGATCAAGCTGGCGCCCCGGTGGT‑AACAGAAAAACCCCCTTTGCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAAt                            <  2:46360/250‑1 (MQ=255)                                                                                                                                                                                                         ccccGCGTGATCATGCTGGCTCCCCTGTGGT‑AACAGATAATCCCCCTTTGCGCTTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAAttn                          <  2:340612/250‑2 (MQ=255)                                                                                                                                                                                                             gCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGcgcg                      >  2:676486/1‑251 (MQ=255)                                                                                                                                                                                                               ttGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTgcgtgg                     >  2:688054/1‑247 (MQ=255)                                                                                                                                                                                                                    cAAGCTTGCGTCCCTGTGGT‑AACAGAAAAATCGCCTTTGCCCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACAt                <  1:408077/250‑1 (MQ=255)                                                                                                                                                                                                                          tGCGTCCCTGTGGT‑GACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAgg          <  1:523464/250‑1 (MQ=255)                                                                                                                                                                                                                             gTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAAc      >  1:579271/1‑251 (MQ=255)                                                                                                                                                                                                                                 cTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGc  >  1:529639/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        GAGGCTTCATGTCCCAGTACTCCCGGAATAAACCGGGAAACGGAGTCAATCAGCGTCATTGCTGGTAACTCACCACCACTGAGAACGTAATCGCCGATTGACCATTCTTCGTCAATTTCGGTTTGGATCACGCGCTCATCTATACCTTCGTAGCGACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGT‑GACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGC  >  NC_000913/2744846‑2745320

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