Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,495,194	A→T	F339Y (TTT→TAT)	ahr ←	NADPH‑dependent aldehyde reductase Ahr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,495,194	0	A	T	100.0%	149.6 / NA	46	F339Y (TTT→TAT)	ahr	NADPH‑dependent aldehyde reductase Ahr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (29/17);  total (29/17)

AAACTATTCTCGCGCGAATGCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAAAAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAAT  >  NC_000913/4494951‑4495440                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       aaaCTATTCTCGCGCGAATGCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGc                                                                                                                                                                                                                                                 >  1:293042/1‑251 (MQ=255)  aaCTATTCTCGCGCGAATGCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCt                                                                                                                                                                                                                                                <  1:51360/251‑1 (MQ=255)                            gggCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCtt                                                                                                                                                                                                                      >  1:187664/1‑251 (MQ=255)                              gCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGGGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATTGGCTTTCAACACCACGCGGTAACGCGCCTTAc                                                                                                                                                                                                                    <  1:53123/251‑1 (MQ=255)                               cTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTAcc                                                                                                                                                                                                                   >  1:37829/1‑251 (MQ=255)                               cTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATAGGCTTTCAACACCACGCGGTAACGCGCCTTAcc                                                                                                                                                                                                                   >  1:198645/1‑251 (MQ=255)                                  ggCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGAGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTc                                                                                                                                                                                                                >  1:479446/1‑251 (MQ=255)                                          aTGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGcaca                                                                                                                                                                                                         >  1:296081/1‑250 (MQ=255)                                               ggTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTgg                                                                                                                                                                                                   <  1:349547/251‑1 (MQ=255)                                                    cGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATggg                                                                                                                                                                                              >  1:333549/1‑250 (MQ=255)                                                        aCTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGcgtc                                                                                                                                                                                           <  2:259094/250‑1 (MQ=255)                                                          tAAAAAGTCAGTTTGAGGAACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAAAGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGtt                                                                                                                                                                                        >  1:79272/1‑251 (MQ=255)                                                            aaaaGTCAATTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCTCATTGATATCAACCATGATATTGCCAATGTCAACGAACACTGTCGGCAGGCTGAGCTGGCGATACGACCAAAAAGAATATGTTCGAAAGAGGGAAGCGCCACAGATCAGCGCTGCGAATGATTTTCACTAATCGGCTTTCAACACCACGCGGTAACGCGACTTACCGTCGCTCCCATGCTGGATGGcgtcgatat                                                                                                                                                                                      >  2:51585/1‑247 (MQ=255)                                                               aGTCATTTTGAGGCACTGGAGCGGCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAAttt                                                                                                                                                                                   <  2:342189/251‑1 (MQ=255)                                                                 tCAGTTTGAGGCACTGGAGCGTCCACAAGCTGATGAACAGGATATTGTGCGCATTGATATCAACCATGATAGTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACATAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAAGCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATTGCTTCGTTAAtttt                                                                                                                                                                                  <  2:90094/250‑1 (MQ=255)                                                                   aGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACAACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCg                                                                                                                                                                                >  1:186274/1‑250 (MQ=255)                                                                        gAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCg                                                                                                                                                                          <  1:26849/251‑1 (MQ=255)                                                                                ggAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGc                                                                                                                                                                                              >  2:338455/1‑223 (MQ=255)                                                                                ggAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGc                                                                                                                                                                                              <  1:338455/223‑1 (MQ=255)                                                                                          caAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCg                                                                                                                                                        <  1:383578/251‑1 (MQ=255)                                                                                                 aTGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAAc                                                                                                                                                 <  2:102572/251‑1 (MQ=255)                                                                                                       aGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTtgct                                                                                                                                           >  1:89576/1‑251 (MQ=255)                                                                                                              tgtgCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTgcgggcg                                                                                                                                    >  1:216628/1‑251 (MQ=255)                                                                                                                     aTTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAAc                                                                                                                             <  1:169272/251‑1 (MQ=255)                                                                                                                          tatCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGAGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATc                                                                                                                        >  2:7784/1‑251 (MQ=255)                                                                                                                           atCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATcc                                                                                                                       >  2:296879/1‑250 (MQ=255)                                                                                                                                 cATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGCCATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTAc                                                                                                                 >  1:281702/1‑251 (MQ=255)                                                                                                                                     atatTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGc                                                                                                               >  1:349398/1‑249 (MQ=255)                                                                                                                                      tatTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGAGCGAAAGAGGGCCGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCCCATGCTGGATGGCGACGTTAATTTTCGACATCGGGAACCGTTCGGTGGTCGGGGCAACCTTGCTGCGGGCGGCAAAACGCATCAGcctcagcagc                                                                                                            >  2:211401/1‑242 (MQ=255)                                                                                                                                          gCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCa                                                                                                         <  1:40473/250‑1 (MQ=255)                                                                                                                                                gTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGc                                                                                                   <  1:296879/250‑1 (MQ=255)                                                                                                                                                                                        aaaaCAGAATATGTGCGAAAGCGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAAt                                                          >  2:438253/1‑251 (MQ=255)                                                                                                                                                                                        aaaaCAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATTATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCTCCTTACCGGCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAAt                                                          <  2:79052/251‑1 (MQ=255)                                                                                                                                                                                        aaaaCAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCACAAGGCGTGCCGGTAGCAGAACCAGCGACGCTGCGATCGCCCGCAAt                                                          >  2:159419/1‑251 (MQ=255)                                                                                                                                                                                               aaTATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACAAGAGACGCTGCGCTCGCCCGCAATTAACga                                                    >  2:425613/1‑249 (MQ=255)                                                                                                                                                                                                           agagGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAac                                        >  2:195484/1‑250 (MQ=255)                                                                                                                                                                                                               ggCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAacagaca                                   <  2:37829/251‑1 (MQ=255)                                                                                                                                                                                                                    cgcCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTGCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAgcggc                              <  2:296921/251‑1 (MQ=255)                                                                                                                                                                                                                     gcCTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGTTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAgcggcg                             >  1:480083/1‑251 (MQ=255)                                                                                                                                                                                                                       cTCAGATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTg                           >  2:231670/1‑251 (MQ=255)                                                                                                                                                                                                                           gATCAGCGCTGCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTgaga                        >  1:16363/1‑250 (MQ=255)                                                                                                                                                                                                                              cAGCGCTGCGAATGATTTTCAATAATAGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCTCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCAGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGCACCAGAGACCCTGCGATCGCCCGCAATTAACGTAAATGCCGGAACAGCCAGCGGCGTGAGCAccc                    >  2:400120/1‑250 (MQ=255)                                                                                                                                                                                                                                     gCGAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCaccgac             >  1:345778/1‑251 (MQ=255)                                                                                                                                                                                                                                       gAATGATTTTCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCaccgaccg           >  1:21545/1‑251 (MQ=255)                                                                                                                                                                                                                                              tttCAATAAGCGGGTTGCAACACCACGCGGTAACGCGCCTTACCGGCGCGCACATGCTGGATGGCGTCGTGAATTTTCGACATCGGGAACAGTGCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCGGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGa     <  2:281702/250‑1 (MQ=255)                                                                                                                                                                                                                                                tCAATAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAAt  >  1:370345/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       AAACTATTCTCGCGCGAATGCAGCGTCGGGCTGGGCATTTTATGCCGGTAGCGTTACTAAAAAGTCAGTTTGAGGCACTGGAGCGTCCACAAGCAGATGAACAGGATATTGTGCGCATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAAAAATCGGCTTTCAACACCACGCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAAT  >  NC_000913/4494951‑4495440

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