Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012971	875,191	0	G	T	62.5%	18.2 / 31.3	32	pseudogene (4/1116 nt)	ECD_RS04405	replication endonuclease
Reads supporting (aligned to +/- strand):  ref base G (6/6);  new base T (6/14);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGCCTGTCGTGCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTAGCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATCCCTGCGTGATAACGCTGAATACGTTAATGCCT  >  NC_012971/874897‑875488                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                          taaGCCTGTCGTGCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAATAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggc                                                                                                                                                                                                                                                                                                      <  2:184256/299‑3 (MQ=255) gaAGCCTGTCGTGCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGNGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggc                                                                                                                                                                                                                                                                                                      <  1:573196/300‑3 (MQ=255) gaAGCCTGTCGTGCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggc                                                                                                                                                                                                                                                                                                      <  1:340428/300‑3 (MQ=255)  aaGCCTGTCGTGCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggc                                                                                                                                                                                                                                                                                                      <  2:90479/299‑3 (MQ=255)   agcttGTCGTGCGTTCGCAGTGCGCCACAATCTTGTGCGACTGGCAGACCGTGCTGGCATGAATGGGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTAGTAATCGCACCAGAAATCTTGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCGGCCATGCGTACTGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTTTTGCTGCAGGTGCGGTATCTggcga                                                                                                                                                                                                                                                                                                    <  2:348128/296‑5 (MQ=255)   aGCCTGTCGTGCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggaga                                                                                                                                                                                                                                                                                                    >  2:497850/1‑296 (MQ=255)   aGCCTGTCGTGCATTCGCATTGCGCCACAACCTGGTGCAACGGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCATAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAGCGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcga                                                                                                                                                                                                                                                                                                    <  2:30430/300‑5 (MQ=255)   aGCCTGTCGTGCATTCGCAGTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCGGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAATAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcga                                                                                                                                                                                                                                                                                                    <  2:348111/300‑5 (MQ=255)      cTGTCGTGCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcgatgt                                                                                                                                                                                                                                                                                                 >  2:534340/1‑293 (MQ=255)            tGCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGAGTGTTGCTGCAGGTGCGGTATCTggcgatgtaaaaac                                                                                                                                                                                                                                                                                           >  2:732515/1‑287 (MQ=255)            tGCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTACGGAACAAGTTGAATCCAGCGCAACCTAATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCCGAATGAAGATTCAAAGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGCTTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGAGCGTGTTGCTGCAGGTGCGGTATCTggcgatgta                                                                                                                                                                                                                                                                                                >  2:732551/1‑287 (MQ=255)             gCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAATGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcgatgtaaaaacc                                                                                                                                                                                                                                                                                          <  2:281522/300‑15 (MQ=255)             gCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTTGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcgatgtaaaaacc                                                                                                                                                                                                                                                                                          <  1:452134/300‑15 (MQ=255)             gCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcgatgtaaaaacc                                                                                                                                                                                                                                                                                          <  1:106490/300‑15 (MQ=255)                   gCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcgatgtaaaaaccagtgc                                                                                                                                                                                                                                                                                     <  1:307369/299‑20 (MQ=255)                       tGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGNTTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcgatgtaaaaaccagtgcaggtc                                                                                                                                                                                                                                                                                <  2:581748/300‑25 (MQ=255)                       tGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcgatgtaaaaaccagtgcaggtc                                                                                                                                                                                                                                                                                <  2:467792/300‑25 (MQ=255)                        gcgcCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATTAGGTGGCAAAATAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcgatgtaaaaaccagtgcaggtc                                                                                                                                                                                                                                                                                <  2:457874/299‑25 (MQ=255)                        gcgcCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcgatgtaaaaaccagtgcaggtc                                                                                                                                                                                                                                                                                >  1:301687/1‑275 (MQ=255)                           ccACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACNGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTATCTggcgatgtaaaaacaagtgcaggtcgtc                                                                                                                                                                                                                                                                             >  2:635433/1‑272 (MQ=255)                                                                                                                                                                                                                                                                     taaagccgagttgcgctggtatggcgtgcgggTAGCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATc                                  <  2:499858/270‑1 (MQ=255)                                                                                                                                                                                                                                                                     taaagccgagttgcgctggtatggcgtgcgggTAGCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATc                                  <  1:201695/270‑1 (MQ=255)                                                                                                                                                                                                                                                                            gagttgcgctggtatggcgtgcgggTAGCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAGACGGGTAAATCCCTGCg                            <  2:524572/276‑1 (MQ=255)                                                                                                                                                                                                                                                                             agtttcgctggtatggcgtgcgggTAGCTTAGCCGCATCATGACGGTACTGTGCACTGGCAGCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATCCCGGCg                            <  2:524573/276‑1 (MQ=255)                                                                                                                                                                                                                                                                                    ctggtatggcgtgcgggTAGCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGAAGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATCCATGCGTGATAACGc                   >  2:749757/16‑300 (MQ=255)                                                                                                                                                                                                                                                                                      ggtatggcgtgcgggTAGCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACNTCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATCCCTGCGTGATAACGCt                  >  1:168665/14‑299 (MQ=255)                                                                                                                                                                                                                                                                                      ggtatggcgtgcgggTAGCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCATGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATCCCTGCGTAATAACGCt                  >  2:8779/14‑299 (MQ=255)                                                                                                                                                                                                                                                                                      ggtatggcgtgcgggTAGCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTAATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATCCCTGCGTGATAACGCt                  >  1:425029/14‑299 (MQ=255)                                                                                                                                                                                                                                                                                      ggtatggcgtgcgggTAGCTGAGCCACATCATGACGGTACTGTGCACTGGCATCTCATATGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAACCGGGTAAATCCATGCGTGATAACGCt                  >  1:757748/14‑299 (MQ=255)                                                                                                                                                                                                                                                                                           ggcgtgcgggTAGCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTAATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGAGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCAAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTaaa                                    >  2:320027/9‑276 (MQ=255)                                                                                                                                                                                                                                                                                                   gttAGCTGAGCCGCATCATGACGGTACTGTACACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATCCCTGCGTGATAACGCTGAATACGTTAATGc    <  2:94030/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                                   ggTAGCTGATCCGCATCATGACGGTACTATGTACCGGCATCTCATTTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGATGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGGCTGGTTGGGGAGATCAGCAAGGAAACGGGTAAATCCCTGCGTGATAACGCTGAATACGTTAATGc    <  2:94375/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       gCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATCCCTGCGTGATAACGCTGAATACGTTAATGCCt  >  2:731403/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                          GAAGCCTGTCGTGCATTCGCATTGCGCCACAACCTGGTGCAACTGGCAGAACGTGCTGGCATGAATGTGCAGATTCTGCGGAACAAGTTGAACCCAGCGCAACCTCATTTATTAACCGCACCAGAAATCTGGCTGCTTACCGATCTGACTGAAGATTCAACGCTGGTAGATGGTTTTCTGGCTCAGATTCACTGCCTGCCATGCGTACCGATTAATGAGGTGGCAAAAGAGAAACTGCCACATTACGTCATGAGTGCAACCGCAGAGATCGGGCGTGTTGCTGCAGGTGCGGTAGCTGAGCCGCATCATGACGGTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGCGTAAGTTTGCCATCCGTGAAGACCGCGAGGAGCTGGGCAATAACACGGGACCACGCTTTAAGTCTGAGCTGATAAACCCGCGCAAAGGTACGCCAACAAGCTATATCGCGAAATACATCAGTAAGAACATTGACGGGCGTGGTCTGGCTGGGGAGATCAGCAAGGAAACGGGTAAATCCCTGCGTGATAACGCTGAATACGTTAATGCCT  >  NC_012971/874897‑875488

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