Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,407,045	C→T	A131V (GCT→GTT)	rnr →	RNase R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,407,045	0	C	T	96.5%	185.7 / NA	59	A131V (GCT→GTT)	rnr	RNase R
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/22);  minor base G (0/2);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCT‑TTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCGGTGGGTAAAATCGTCGA  >  NC_000913/4406808‑4407288                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      ttGAAGGCTTGCATCGCCGCATGCGCGCGATGGAGCCCGATGGTCAACTGGTCTTCCCTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGCAAGGCACCGTTATTGGCCACCGGGATGGCTACGGCTTTCTGCGGGTTGTAGGTCGTAAAGATGATTTGTATCTCTCCCGCGAGCAGATGAATACCTGCATTCATGGCGATCCGGTGCTGGCTCAGCCGCTGGGTG‑TGTACCGTAAAGgt                                                                                                                                                                                                                                         >  2:85139/1‑251 (MQ=255)  tGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGgtc                                                                                                                                                                                                                                        >  1:402654/1‑251 (MQ=255)              tCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGAAAAGGTCGTCGTGAAgcgc                                                                                                                                                                                                                            >  1:219019/1‑251 (MQ=255)                   gccTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTAtt                                                                                                                                                                                                                       <  1:135078/251‑1 (MQ=255)                    ccTGCGCCCGATGGAGCGCGATGGACACCTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGCCCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAACGATGATTTGTATCTCTCCAGCGAGCAGATGAAACCCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACAGTAAAGGTagtcgtaccacgcgtata                                                                                                                                                                                                                       >  2:462496/1‑238 (MQ=255)                          cgcgATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTCCGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCcgcg                                                                                                                                                                                                                >  2:384309/1‑251 (MQ=255)                          cgcgATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAACCCTGCATTAATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGACGTCGTGAAGCGCGTATTGTCcgcg                                                                                                                                                                                                                >  2:371105/1‑251 (MQ=255)                                 gAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGt                                                                                                                                                                                                         <  1:416973/251‑1 (MQ=255)                                   gcgcGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGc                                                                                                                                                                                                       <  1:244295/251‑1 (MQ=255)                                    cgcgATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGTGTGTAAAGATGATTTGTATCTCTCCAGCGAGCAGAGGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGcc                                                                                                                                                                                                      <  1:446622/251‑1 (MQ=255)                                          ggTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAc                                                                                                                                                                                                >  2:171982/1‑251 (MQ=255)                                           gTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAcc                                                                                                                                                                                               >  1:437050/1‑251 (MQ=255)                                                 tGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAccagcc                                                                                                                                                                                           >  1:156955/1‑249 (MQ=255)                                                     cTTCACTCGTCGTCAGTGCTATGCGCTGCCGGACCGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGATGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTAATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGAACTGGTGCCAAAAACCAGCCAGAttg                                                                                                                                                                                     >  1:381309/1‑251 (MQ=255)                                                               cgtcAGTGCTATGCGCTGCCGTAACGCCGCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTGCGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGAGAGCAGATGAAATCCTGCATTCATGGCGAGCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTAGTGGCCGCGTACTGGTTCCAAAAACCAGCCAGAATGTTGGTCGCTa                                                                                                                                                                           <  2:371394/251‑1 (MQ=255)                                                                         aTGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAccagcc                                                                                                                                                                                           <  1:231350/225‑1 (MQ=255)                                                                         aTGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAccagcc                                                                                                                                                                                           >  2:231350/1‑225 (MQ=255)                                                                          tGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGGGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCgc                                                                                                                                                                  >  1:25818/1‑248 (MQ=255)                                                                             gcTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCgg                                                                                                                                                             >  1:161918/1‑251 (MQ=255)                                                                                    gAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCgg                                                                                                                                                      >  1:218123/1‑251 (MQ=255)                                                                                       cGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCtt                                                                                                                                                   <  1:361935/251‑1 (MQ=255)                                                                                        gCCTCGCCCTGGTGAAAGGTCACGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGTGCGTAAAGATGATTGGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCTCGTATTGTCCGCGTACTGGTGCCAAAATCCAGCCAGATGGTTGGTCGCTACTTTACCGAAGCGggctttggcttt                                                                                                                                                  <  2:2845/251‑1 (MQ=255)                                                                                           ccGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTtgtg                                                                                                                                               <  2:156955/250‑1 (MQ=255)                                                                                            cGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGATGCGGGCGTCGGTTTtgtg                                                                                                                                               <  2:405607/250‑1 (MQ=255)                                                                                                cTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTcc                                                                                                                                          >  2:189547/1‑251 (MQ=255)                                                                                                  ggTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTg                                                                                                                                        >  1:200449/1‑251 (MQ=255)                                                                                                     gAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTgacg                                                                                                                                     <  2:187115/251‑1 (MQ=255)                                                                                                     gAAAGGTACAGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTgacg                                                                                                                                     >  1:388953/1‑251 (MQ=255)                                                                                                          gTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGc                                                                                                                                <  1:388968/251‑1 (MQ=255)                                                                                                            aCCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCg                                                                                                                              >  1:452881/1‑251 (MQ=255)                                                                                                             ccGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGt                                                                                                                             >  1:335837/1‑251 (MQ=255)                                                                                                                 tATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGGCGGCTTTGTGGTTCCTGATGTCAGCCGCCTgg                                                                                                                         >  2:362926/1‑250 (MQ=255)                                                                                                                     ggCCACCGTGATGGCTACGGCTTTCTGCGGGTTGACGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAcca                                                                                                                                                                                              >  2:104141/1‑178 (MQ=255)                                                                                                                     ggCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAAcca                                                                                                                                                                                              <  1:104141/178‑1 (MQ=255)                                                                                                                     ggCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCt                                                                                                                      <  2:80100/250‑1 (MQ=255)                                                                                                                       ccACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCg                                                                                                                   >  2:112080/1‑251 (MQ=255)                                                                                                                        cACCGCGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCAAGCAGAAGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTCGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGAGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCATGACGACAGCCGTCTGAGCTTCg                                                                                                                   >  2:305699/1‑250 (MQ=255)                                                                                                                               aTGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTGTACCGAAGCGGGCGTCGGCTGTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTa                                                                                                           <  2:301135/251‑1 (MQ=255)                                                                                                                                  gCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATc                                                                                                        >  1:153037/1‑251 (MQ=255)                                                                                                                                       ggCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTACATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGGTG‑ACCGTAAAGGTCGTCGTGAAGCGCGGATTGTCCGCGTTCTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCTTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATcccgcc                                                                                                   <  2:405450/251‑1 (MQ=255)                                                                                                                                               gCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCgatca                                                                                             <  1:207190/249‑1 (MQ=255)                                                                                                                                                cGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGCTGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGACGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCccgatccgat                                                                                          >  2:156326/1‑251 (MQ=255)                                                                                                                                                      gAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATggg                                                                                    >  1:287021/1‑251 (MQ=255)                                                                                                                                                            cGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGcgcgcg                                                                              >  1:430137/1‑251 (MQ=255)                                                                                                                                                                    tgatTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCt                                                                      >  1:108969/1‑251 (MQ=255)                                                                                                                                                                     gatTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCt‑t                                                                    >  1:117244/1‑251 (MQ=255)                                                                                                                                                                     gatTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGGCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCt                                                                      <  2:388953/250‑1 (MQ=255)                                                                                                                                                                                   cagggAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGGTCAGCCGCGGGGGGTTG‑AGCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCT‑TTGTGGTCGTAGTCg                                                      <  2:129058/247‑1 (MQ=255)                                                                                                                                                                                   cAGCGAGCAGATGAAAAACTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGCTGGGCT‑TTGTGGTCGTAgacg                                                      >  2:20072/1‑248 (MQ=255)                                                                                                                                                                                        aGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCTCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCT‑TTGTGGTCGTAGTCGAACTg                                                 <  1:326093/251‑1 (MQ=255)                                                                                                                                                                                               gcaaaCCTGCAGTCATGGCGATCGGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGCGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCAAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTGTGAGGTTCCTGACGACAGCCGTCTGAGCTTCGTTAGCTTAATACAGCCAGGTCAGATCATGGGCGCGCGGATGGGCTGTTG‑GGTCGTAGTCGCAATTACTCAGc                                          >  2:16056/3‑251 (MQ=255)                                                                                                                                                                                                aaaaCCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCT‑TTGTGGTCGTAGTCGAACTTACTCAGCg                                         <  1:136199/251‑1 (MQ=255)                                                                                                                                                                                                   aCCTGCATGCATGGCGATCAGGTGCGGGCGCAGCCGCTGGTTTGTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTAGTGGCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACGTTACCGAAGCGGGCGTCGGCGTTGTGGTTCCTGTCGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCT‑TTGTGGTCGGAGTCGAACTGACTCAGCGTcc                                      <  2:218123/251‑1 (MQ=255)                                                                                                                                                                                                        cATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGCTCAGATCATGGGCGCGCGGTTGGGCT‑TTGTGGTCGTAGTCGAACTGACGCAGCGTCCGACTc                                 >  2:335987/1‑251 (MQ=255)                                                                                                                                                                                                            cATGGCGATCAGGTGCTGGCTCCGCCGCTGGGTGTTG‑ACCGTAAAGGTCGGCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGAGCATGGGCGCGCGGATGGGCT‑TGGTGGTCGTAGTCGAACTGACTCAGCGTCCTACTcgcc                              <  2:75131/250‑1 (MQ=255)                                                                                                                                                                                                                             ggCTCAGCCGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCT‑TTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCGGTGGGTa            <  2:153037/251‑1 (MQ=255)                                                                                                                                                                                                                                    ccGCTGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCT‑TTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCGGTGGGTAAAAtcgg     >  2:180474/1‑250 (MQ=255)                                                                                                                                                                                                                                       cTGGGTGTTG‑ACCATAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCAAAAACCCAGCCAGCTTGGTGGTGGCTACTTTACCGACGCGGACGTCGGCTTTGTGGTTCCTGACGACAGCGGTCTGAGCTTCGATCTCTTAATCCCGCACGTTCAGATCATGGGCGGGCGGACGGGCT‑TTGTGGTCGTAGTCGGACTGGCTCAGCGCCCGCCTCGCTGCCCCTAAGCGGTGGGTatcatctgcgg  >  2:128715/1‑241 (MQ=255)                                                                                                                                                                                                                                        tGGGTGTTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGCTTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGACCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGCTCAGATCATGGGCGCGAGGATGGGCT‑TTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCGGTGGGTAAAAtcgtcgc  >  1:455101/1‑249 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      TTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCACCGTGATGGCTACGGCTTTCTGCGGGTTGAAGGGCGTAAAGATGATTTGTATCTCTCCAGCGAGCAGATGAAAACCTGCATTCATGGCGATCAGGTGCTGGCTCAGCCGCTGGGTGCTG‑ACCGTAAAGGTCGTCGTGAAGCGCGTATTGTCCGCGTACTGGTGCCAAAAACCAGCCAGATTGTTGGTCGCTACTTTACCGAAGCGGGCGTCGGCTTTGTGGTTCCTGACGACAGCCGTCTGAGCTTCGATATCTTAATCCCGCCCGATCAGATCATGGGCGCGCGGATGGGCT‑TTGTGGTCGTAGTCGAACTGACTCAGCGTCCGACTCGCCGCACCAAAGCGGTGGGTAAAATCGTCGA  >  NC_000913/4406808‑4407288

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