Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,585,517	G→A	A436V (GCC→GTC)	ggt ←	glutathione hydrolase proenzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,585,517	0	G	A	98.2%	180.8 / NA	55	A436V (GCC→GTC)	ggt	glutathione hydrolase proenzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/25);  minor base C (1/0);  total (29/26)
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.97e-01

CACTGATGGTGGAAACGCGGCGCATTGGTCGCTTCGGCGACGTTCAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGATTTGGCATAGGCTTT  >  NC_000913/3585271‑3585755                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               cACTGATGGTGTAAACGCGGCGAATTGGTCGCTTCGGCGACGTTCAAGCCATAATCGAGGCTATTCACCTCCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCTGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGga                                                                                                                                                                                                                                            <  1:76744/251‑1 (MQ=255)                 gcggcgCATTGGTCGCTTCGGCGACGTTCAAGCCATAGTCGATGCTATTCACCACCATTTGCAGCACTGCAGTGCTGATCAGGCTACCGCCTGGGCTACCGTTAACCAGCCAGGTATTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCAGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTCCGGTACGCCGGGTTTGACGGAGAACTCATCCATCTgg                                                                                                                                                                                                                            >  1:123309/1‑251 (MQ=255)                           ggTCGCTTCTGCGACGTTTAAGCCATAATCGATGCTATTCACCACCATTTGCAGCTCTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGTTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCTGGCGTTTGTTCGGCCCGACGGCGTTGGCATCCCCGCCCACCAGCCCGTAAACGTTCTGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCa                                                                                                                                                                                                                  <  1:184820/251‑1 (MQ=255)                               gCTTCGGCGACGTTCAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAAt                                                                                                                                                                                                              >  1:304182/1‑251 (MQ=255)                                          gTTCAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCGGAATACCGCTCTCGc                                                                                                                                                                                                   >  2:342221/1‑251 (MQ=255)                                           ttCAAGCCATAATCGCTGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGCCGGAGAAATCATCCCTCTGGTTATTAAGCAGAATACCGCTCTCGcc                                                                                                                                                                                                  >  2:290033/1‑251 (MQ=255)                                            tCAAGACATAATCGCTGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGCCGGCGTTGGCACCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCAGCTGGATATTAAGCAGAATACCGATCTCGccc                                                                                                                                                                                                 >  2:78178/1‑251 (MQ=255)                                                         gATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCt                                                                                                                                                                                                                              >  2:83920/1‑209 (MQ=255)                                                         gATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCt                                                                                                                                                                                                                              <  1:83920/209‑1 (MQ=255)                                                           tGCTATTCACCACCATTTGCAGCACTGAAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGtc                                                                                                                                                                                  >  2:160012/1‑250 (MQ=255)                                                               aTTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCg                                                                                                                                                                               >  1:340252/1‑250 (MQ=255)                                                                      accaTTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAggtggg                                                                                                                                                                       >  1:364644/1‑250 (MQ=255)                                                                        cattggaagcacggttgggatgatCCGGCTGCCGCCGGGGCTACCGGTAACCAGCCATGTTTTACCGGCTTGCGCCACAAGGGTCGGCGACATCGACGACAGCGGGCGTTTGTGCGGCCCGTCGTCGTTGGCATCACCGCCCACCAGCCCGGTAACGTTCGGTGCGCCCGGGATGTCGGAGAAAACATCCATCTGGTTATTAAGCAGAATACGGCTGTCTCCCGGGACAATGCCCGGACCGAAGGTGgtgt                                                                                                                                                                     <  2:185370/233‑1 (MQ=255)                                                                            tGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCa                                                                                                                                                                  <  1:399413/250‑1 (MQ=255)                                                                            tGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTCCGGCCAGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTc                                                                                                                                                                   >  1:392264/1‑249 (MQ=255)                                                                                  ccTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCGCCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTGCCGAATGTGGTGTGCCGCGtata                                                                                                                                                           <  2:304182/250‑1 (MQ=255)                                                                                       aGTGAGGAGCCGGCTACCGCCTGGGGTACCTGTAACCAGCCAGGTTTTGCCGCCTTTCACCACAATGGTCGGCGACATCGACGACAGCGTGCGGTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCTGTTTGACGGAGAAATCATCCATCTGTTTATGAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGTTGTTCAGCGTCTAGGTc                                                                                                                                                       <  2:401806/250‑1 (MQ=255)                                                                                          gatgatCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACAGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGCGGTGTTCAGCGTATAGGTcaccg                                                                                                                                                   >  2:56016/1‑251 (MQ=255)                                                                                                             gggCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCt                                                                                                                                  <  1:323019/249‑1 (MQ=255)                                                                                                                cTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTACGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGATACCATCTTTATc                                                                                                                             >  1:32944/1‑251 (MQ=255)                                                                                                                    cGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATccacc                                                                                                                         <  1:347678/251‑1 (MQ=255)                                                                                                                            agccagGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTa                                                                                                                 <  1:173356/251‑1 (MQ=255)                                                                                                                             gccagGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACAGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTa                                                                                                                 >  1:334586/1‑250 (MQ=255)                                                                                                                              ccagGTTTTACCGTCTTTCACCCCAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAAt                                                                                                               <  1:31825/251‑1 (MQ=255)                                                                                                                               caggtaagaccgaatttgaccacAATGGTCGGCGACAGCGACGAAAGGGGGAGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCGCGTAAACGTTCGGTACGCGCGGATTGACGGAGAAATCAACCAGCTGGTTATTACGCAGAATACCGCTCTCGCCCGGGACAAAGCCAGTAACGAAAGCGGTGTTCAGCGTATAGGTGACCGCCACCGCGTTACCGTCTTTGTCAACCACTGAGTAATg                                                                                                              >  2:22102/19‑251 (MQ=255)                                                                                                                               cagGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGAAATg                                                                                                              >  1:259132/1‑251 (MQ=255)                                                                                                                                     ttACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGCGTAATGGGTAGt                                                                                                        >  2:386910/1‑251 (MQ=255)                                                                                                                                         cGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGa                                                                                                    <  2:140838/251‑1 (MQ=255)                                                                                                                                          gTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGAt                                                                                                   >  1:146110/1‑251 (MQ=255)                                                                                                                                                caccacAATGGTCGGCGACATCGCCGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTActct                                                                                             <  2:241466/251‑1 (MQ=255)                                                                                                                                                caccacAATGGTCGGCGAAATCGACGACAGCGGGCGTTTGTGCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTActct                                                                                             <  2:440575/251‑1 (MQ=255)                                                                                                                                                 accacAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACtctc                                                                                            <  1:129379/251‑1 (MQ=255)                                                                                                                                                      aaTGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAgg                                                                                                                                                         <  1:104804/185‑1 (MQ=255)                                                                                                                                                      aaTGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAgg                                                                                                                                                         >  2:104804/1‑185 (MQ=255)                                                                                                                                                           tCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTa                                                                                                                                                                                  >  1:11903/1‑155 (MQ=255)                                                                                                                                                             ggCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCCTAAGGCGCa                                                                                  >  1:86223/1‑249 (MQ=255)                                                                                                                                                                  cATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCAACACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGc                                                                           >  1:268276/1‑251 (MQ=255)                                                                                                                                                                             gCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAAtt                                                                >  1:285203/1‑251 (MQ=255)                                                                                                                                                                                   atttggacGGCACGACGGCGTTGGCATCACCGCCCACCAGCAAGTAAACGTACGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGAGAAAATGCCCGTAACGACGGTGGTGTTCAGCGTTTAGGTCACCGCCACCGCGTTACAATCGTTATTCACCACTGCGTAATGGGTAGTTTGATTAATCTCATACTGAGCAAGGTTGCCGGGGCGACTTTCgttg                                                          >  2:261695/8‑248 (MQ=255)                                                                                                                                                                                     ttGTTCGGCCCGACGGCGTTGGCATCCCCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCCCCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCGTGCCGGGGCGAATGTCGCTgg                                                         <  2:29575/250‑1 (MQ=255)                                                                                                                                                                                       gTTCGGCCCGACGGCGTTGGCATCACCGCCCTCCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTTTCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGTATTTCTCTGGAt                                                       <  2:52836/250‑1 (MQ=255)                                                                                                                                                                                            gCCCGACGGAGTTGGCATCACCGCCCACAAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCAACACTGAGTAATGGGTAGTTTGATTAATCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTc                                                 >  2:239739/1‑251 (MQ=255)                                                                                                                                                                                               cGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCt                                              >  1:199441/1‑251 (MQ=255)                                                                                                                                                                                                   ggCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTAt                                          >  1:46568/1‑251 (MQ=255)                                                                                                                                                                                                     cGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCTGGGCGAATTTCGCTGGATGTCTTCGCTTTATTg                                        <  1:390024/251‑1 (MQ=255)                                                                                                                                                                                                     cGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTg                                        >  1:400103/1‑251 (MQ=255)                                                                                                                                                                                                        tGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGat                                      >  1:407864/1‑250 (MQ=255)                                                                                                                                                                                                         ggCATCACCGCCCACTAGCCCGTAAACGTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTAGCCCGGGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGACGGCTTTGCTTTATTGatat                                    <  2:38414/251‑1 (MQ=255)                                                                                                                                                                                                                       ccaGCCCGTAAACGTTCGGGACGCCCGGTTTGACGGAGAAATCATCCAGCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTGTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCa                      <  2:46568/251‑1 (MQ=255)                                                                                                                                                                                                                             cGTAAACGTTCGGCACGCCCGGTTTGACGGAGAAATCATACATCTGGTTCTTAAGCAGAATCCCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGAGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCt                                                                                                                                  >  2:415901/1‑137 (MQ=255)                                                                                                                                                                                                                             cGTAAACGTGCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGTTTATTAAGCAGAATACCGCTCTCTCCCGCGACAATGCCCGTACCTAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCt                                                                                                                                  <  1:415901/137‑1 (MQ=255)                                                                                                                                                                                                                                    gTTCGGTACGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGATTTGGCAt         <  1:342221/251‑1 (MQ=255)                                                                                                                                                                                                                                        ggTACGCCCGGTTTGACGTAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTTGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGATTTGGCTTAGGc     <  2:212330/251‑1 (MQ=255)                                                                                                                                                                                                                                           aCGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCAGCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGATGTGGCATAGGCttt  <  2:215648/251‑1 (MQ=255)                                                                                                                                                                                                                                           aCGCCCGGTTTGACGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTGCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGATTTGGCATAGGCtt   <  2:219559/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               CACTGATGGTGGAAACGCGGCGCATTGGTCGCTTCGGCGACGTTCAAGCCATAATCGATGCTATTCACCACCATTTGCAGCACTGTAGTGATGATCCGGCTACCGCCTGGGCTACCGGTAACCAGCCAGGTTTTACCGTCTTTCACCACAATGGTCGGCGACATCGACGACAGCGGGCGTTTGTTCGGCCCGACGGCGTTGGCATCACCGCCCACCAGCCCGTAAACGTTCGGTACGCCCGGTTTGGCGGAGAAATCATCCATCTGGTTATTAAGCAGAATACCGCTCTCGCCCGCGACAATGCCCGTACCGAAGGTGGTGTTCAGCGTATAGGTCACCGCCACCGCGTTACCATCTTTATCCACCACTGAGTAATGGGTAGTTTGATTACTCTCATAAGGCGCAAGCTTGCCGGGGCGAATTTCGCTGGATGGCTTCGCTTTATTGATATCAATTTGATCGGCAATAGATTTGGCATAGGCTTT  >  NC_000913/3585271‑3585755

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