Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,139,511	A→C	S4R (AGT→CGT)	unknown →	Molybdenum ABC transporter, substrate‑binding protein ModA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,139,511	0	A	C	100.0%	206.4 / NA	71	S4R (AGT→CGT)	unknown	Molybdenum ABC transporter, substrate‑binding protein ModA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (34/37);  total (34/37)

TGATGTGTTCCCGCTCGCCAAGCCGACCACAAACCTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTAGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTGAATGGGGATGTGGATGCGGCAGTTGT  >  1/1139274‑1139737                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                   tGATGTGTTCCCGCTCGCCAAGCCGACCACAAGCCTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAg                                                                                                                                                                                                                       <  1:880227/251‑1 (MQ=255) tGATGTGTTCCCGCTCGCCAAGCCGACCACAAACCTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAg                                                                                                                                                                                                                       <  1:582433/251‑1 (MQ=255)       gtTCCCGCTCGCCAAGCCGACCACAAACCTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCgg                                                                                                                                                                                                                 >  1:141071/1‑251 (MQ=255)                cGCCAAGCCGACCACAAACCTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGAtt                                                                                                                                                                                                         >  1:74769/1‑250 (MQ=255)                cGCCAAGCCGACCACAAACCTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTc                                                                                                                                                                                                        <  1:270989/251‑1 (MQ=255)                cGCCAAGCCGACCACAAACCTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTc                                                                                                                                                                                                        <  1:83877/251‑1 (MQ=255)                    aaGCCGACCACAAACCTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCtttt                                                                                                                                                                                                    <  1:776908/251‑1 (MQ=255)                          aCCACAAACCTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGc                                                                                                                                                                                               >  1:394644/1‑250 (MQ=255)                                  ccTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGa                                                                                                                                                                                      >  1:1967/1‑251 (MQ=255)                                    ttGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGaaa                                                                                                                                                                                    >  1:83215/1‑251 (MQ=255)                                          cacaAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCaaa                                                                                                                                                                               <  1:519232/250‑1 (MQ=255)                                          cacaAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAAt                                                                                                                                                                              <  1:754888/251‑1 (MQ=255)                                              aGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGcc                                                                                                                                                                                              >  1:131241/1‑231 (MQ=255)                                                                      gCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCGCCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGt                                                                                                                                                                         <  1:245088/228‑1 (MQ=255)                                                                           aTGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATgg                                                                                                                                             >  1:657679/1‑251 (MQ=255)                                                                                aTTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACg                                                                                                                                                                          >  1:702310/1‑217 (MQ=255)                                                                                        tcaGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAgg                                                                                                                                <  1:524891/251‑1 (MQ=255)                                                                                        tcaGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAgg                                                                                                                                <  1:446909/251‑1 (MQ=255)                                                                                             cgcAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCgg                                                                                                                           <  1:784929/251‑1 (MQ=255)                                                                                             cgcAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCgg                                                                                                                           <  1:184304/251‑1 (MQ=255)                                                                                               cAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGgag                                                                                                                         <  1:663187/251‑1 (MQ=255)                                                                                                     ccGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATc                                                                                                                   <  1:678086/251‑1 (MQ=255)                                                                                                     ccGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATc                                                                                                                   <  1:722127/251‑1 (MQ=255)                                                                                                        gtgtTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCttt                                                                                                                >  1:72950/1‑251 (MQ=255)                                                                                                          gtTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGagaaga                                                                                                                     >  1:119724/1‑244 (MQ=255)                                                                                                          gtTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGa                                                                                                              >  1:135111/1‑251 (MQ=255)                                                                                                          gtTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGa                                                                                                              >  1:681760/1‑251 (MQ=255)                                                                                                          gtTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGa                                                                                                              >  1:224545/1‑251 (MQ=255)                                                                                                          gtTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGa                                                                                                              >  1:6222/1‑251 (MQ=255)                                                                                                           ttttATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGaa                                                                                                             <  1:650044/251‑1 (MQ=255)                                                                                                             ttATGGCGTGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAAtt                                                                                                           >  1:52112/1‑251 (MQ=255)                                                                                                             ttATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAgg                                                                                                                                                                    >  1:701047/1‑194 (MQ=255)                                                                                                             ttATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAAtt                                                                                                           >  1:490669/1‑251 (MQ=255)                                                                                                             ttATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAAtt                                                                                                           >  1:200751/1‑251 (MQ=255)                                                                                                                 ggcatggcGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCg                                                                                                       >  1:776335/1‑251 (MQ=255)                                                                                                                       gcgcAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCg                                                                                                                                                            <  1:808301/192‑1 (MQ=255)                                                                                                                       gcgcAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTAt                                                                                                 >  1:323039/1‑251 (MQ=255)                                                                                                                                      ctgctTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAAt                                                                                  >  1:818917/1‑251 (MQ=255)                                                                                                                                        gctTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGc                                                                                                                                                       <  1:540441/180‑1 (MQ=255)                                                                                                                                           ttCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACaa                                                                             >  1:628762/1‑251 (MQ=255)                                                                                                                                                 gTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAgtggg                                                                       >  1:364220/1‑250 (MQ=255)                                                                                                                                                  tATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAgtgg                                                                        <  1:240140/249‑1 (MQ=255)                                                                                                                                                        tttGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTgaga                                                                <  1:798408/251‑1 (MQ=255)                                                                                                                                                                   gCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGgaga                                                                                                                        <  1:667478/184‑1 (MQ=255)                                                                                                                                                                   gCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGt                                                     >  1:615620/1‑251 (MQ=255)                                                                                                                                                                   gCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGt                                                     >  1:208541/1‑251 (MQ=255)                                                                                                                                                                         gATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAAcagc                                                                                          >  1:149621/1‑208 (MQ=255)                                                                                                                                                                          aTCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAAc                                                <  1:335196/249‑1 (MQ=255)                                                                                                                                                                          aTCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAAc                                                <  1:693094/249‑1 (MQ=255)                                                                                                                                                                           tCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCg                                                                                                                                                                                        >  1:215883/1‑112 (MQ=255)                                                                                                                                                                              aGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAgtgg                                                                        >  1:357290/1‑221 (MQ=255)                                                                                                                                                                                 cATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATg                                        <  1:634634/250‑1 (MQ=255)                                                                                                                                                                                         gCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTatat                                  <  1:328781/248‑1 (MQ=255)                                                                                                                                                                                         gCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTatat                                  <  1:469099/248‑1 (MQ=255)                                                                                                                                                                                            aaaCTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTg                            >  1:844253/1‑251 (MQ=255)                                                                                                                                                                                                     ccACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAAc                                                                                             >  1:527918/1‑177 (MQ=255)                                                                                                                                                                                                     ccACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTGAATGGGGa                    <  1:842248/250‑1 (MQ=255)                                                                                                                                                                                                     ccACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTGAATGGGGa                    <  1:900124/250‑1 (MQ=255)                                                                                                                                                                                                      caccgcaccgGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTGAATGGGGAtg                  <  1:380623/251‑1 (MQ=255)                                                                                                                                                                                                           caccgGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCATAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTGAATGGGGATGTGGa              <  1:94712/250‑1 (MQ=255)                                                                                                                                                                                                           caccgGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTGAATgg                       <  1:450875/241‑1 (MQ=255)                                                                                                                                                                                                             ccgGTGGTGTTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTAtattt                                <  1:858972/230‑1 (MQ=255)                                                                                                                                                                                                                   gtgtTGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAAAGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTGAATGGGGATGTGGa              >  1:862772/1‑242 (MQ=255)                                                                                                                                                                                                                       tGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGt                                                     <  1:636154/199‑1 (MQ=255)                                                                                                                                                                                                                       tGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGtt                                      <  1:637105/214‑1 (MQ=255)                                                                                                                                                                                                                       tGCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTGAATGGGGATGTGGATGCGGCAGTTGt  >  1:525591/1‑250 (MQ=255)                                                                                                                                                                                                                        gCATATTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTGAATGGGGATGTGGATGCGGCAGTTg   >  1:145311/1‑248 (MQ=255)                                                                                                                                                                                                                          atatTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAg                                              <  1:845280/203‑1 (MQ=255)                                                                                                                                                                                                                            atTCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATc                                                                                                                   <  1:290752/132‑1 (MQ=255)                                                                                                                                                                                                                             ttCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAAcagc                                                                                          <  1:117750/156‑1 (MQ=255)                                                                                                                                                                                                                             ttCCGGTGATGGCAATTCGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTAtattt                                >  1:605896/1‑214 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                   TGATGTGTTCCCGCTCGCCAAGCCGACCACAAACCTTGCTGCACAAGTAGAAAAAATTTTTCCGCAACTGCCGAATGTGATTTTTCATCAGGCGCAAGTGCCGGTGTTTTATGGCATGGCGCAAATGGTCACGCTGCTTTCCGAGTATGAATTTGATGAACAGCCGCAGATCCAGGCATGGCAGGCGAAACTAAAGCCACCGCACCGGTGGTGTTGCATATTCCGGTGATGGCAATTAGTCGGGATAAAAGCGCCGGCATTGATTCTTTTAAAGCCTTAGCGGAAAGCAAATTACGTTTAGGCGTTGGCGATAGCAAAGCGATGGCGTTGGGCAAAGGTGCGGAGAAGATCTTTGAATTGTCGGGCTATCAACAGCAGTTGAATGACAAAGTGGTGGTGAGAGCCGCGACGGTAAAACAGTTGATGTTATATTTATTGAATGGGGATGTGGATGCGGCAGTTGT  >  1/1139274‑1139737

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