Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,817,609	T→C	I6M (ATA→ATG)	unknown ←	Ribonucleotide reductase of class Ia (aerobic), alpha subunit (EC 1.17.4.1)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,817,609	0	T	C	100.0%	165.7 / NA	48	I6M (ATA→ATG)	unknown	Ribonucleotide reductase of class Ia (aerobic), alpha subunit (EC 1.17.4.1)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/29);  total (19/29)

ATGGAACACCGCAAGGCGCGCGGCAAGGTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACTATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAATGTGTTGGATGTGAAAATTCAAGTACTAAATTTTTTACTTGAATTATTGACATAAAATAACTGCTT  >  1/1817363‑1817767                                                                                                                                                                                                                                                       |                                                                                                                                                               aTGGAACACCGCAAGGCGCGCGGCAAGGTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCa                                                                                                                                                            >  1:37881/1‑251 (MQ=255)        acCGCAAGGCGCGCGGCAAGGTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCtt                                                                                                                                                     >  1:542510/1‑251 (MQ=255)         ccGCAAGGCGCGCGGCAAGGTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTg                                                                                                                                                    >  1:95282/1‑251 (MQ=255)            cAAGGCGCGCGGCAAGGTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAACGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTc                                                                                                                                                 <  1:431535/251‑1 (MQ=255)                 cgcgcgGCAAGGTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGa                                                                                                                                            <  1:56618/251‑1 (MQ=255)                 cgcgcgGCAAGGTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGa                                                                                                                                            <  1:176653/251‑1 (MQ=255)                 cgcgcgGCAAGGTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGa                                                                                                                                            <  1:31939/251‑1 (MQ=255)                        cAAGGTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATAcc                                                                                                                                     <  1:116108/251‑1 (MQ=255)                           ggTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCt                                                                                                                                                      <  1:299200/231‑1 (MQ=255)                            gTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTc                                                                                                                                                 <  1:543497/235‑1 (MQ=255)                                      tCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAat                                                                                                                       <  1:899487/251‑1 (MQ=255)                                          gggAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATc                                                                                                                                                             <  1:550829/209‑1 (MQ=255)                                              aTCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAatat                                                                                                                     >  1:236820/1‑245 (MQ=255)                                                       gATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGa                                                                                                                                            <  1:703532/213‑1 (MQ=255)                                                            gATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAg                                                                                                 >  1:686145/1‑251 (MQ=255)                                                             aTCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGt                                                                                                <  1:484829/251‑1 (MQ=255)                                                               ccgccgccgCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTg                                                                                              <  1:393666/251‑1 (MQ=255)                                                               ccgccgccgCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTg                                                                                              <  1:411629/251‑1 (MQ=255)                                                               ccgccgccgCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTg                                                                                              <  1:635138/251‑1 (MQ=255)                                                               ccgccgccgCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTg                                                                                              <  1:45850/251‑1 (MQ=255)                                                                        cTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATAcc                                                                                                                                     <  1:240381/203‑1 (MQ=255)                                                                          ttAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGa                                                                                                                                            <  1:65128/194‑1 (MQ=255)                                                                           taataaTGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGa                                                                                                                                            >  1:901431/1‑193 (MQ=255)                                                                            aataatGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAAtt                                                                                 >  1:244419/1‑251 (MQ=255)                                                                                             tGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTc                                                                                                                                                 >  1:116920/1‑170 (MQ=255)                                                                                               tCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGa                                                                                                                                         <  1:415712/176‑1 (MQ=255)                                                                                                   aGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAATGTGTTGGa                                                           <  1:114815/250‑1 (MQ=255)                                                                                                       gCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAAtg                                                                  >  1:736504/1‑239 (MQ=255)                                                                                                          aaTCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAATGTGTTGGATGTGaaaa                                                   >  1:271618/1‑251 (MQ=255)                                                                                                               cTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTtc                                                                                                                               <  1:172211/170‑1 (MQ=255)                                                                                                                 tCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGAt                                                                                                                                        <  1:745053/159‑1 (MQ=255)                                                                                                                  cATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATAcc                                                                                                                                     <  1:449643/161‑1 (MQ=255)                                                                                                                   aTAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGcc                                                                                                                                                       <  1:76987/142‑1 (MQ=255)                                                                                                                   aTAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGa                                                                                                                                            >  1:306996/1‑153 (MQ=255)                                                                                                                       aaCTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGa                                                                                                                                            >  1:286773/1‑149 (MQ=255)                                                                                                                                tgAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAATGTGTTGGATGTGAAAATTCAAGTACTAAATTTTTTAc                              <  1:726585/250‑1 (MQ=255)                                                                                                                                              cACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATAcc                                                                                                                                     <  1:550919/133‑1 (MQ=255)                                                                                                                                               aCTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTAttt                                                                                                                                >  1:750488/1‑137 (MQ=255)                                                                                                                                                     gaCACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGaaaa                                                                                                              >  1:690530/1‑149 (MQ=255)                                                                                                                                                         cGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAATGTGTTGGATGTGAAAATTCa                                               >  1:73011/1‑208 (MQ=255)                                                                                                                                                         cGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAATGTGTTGGATGTGAAAATTCAAGTACTAAATTTTTTACTTGAATTATTg                   >  1:495549/1‑236 (MQ=255)                                                                                                                                                           gACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAATGTGTTGGATGTGAAAATTCAAGTACTAAATTTTTTACTTGAATTATTGACATAAAATAACTGCtt  >  1:160668/1‑251 (MQ=255)                                                                                                                                                               caTTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGt                                                                                                <  1:179656/153‑1 (MQ=255)                                                                                                                                                               caTTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAAtgt                                                                 >  1:562892/1‑184 (MQ=255)                                                                                                                                                                   gTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACatat                                                                                                   >  1:322895/1‑146 (MQ=255)                                                                                                                                                                                        aGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGa                                                                                                                                            <  1:473130/84‑1 (MQ=255)                                                                                                                                                                                        aGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAATGTGTTGGATGAGAAAATTCAAGTACTAAATTTTTTACTTGAATTATTGAc                 <  1:850862/207‑1 (MQ=255)                                                                                                                                                                                         gTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACCATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAAtg                                                                  <  1:492386/157‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                               ATGGAACACCGCAAGGCGCGCGGCAAGGTATTGATAATCCGGGGAATCTTTGCTGATCAGATCCGCCGCCGCTTTAATAATGGTTTCGTGAATGTCGGAGGTGCGAATCCCTTCATAAAACTGAATGTGAGAACGCAATTCCACTTGAGACACGGACACATTGTCCAACCCTTCCGCCGCCCAAGTGATCACGCGGTGAATTTTATCAAGATTAATCGGTTCCGGTTTTCCGTCACGTTTGGTGACTATCAATGCCTTGTTCATGGAAGATACCTATTTCTCTAAAATATATGAAAAACACAACATATAGTTGCGTTATAAAAATTTTGCACAAGATAATGTGTTGGATGTGAAAATTCAAGTACTAAATTTTTTACTTGAATTATTGACATAAAATAACTGCTT  >  1/1817363‑1817767

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