Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,116,975	C→T	R57* (CGA→TGA)	trhO →	tRNA U34 hydroxylase TrhO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,116,975	0	C	T	97.7%	140.8 / NA	44	R57* (CGA→TGA)	trhO	tRNA U34 hydroxylase TrhO
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/18);  minor base G (1/0);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCG  >  NC_000913/1116737‑1117220                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      tGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTgg                                                                                                                                                                                                                                           >  1:17539/1‑251 (MQ=255)       ttATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATg                                                                                                                                                                                                                                     >  2:7843/1‑251 (MQ=255)          tCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAgg                                                                                                                                                                                                                                  >  2:203814/1‑251 (MQ=255)               aaTTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACGACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCAtt                                                                                                                                                                                                                             >  1:245410/1‑251 (MQ=255)                                      ttttACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACGACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTg                                                                                                                                                                                                       <  2:245410/250‑1 (MQ=255)                                       tttACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCg                                                                                                                                                                                                     >  1:211228/1‑251 (MQ=255)                                          aCTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAg                                                                                                                                                                                                   <  2:211228/250‑1 (MQ=255)                                             ttCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAAt                                                                                                                                                                                               >  1:253243/1‑251 (MQ=255)                                                ccATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATGTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGtt                                                                                                                                                                                            <  1:341793/251‑1 (MQ=255)                                                 cATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTg                                                                                                                                                                                           >  2:147416/1‑251 (MQ=255)                                                    ctcGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGaa                                                                                                                                                                                         >  1:394543/1‑250 (MQ=255)                                                         cggAACAGTACACCATGCCAGTGTTACACAACCGCCTTTCCAACGACACGCTAAAATCCACAATGTTGGCTGAGAACGAACCGCGAACCACCATTTCGCTTTACAAGTATTTCCACATAGCCGATCCTAAGGCGACCCGTGCCGCTTTCTATCAGCTGTTTAACGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACATGCGAGCACTGTTGAAACAttt                                                                                                                                                                                   >  2:382195/2‑251 (MQ=255)                                                                gTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACAATTCgcgcgc                                                                                                                                                                            >  2:620187/1‑251 (MQ=255)                                                                       aTGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGctct                                                                                                                                                                      >  2:325182/1‑250 (MQ=255)                                                                          ccAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGc                                                                                                                                                                  >  2:39125/1‑251 (MQ=255)                                                                               gtTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTc                                                                                                                                                              <  1:324897/250‑1 (MQ=255)                                                                                   cacaACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGAccc                                                                                                                                                         <  1:472628/251‑1 (MQ=255)                                                                                                 cAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAgtgt                                                                                                                                                                                                                                                 <  1:13615/149‑1 (MQ=255)                                                                                                 cAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAgtgt                                                                                                                                                                                                                                                 >  2:13615/1‑149 (MQ=255)                                                                                                            tAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGGTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCAGGCTTTAGAGGGTTTACGCCTGAat                                                                                                                                >  1:632251/1‑251 (MQ=255)                                                                                                               aaGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACACGTCTTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGATTACGCCTGAATATc                                                                                                                             >  1:77754/1‑251 (MQ=255)                                                                                                                aGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATcg                                                                                                                            <  2:295134/251‑1 (MQ=255)                                                                                                                 gCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATcgc                                                                                                                           <  1:180630/251‑1 (MQ=255)                                                                                                                    aaaaTGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGtt                                                                                                                        >  1:187233/1‑251 (MQ=255)                                                                                                                          ttGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGgatg                                                                                                                   >  1:482663/1‑250 (MQ=255)                                                                                                                                 agagCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACggg                                                                                                           <  1:104512/251‑1 (MQ=255)                                                                                                                                   agCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGaa                                                                                                         <  1:203814/251‑1 (MQ=255)                                                                                                                                               accaccATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGCATATCGCGTTGGCTGATGACGGGCAATCCTTCTGGGt                                                                                             >  2:453775/1‑251 (MQ=255)                                                                                                                                                                      tttCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATc                                                                      >  1:13246/1‑251 (MQ=255)                                                                                                                                                                         ccACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGc                                                                    <  1:335430/250‑1 (MQ=255)                                                                                                                                                                           acaTCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGCAtt                                                                 >  1:429869/1‑251 (MQ=255)                                                                                                                                                                               cGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGCAttgtt                                                              <  1:39125/250‑1 (MQ=255)                                                                                                                                                                                gCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGt                                                                                                                                            >  1:112700/1‑171 (MQ=255)                                                                                                                                                                                gCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGt                                                                                                                                            <  2:112700/171‑1 (MQ=255)                                                                                                                                                                                      ccTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGCATTGTTGCCGACgg                                                      <  1:249609/251‑1 (MQ=255)                                                                                                                                                                                         aaGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGTATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGATCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGCGAAATCCTTCTTGGTACTGCGCATGAAGGTACGCGATCGCATTGTTGCCGACGGt                                                     <  1:382195/249‑1 (MQ=255)                                                                                                                                                                                              gccccGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTGTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACg                                              <  2:459143/249‑1 (MQ=255)                                                                                                                                                                                                  cGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATcc                                          <  1:644392/251‑1 (MQ=255)                                                                                                                                                                                                   gTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCATTCTGGGTACTGCGCATGAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCt                                         >  1:301348/1‑251 (MQ=255)                                                                                                                                                                                                                        ttACCGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGCTCGCATTGTTGCCGACGGTATTGACGATCTTCACTTTGATGCCAGCAATGtt                    >  2:172234/1‑251 (MQ=255)                                                                                                                                                                                                                           ccGCGCTGAATGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGt                 <  1:157621/251‑1 (MQ=255)                                                                                                                                                                                                                                     tGTTTTTGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGCTGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGc       >  2:645578/1‑251 (MQ=255)                                                                                                                                                                                                                                     tGTTTTTGGGGGAGTGTATCTGGAGCATGAGGGCATTACCGCGCAAATCAGCGAACCAGAGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCATTCGACCCGGCTTTAAAGGGTTTACGCCTGCATATCGCGTTGGATGATGCCGGGAAATCCTTCTGGGTACTGCGCATGCAGGGACGCGATCGCATATTTGCCGACGGTATTGACGACCCTCACTTTGATGCCAGCAATGTTGGTGAGAATCTGc       >  2:375040/1‑251 (MQ=255)                                                                                                                                                                                                                                          ttGGGTGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAgcg  >  1:138577/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      TGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGGCATTAACGCGCAAATCAGCGTACCTGCGAGCAATGTTGAAACATTTCGCGCGCAGCTCTATGCCTTCGACCCGGCTTTAGAGGGTTTACGCCTGAATATCGCGTTGGATGATGACGGGAAATCCTTCTGGGTACTGCGCATGAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCG  >  NC_000913/1116737‑1117220

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