Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,647,428	C→T	A377V (GCG→GTG)	gor →	glutathione reductase (NADPH)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,647,428	0	C	T	91.5%	129.0 / NA	47	A377V (GCG→GTG)	gor	glutathione reductase (NADPH)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/24);  minor base G (0/4);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGA  >  NC_000913/3647208‑3647668                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 aaCACGGGGGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGcc                                                                                                                                                                                                                     <  1:174159/250‑1 (MQ=255)          gCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCg                                                                                                                                                                                                           >  1:165813/1‑251 (MQ=255)                        acCGGTGGCAGTTGCAGGGGGTCGCCGTCGCTCTGAACGCCTGTTAAATAACAAGCTGGATGAGCATCTGGATTACAGCAACATTCCGCCCGTGGTCTCCCGCCATCCGCCGATTGGTCCTGTTGGTTTAACGGAACCGCAGTCGCGCGAGCAGTATGGCGACGCTCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGTTGTATACCGCCGTCACCACTCCCCGCCAGCCGTGCCGCATGAAGCTGgtgt                                                                                                                                                                                             <  1:236675/251‑1 (MQ=255)                             tGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGATGGTTTAACGGAACCGCAGTCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGtt                                                                                                                                                                                        >  1:51545/1‑251 (MQ=255)                              ggCAGTTGCAGCGGGTCGTCGTCTCTCTTAACGCCTGTTTAATAACAAGCCGGATGAGCAACTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCTTCCGCCGATGGTTCCTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTAGGGCGACGATCAGGTGAAAGTGTATAAATCCCCTTTCACCGTGAGGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGtt                                                                                                                                                                                        <  1:320570/250‑1 (MQ=255)                                  gTTGAAGCGGGGCGCCGGCTCTCTGAACGCCTGTTTAAAAACAATCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTAGGGCGACGATCAGGTGAAAGGGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATc                                                                                                                                                                                   <  2:258699/251‑1 (MQ=255)                                              cGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCACCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagagaaga                                                                                                                                                                        <  1:211884/250‑1 (MQ=255)                                                ccGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTg                                                                                                                                                                     >  1:269200/1‑251 (MQ=255)                                                   tctctcTGAACGCCTGTTTAATAACAAGCCGGATGAGCCTCTGGATTACAGCAACATTCCGAACGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGGGTGCGTTGGATCGgaag                                                                                                                                                                              >  2:34001/1‑239 (MQ=255)                                                   tctctcTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaag                                                                                                                                                                              <  1:34001/239‑1 (MQ=255)                                                   actctctGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGCTTACAGCAACATTCCGACAGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGTACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCATCTTTCACCGTGATGTGTACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCg                                                                                                                                                                  >  2:14154/2‑251 (MQ=255)                                                    ctctctGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCgg                                                                                                                                                                 >  1:281679/1‑251 (MQ=255)                                                        ctGAACGCCGGTTTAAAAACAAGCCGGATGAGCATTGGGATTACAGCAACAGTCCGCCCGTGGTCTCCAGCCACCCGCCGATGGGTACCGTTGGGTTAATGGAACCGCAGCCGCGCGAGCAGTATGGCGACGATCAGGAGAAAGTGTATAAATCCTCTGTCACGGTGATGTAGACCGCCGTCACCACTCACCGCCAGCCGCGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTa                                                                                                                                                               <  1:143842/249‑1 (MQ=255)                                                              gCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTCTGGCGCCGATCAGGTGAAAGTGCATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGAGCGTTGGATCGGAAGAGCAGATTGTCGGTCTTCCCGGc                                                                                                                                                       >  2:122294/1‑251 (MQ=255)                                                                 tGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCCCCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGTAGATTGTCGGTATTCACGGCAtt                                                                                                                                                    >  1:302103/1‑251 (MQ=255)                                                                     taataaCAAGCCGGATGAGCATCTGGATTACAGCAACTGTCAGACCGTGGCCTTCAGCCAGCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGACGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATCCCGCCGTCACCACGCACCGCCAGCCGTGCCTCATTAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCTCGTCATTGGc                                                                                                                                                 <  2:286770/250‑1 (MQ=255)                                                                     aaataaAAAGCCGGATGAGCATTTGGATTACAGAACCTTTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGGTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAATGGTATAAATCCTCTGTCACCGTGATTTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCt                                                                                                                                                <  2:10362/250‑1 (MQ=255)                                                                           cAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagagaa                                                                                                                                                                          >  2:22855/1‑219 (MQ=255)                                                                           cAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagagaa                                                                                                                                                                          <  1:22855/219‑1 (MQ=255)                                                                            aaGCCGGATGAGCATCGGGATTACAGAAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTGGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTAt                                                                                                                                         <  2:136386/251‑1 (MQ=255)                                                                             aGCCGGATGAGCATCGGGATTACCGCAACATTCCGACCGTGTTCTTCAGCCACCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTGCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATg                                                                                                                                        <  2:252752/251‑1 (MQ=255)                                                                                          tCTGGATTACAGCAACATTCCGCCCGTGGTCTTCAGCCATCCGCCGATTGGAACTGTTGGTTTACCGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGc                                                                                                                           >  2:65629/1‑251 (MQ=255)                                                                                               aTTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGc                                                                                                                      >  1:167192/1‑251 (MQ=255)                                                                                                       aaCATTCCGACCGCGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGAATGGCGACGATCAGGTGCAAGCGTATAAATCCTCTCTCACCGTGATGCATACCGCCGTCAACACTCACCGCCAGCCGt                                                                                                                                                                                                               >  1:167634/1‑154 (MQ=255)                                                                                                            tCCGACTGTGGTCTGCTGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGGGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCGCGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGgcg                                                                                                          <  2:21006/250‑1 (MQ=255)                                                                                                                 ccGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGAGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGAAGGGCTTCGCGGTGGCGCTGa                                                                                                      >  2:317469/1‑249 (MQ=255)                                                                                                                   gTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGa                                                                                                      >  1:247104/1‑247 (MQ=255)                                                                                                                      gTCTTCAGCCACCCGCCGATTGGTACTGTTGGTTTACCGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATTCCGCCGTCACCACCCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATggg                                                                                               <  2:274470/251‑1 (MQ=255)                                                                                                                        cTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATggggg                                                                                             >  2:4131/1‑251 (MQ=255)                                                                                                                         ttCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGc                                                                                            >  2:1306/1‑251 (MQ=255)                                                                                                                          tCAGCCACCCGCCGATTGGTACGGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCCCCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATggggg                                                                                             <  2:247104/249‑1 (MQ=255)                                                                                                                                        aTTGGTACTGTTGGTTTAACGGAACCGCATGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAACCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTGTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTTAAGATGGGGGCAACCAAAAAAGACtt                                                                             <  2:309054/251‑1 (MQ=255)                                                                                                                                                 gTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATAc                                                                    <  2:122946/251‑1 (MQ=255)                                                                                                                                                        tAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTAGcc                                                             >  1:128690/1‑251 (MQ=255)                                                                                                                                                                    aGGCGCGCGAGCGGTATGGCGACAATCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTGCGACAATACCGTCGCCATTCACCCAAcg                                                 <  2:206381/251‑1 (MQ=255)                                                                                                                                                                           cgAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGGGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCag                                          <  1:241567/251‑1 (MQ=255)                                                                                                                                                                            gAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGGGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCaga                                         <  2:105245/251‑1 (MQ=255)                                                                                                                                                                                         acgaTCAGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCCCCCAACGGCGGCAGAAGAGTTCGTGAc                             >  2:138695/1‑250 (MQ=255)                                                                                                                                                                                           gaTCAGGTGAAAGTGTATAAATCCCCTTTCACCGTGATGTATACCGCCGTCACCACACACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGTGCAACCAAAAAAGACTTCGACAATACTGTCGCCATTCGCCCAACGGCGGCAGAAGAGTTCGTGACAAt                          <  1:253945/251‑1 (MQ=255)                                                                                                                                                                                               aGGTGAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGt                      <  2:167192/251‑1 (MQ=255)                                                                                                                                                                                               aGGTGAAAGTGTATAAATCCTCTTTCACCGGGAGGTATACCGCCGGCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGt                      <  2:67713/251‑1 (MQ=255)                                                                                                                                                                                                   gAAAGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTaa                   >  1:221954/1‑250 (MQ=255)                                                                                                                                                                                                     aaGTGTATAAATCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATg                <  1:309540/251‑1 (MQ=255)                                                                                                                                                                                                      aGTGTATAAATCCTCTTGCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGt               <  1:246807/251‑1 (MQ=255)                                                                                                                                                                                                        tgtATAAATCCTCTTTCACCGTGATGTATACCGCCGTCCCCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGtt              <  1:273817/250‑1 (MQ=255)                                                                                                                                                                                                             aaaTCCTCTTTCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAgg         <  1:307197/250‑1 (MQ=255)                                                                                                                                                                                                                      ttCACCGTGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAaga  <  1:179212/248‑1 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 AACACGGGTGCAGTGGAGCTGACACCGGTGGCAGTTGCAGCGGGTCGCCGTCTCTCTGAACGCCTGTTTAATAACAAGCCGGATGAGCATCTGGATTACAGCAACATTCCGACCGTGGTCTTCAGCCATCCGCCGATTGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCACCCAACGGCGGCAGAAGAGTTCGTGACAATGCGTTAAATGTTAAAGGGCTAAGA  >  NC_000913/3647208‑3647668

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