Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,883	+C	coding (667/687 nt)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,883	1	.	C	95.2%	204.2 / 1.9	63	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (42/18);  total (44/19)
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.75e-01

CCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑C‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGC  >  NC_000913/3815655‑3816130                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAggg                                                                                                                                                                                                                                    >  2:140359/1‑249 (MQ=255)      cGTTGAAGAAATCGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGt                                                                                                                                                                                                                                >  1:69836/1‑250 (MQ=255)       gTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAggtgcgt                                                                                                                                                                                                                                >  2:289179/1‑243 (MQ=255)          gaagaaATAGGGGCTTTCGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGctc                                                                                                                                                                                                                           <  2:201895/251‑1 (MQ=255)              aaaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTc                                                                                                                                                                                                                        >  2:157082/1‑250 (MQ=255)               aaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCAt                                                                                                                                                                                                                      <  2:13309/251‑1 (MQ=255)                  aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGCGTAGCCCTTCATggg                                                                                                                                                                                                                   >  2:368955/1‑251 (MQ=255)                  aGGGGCTTTTGCGCACGGATGTCAGCGTAAAATCGCCAAACTTTAATCCCTGATTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCACCAGTCGCCTTAAAAATCAGTTGGCCAGCGCCGCCTACTGCGTCGCTACAATGGATTCGCGTG‑CCC‑C‑TCGGGCCAGCTCCAATACGGTGAGTAGTTCTTCaggg                                                                                                                                                                                                                    <  2:65728/250‑4 (MQ=255)                              gcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGGGAACGGCTCGc                                                                                                                                                                                                       >  2:208982/1‑251 (MQ=255)                                 cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                                                                                                                                                    <  1:434253/251‑1 (MQ=255)                                     tttCAGCGTAAACTCGCCAAACTTTAACACCTGCCTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                                                                                                                                                                                                >  2:436628/1‑251 (MQ=255)                                        cAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGt                                                                                                                                                                                             <  1:140359/251‑1 (MQ=255)                                                   cGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGACCCCTGCACtt                                                                                                                                                                                  >  2:27281/1‑251 (MQ=255)                                                    gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATTACCC‑C‑TCGTGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                                                                                                                 >  1:301168/1‑251 (MQ=255)                                                              aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcg                                                                                                                                                                        <  2:495371/250‑1 (MQ=255)                                                                      cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGATACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCACCACGGTGAGTAGCGCTTCATGGGTGCACGGCTCGCCTTCTGCCGTCACCTGCACTTCAATGATTCGCCCGTCttt                                                                                                                                                               >  2:58884/1‑250 (MQ=255)                                                                         gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGACCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                                                                            >  1:283450/1‑251 (MQ=255)                                                                            aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                                                                         >  1:44557/1‑251 (MQ=255)                                                                               cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                                                                      >  1:128520/1‑251 (MQ=255)                                                                                   aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGGGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                                                                                                                  >  2:371363/1‑251 (MQ=255)                                                                                     ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                                                                                >  1:84153/1‑251 (MQ=255)                                                                                            aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                                                                         >  1:489228/1‑251 (MQ=255)                                                                                                ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTACTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct                                                                                                                                     <  1:406507/251‑1 (MQ=255)                                                                                                     tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                                                                                >  1:464361/1‑251 (MQ=255)                                                                                                        tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGGCGGTCTCTGCGGct                                                                                                                             >  1:225940/1‑250 (MQ=255)                                                                                                               tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAACAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCtt                                                                                                                      >  2:188173/1‑251 (MQ=255)                                                                                                                cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCtt                                                                                                                      <  1:368955/250‑1 (MQ=255)                                                                                                                cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCtt                                                                                                                      <  1:208982/250‑1 (MQ=255)                                                                                                                 aTGCCTTCGCGACTCATCTTACTTTTCGACAGACGAAAAAAAGGCGGCTCACCAGTCGCCTTAAAAAACAGTTTGCCAGCGCCGCCATCTGCGTCCCTTCAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGACAACGCGGTGAGTAGCTCTTCATGGGTAAACGGCTCGCCTTCGGACGTCGCCTGCTCTTCAACTATGCGGGCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCtgca                                                                                                                    <  2:438828/251‑4 (MQ=255)                                                                                                                   gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCaa                                                                                                                   <  1:371363/250‑1 (MQ=255)                                                                                                                               cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTcc                                                                                                           >  2:185227/1‑246 (MQ=255)                                                                                                                                    tACTTTTCTACAGACACAAAAAAGGCGACTCATCAGTCGACATCAAAATCAGTTTGCCAGCACCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTTAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTTCCCTGCACTTCAATGATGCGCCCGTCTTCTGACATCACTAAGTTCCTGTCGGTCGCTGACACTGCATCTTCAACGTATTCCCGATCGCaa                                                                                                  >  2:483022/1‑250 (MQ=255)                                                                                                                                     aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaat                                                                                                >  1:191233/1‑250 (MQ=255)                                                                                                                                          tCTACAGACAAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaacttgctt                                                                                           >  2:60309/1‑244 (MQ=255)                                                                                                                                                 aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGtt                                                                                    >  1:430392/1‑251 (MQ=255)                                                                                                                                                  cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGttt                                                                                   >  1:302322/1‑250 (MQ=255)                                                                                                                                                                 tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAccga                                                                     >  2:316056/1‑247 (MQ=255)                                                                                                                                                                  catccttCGCCTTAAAAACCATTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGGCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGa                                                                     <  2:389351/243‑1 (MQ=255)                                                                                                                                                                  catcaGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTATGCGTCGCTACAATGGATTCG‑ATT‑CCCAC‑TCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACAGACAGaa                                                                    >  2:72285/1‑250 (MQ=255)                                                                                                                                                                     caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTg                                                                >  1:268191/1‑251 (MQ=255)                                                                                                                                                                     caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTg                                                                <  1:185227/251‑1 (MQ=255)                                                                                                                                                                      aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGc                                                               >  1:137216/1‑251 (MQ=255)                                                                                                                                                                          gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCt                                                           >  1:234598/1‑251 (MQ=255)                                                                                                                                                                            cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTAc                                                         >  1:46578/1‑251 (MQ=255)                                                                                                                                                                             ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGGGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGTTTACAATACCGACAGAAACTGCGGCTAcc                                                        >  2:9659/1‑251 (MQ=255)                                                                                                                                                                              tAAAAATCATTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGCGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTTCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTAc                                                         <  2:347002/249‑1 (MQ=255)                                                                                                                                                                               aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAACTTGCTTCACCGTTTACAATACCGACCGAACCTGCGGCTACCAt                                                      >  2:364281/1‑251 (MQ=255)                                                                                                                                                                               aaaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCa                                                       <  1:157082/250‑1 (MQ=255)                                                                                                                                                                                 aaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAACTTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATcc                                                    >  1:407459/1‑251 (MQ=255)                                                                                                                                                                                     cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGTGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTGCGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCtt                                                 <  1:203821/250‑1 (MQ=255)                                                                                                                                                                                           gCCAGCGCCGCATTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACAGACAGAAACTGCGGCTACCATCCCTTTAATCgg                                          >  2:279498/1‑251 (MQ=255)                                                                                                                                                                                                cgccgcCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTgg                                     >  1:43670/1‑251 (MQ=255)                                                                                                                                                                                                    cccTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGTCCAGCGCCAACAAGGTGAGTATCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGGCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCGTTCATCGGGTTGGtttt                                 <  2:157601/250‑1 (MQ=255)                                                                                                                                                                                                            cGTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCAGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCCTGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGc                         >  2:408147/1‑251 (MQ=255)                                                                                                                                                                                                             gTCGCTACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGcc                        <  2:268191/251‑1 (MQ=255)                                                                                                                                                                                                                  tACAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGtttt                   >  1:295221/1‑251 (MQ=255)                                                                                                                                                                                                                   aCAATGGATTCG‑GTT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGCGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTNTCATCGGGTTGGTTTTCAGCTTGCCGTTTTc                  >  2:219550/1‑251 (MQ=255)                                                                                                                                                                                                                   aCAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTc                  >  1:136962/1‑251 (MQ=255)                                                                                                                                                                                                                   aCAATGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGATTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTc                  >  1:175914/1‑251 (MQ=255)                                                                                                                                                                                                                     aaTGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTccc                >  2:176517/1‑250 (MQ=255)                                                                                                                                                                                                                      aTGGATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTACACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTccac               >  2:29032/1‑251 (MQ=255)                                                                                                                                                                                                                        ggATTCG‑ATT‑CCC‑CCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGGCCCCTGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTTGTTTTCAGCTTGCCGTTTTccacca             <  2:191233/251‑1 (MQ=255)                                                                                                                                                                                                                           ttcg‑att‑cCC‑CCTCGGGCCAGCGCCAACCAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCCTCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCCCCAGCt          >  2:490643/1‑251 (MQ=255)                                                                                                                                                                                                                                   cccc‑c‑TCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGt     >  2:464810/2‑249 (MQ=255)                                                                                                                                                                                                                                    ccc‑c‑TCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgc  >  1:201499/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        CCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCG‑ATT‑CCC‑C‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGC  >  NC_000913/3815655‑3816130

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