Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,828,654	C→T	G361R (GGG→AGG)	astD ←	aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,828,654	0	C	T	95.9%	138.4 / ‑5.2	49	G361R (GGG→AGG)	astD	aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (26/21);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCC‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCA  >  NC_000913/1828409‑1828884                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         gAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAACCCCCGCCCGCGCCTCCAGCAACAGCTGATCGCACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTAGCCATTCGAATCGCTTCATCGAACGCCTCATAACGCCAGACGCACAATAACGGTCCGCACACCTCTTCATCTGGTACGCCCGCCACGCCTGTCATTTCAATGATCCT‑CGGCg                                                                                                                                                                                                                                    >  1:248606/1‑251 (MQ=255)  aaTGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATACCCGCCCGCGCCTCCAGCAACAGTTGAACGAACTTTTCCCGCCCGGGGGAAACCAGACCGCAAGAGGGTCCGAAGCGTGTGTTATTCGCCATTCAAATCGCTTCATAGAAAGGATCATCACGCCAGACGCGCAATAACGGTCCGACCACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCgc                                                                                                                                                                                                                                   >  1:563395/1‑250 (MQ=255)       cgcgGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGc                                                                                                                                                                                                                               >  1:192907/1‑250 (MQ=255)         cgGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAAc                                                                                                                                                                                                                            <  1:44424/251‑1 (MQ=255)          ggTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACg                                                                                                                                                                                                                           <  2:282240/251‑1 (MQ=255)           gTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCCACg                                                                                                                                                                                                                           >  1:67099/1‑250 (MQ=255)                   agcaCCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTg                                                                                                                                                                                                                  >  1:435429/1‑251 (MQ=255)                      cccGGTAAGCGGGTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCtt                                                                                                                                                                                                               <  1:496482/250‑1 (MQ=255)                                     gTTCCAGTTAACCACCCCCGCCCGCGCCTCCCGCCACAGTTGATCGAACGTGTCCCGCCCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCGTCAGCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACAACTCTTCATCTGGTACGCCAGCAACGCCTGTCATGTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCgctgcg                                                                                                                                                                                                <  2:417074/251‑1 (MQ=255)                                          aGTTAACAATCCCCGCCCGCGCCTCCAGCAAAAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGc                                                                                                                                                                                           >  2:577070/1‑251 (MQ=255)                                             tAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCACGCAgg                                                                                                                                                                                        >  2:175288/1‑251 (MQ=255)                                              aaCAATCCCCGCCCTCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAATCGCGGCGCAAGCAggg                                                                                                                                                                                       <  1:485472/251‑1 (MQ=255)                                                 aaTCCCCGCCCACGCCTCCAGCAACAGTTGATCGAACGTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTTGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGgtc                                                                                                                                                                                    <  2:284296/251‑1 (MQ=255)                                                   cccccGCCCGCGCCTCCAGAAACAGTTTATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGgtcgt                                                                                                                                                                                  <  2:296608/250‑1 (MQ=255)                                                             cgcCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCg                                                                                                                                                                        >  1:393129/1‑251 (MQ=255)                                                             cgcCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCAGCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCTT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCg                                                                                                                                                                        <  2:450768/251‑1 (MQ=255)                                                                              gATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCAt                                                                                                                                                       <  2:155003/251‑1 (MQ=255)                                                                               aTCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCACGCCGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATg                                                                                                                                                      >  2:134569/1‑251 (MQ=255)                                                                                     cTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTa                                                                                                                                                 >  2:234822/1‑250 (MQ=255)                                                                                        ttCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAacc                                                                                                                                              >  1:6417/1‑250 (MQ=255)                                                                                         tCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCGGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGCTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAacca                                                                                                                                             <  2:540545/250‑1 (MQ=255)                                                                                           ccGCTCGGGGGAAACCAGACCGCAAGAGATCCCGAAGCGAGTGTTATTCGCCACTCGAATCGCTCCATCGAAAGTATCATAACGCCAGACGCGCAATAACGTTCCGAACACCTCGGCGTCTGGTACGCCATCAACGCCTTTCAGTTCAATGACCGT‑CGGCGGCAGCAACGATGTCCCTACTTGTAATAAGCGCGGAGCAAGGAGGGGTCGGGCGCCCACCGCTTCCAGTTGCTGCCATGCAGttgccacca                                                                                                                                          <  2:31797/251‑9 (MQ=255)                                                                                           ccGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACct                                                                                                                                          <  2:243425/251‑1 (MQ=255)                                                                                                cGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGGTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCAGCGCTTCCAGTTGCTGCCATGCAGTAACCACctgctg                                                                                                                                      <  2:292691/250‑1 (MQ=255)                                                                                                         ccAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGt                                                                                                                            <  2:230675/251‑1 (MQ=255)                                                                                                                   agagagTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATc                                                                                                                  >  1:4857/1‑250 (MQ=255)                                                                                                                        gTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGcc                                                                                                              >  2:569765/1‑249 (MQ=255)                                                                                                                                       tATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATACT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGc                                                                                             >  1:518307/1‑251 (MQ=255)                                                                                                                                       tATTCGCCATTCGAATCGCTTCATCGAAAGTATAATAACGCCAGACGCGAAATAACGGTCCGAAACCCTCTCCATCTGGGACGCCAGAAACGCTTGTCATTTCAATGATCCT‑CGGCGTAAGAAAGGATGTCCCTGTTTGTAATAAGCGCGGCGCAAGAAGGGGTGGTCCGCCCACCGCTTCCAGTTGCTGCCATGCAGTAACCCCCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTg                                                                                              <  1:283013/250‑1 (MQ=255)                                                                                                                                         ttCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCgg                                                                                           >  1:492120/1‑251 (MQ=255)                                                                                                                                           cGCCATTCGAATCGCTTCAGCAAAAGTATCATAACGCCAGACGCGAAATAACGGTCCGAACACCTCTTCATCTGTTACGCCGGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGCCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAACCAGCCCGCCAATAAACGACTGCGGtt                                                                                         <  1:8474/251‑1 (MQ=255)                                                                                                                                               aTTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTc                                                                                     >  1:143194/1‑251 (MQ=255)                                                                                                                                                  cGAATCGTTTCAACGAAATTATCATACCGCCAGCCGCGAAATACCGTCCCGAACCCCTCTACATCTGGTCCCCCAGCAACGCCTTTCATTTCAATGATCCC‑CGGCGTCACCAACGATGCCCCTTCTTGAAATAAGCCGGGCGCAAGCGGGGGTCGTCCGCCCATCGCTCACATTTCCCGCAATGCAGAACCCACCTGCTTTGCGGC‑TGTTTAAAAAACCAGCCCGCCAATAAACGGCTGCGGTTCGTCATc                                                                                  <  2:539069/251‑1 (MQ=255)                                                                                                                                                      tCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATccc                                                                                >  2:179858/1‑249 (MQ=255)                                                                                                                                                                gAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGAGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCAAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCACGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCg                                                                    >  1:514674/1‑251 (MQ=255)                                                                                                                                                                 aaaGTATCAAAACGCCAGACGCGCAAAAACGGTCCGAACCCCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGACCCT‑CGGCGTCAGCAACGATGCCCCTGCTTGTAAAAAGCGCGGCGCAAGCGGGGGTCGCCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTACCCACCGGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGt                                                                   <  1:99619/251‑1 (MQ=255)                                                                                                                                                                   attATCAAAACGCCAGACGCGCAATAACGGTCCTACCCCCTCTTCATCGGTTACCCCAGCAACGCTTGTCTTTTCAAGGATCCT‑CGGTTTAAGCAACGATGCCCCTGCTTGAAACAAGCGCGGCGAAAGCGGGGGCCGCCCCCCCATCGCTTCCAGTTTCCGCCATGCATAAACCACCTGCTGGTCGGCCTGTC‑CAAAAAACAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTa                                                                 <  1:537482/249‑1 (MQ=255)                                                                                                                                                                     tATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAAt                                                               >  1:514118/1‑251 (MQ=255)                                                                                                                                                                      aTCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGGTCGTCATCCCCGTTGCCCGGCGTTCATc                                                              >  2:476927/1‑251 (MQ=255)                                                                                                                                                                         atacccacagccgcgCAATAAGGGTCCGAACACCTCTTCTTCTGGTACCCCAGCAACGCTTGTAATTAAAAGGATCCT‑CGGCGGCAACAACGCTGTCCCTGCTTGTAATAAGCGTGGCGCAAACAGGGGTCGTCCGCCCATCGCTTCCTGTTGCTTCCTTGCATGAACCACCAGCGGTGCGGCCTGTC‑CAGAAACCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgct                                                           <  1:56374/240‑1 (MQ=255)                                                                                                                                                                              gCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTCCGCCAGCAACGCCTGTCATTTCACTGCTCCT‑CGGCGTCAGCAACGCTGTCCCTGCTTGTAATCAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctggctg                                                      >  2:200329/1‑251 (MQ=255)                                                                                                                                                                                  gACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACgg                                                  >  2:54806/1‑251 (MQ=255)                                                                                                                                                                                            aaCGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACg                                        >  2:442074/1‑251 (MQ=255)                                                                                                                                                                                            aaCGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGGTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACg                                        >  1:173737/1‑251 (MQ=255)                                                                                                                                                                                                   cGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCa                                    >  2:81540/1‑248 (MQ=255)                                                                                                                                                                                                      acacCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAAc                              >  1:214830/1‑251 (MQ=255)                                                                                                                                                                                                            ctccATCGGGTCCGCCAGCAACGCCTGTCATTCCAAGGATCCT‑GGGCGCCAGCAACGATGCCCCTGTTTGTAATAGGCGCGGCGCAAGCGGGGGTCTTCCCCCCACCCCTCCCAGTTGCTGCCATGCAAAAACCCCCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCAt                          <  1:134569/246‑1 (MQ=255)                                                                                                                                                                                                                              cAACGCCTGTCATTTAAATGACCCT‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCACCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCTTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTc      <  1:577070/251‑1 (MQ=255)                                                                                                                                                                                                                                  cccTGTCTTTTCA‑TGATCCCTGGGCGTCAGCAACGATGCCCCTGGGTTAAAAAAGCGCGGCGCAAGCAGGGGTCGCCCCCCAATCGCTACCTGTTGCTACCACGCAGTACCCACCGGGTGGGGGCCCTGTA‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCa  <  2:435429/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         GAATGGCGCGGTACTGGCAGCACCGGTAAGCGGTTTGTTCCAGTTAACAATCCCCGCCCGCGCCTCCAGCAACAGTTGATCGAACTTTTCCCGCTCGGGGGAAACCAGACCGCAAGAGAGTCCGAAGCGAGTGTTATTCGCCATTCGAATCGCTTCATCGAAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCC‑CGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTT‑CAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTGCGCCCCGCTTTTCAGCA  >  NC_000913/1828409‑1828884

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