Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,207,527	A→C	R102S (AGA→AGC)	unknown →	DNA‑binding domain of ModE / Molybdate‑binding domain of ModE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,207,527	0	A	C	100.0%	183.8 / NA	61	R102S (AGA→AGC)	unknown	DNA‑binding domain of ModE / Molybdate‑binding domain of ModE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (23/38);  total (23/38)

ATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGAATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGAAGTTATGTTGATGTTTAAAGC  >  1/1207296‑1207772                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      attaAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAgg                                                                                                                                                                                                                                    <  1:300599/251‑1 (MQ=255)         aaaTTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGc                                                                                                                                                                                                                              <  1:595960/249‑1 (MQ=255)               aaaaaTGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAg                                                                                                                                                                                                                       <  1:514039/250‑1 (MQ=255)                   aTGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGCTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGc                                                                                                                                                                                                                  <  1:41078/251‑1 (MQ=255)                           ccATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCg                                                                                                                                                                                                            >  1:404475/1‑249 (MQ=255)                           ccATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCg                                                                                                                                                                                                            >  1:235396/1‑249 (MQ=255)                           ccATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCCGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAAAAGTTTGCTTTCCg                                                                                                                                                                                                            >  1:866007/1‑249 (MQ=255)                             aTTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTAc                                                                                                                                                                                                        <  1:358072/251‑1 (MQ=255)                             aTTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTAc                                                                                                                                                                                                        <  1:326955/251‑1 (MQ=255)                             aTTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTAc                                                                                                                                                                                                        <  1:130746/251‑1 (MQ=255)                                                   ccAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGt                                                                                                                                                                                  <  1:32883/251‑1 (MQ=255)                                                         attttattcatggcttcaaggATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGcgca                                                                                                                                                                            >  1:736929/20‑251 (MQ=255)                                                             ttACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGcgcac                                                                                                                                                                           <  1:323705/248‑1 (MQ=255)                                                                      cgcCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGtt                                                                                                                                                                 <  1:886024/249‑1 (MQ=255)                                                                                  tcTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGATTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGtt                                                                                                                                                    >  1:448713/2‑250 (MQ=255)                                                                                  ccTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGtt                                                                                                                                                    >  1:70170/1‑250 (MQ=255)                                                                                  ccTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGtt                                                                                                                                                    >  1:810137/1‑250 (MQ=255)                                                                                  ccTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGtt                                                                                                                                                    >  1:134462/1‑250 (MQ=255)                                                                                  ccTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGATTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGtt                                                                                                                                                    >  1:423960/1‑250 (MQ=255)                                                                                       aaGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGc                                                                                                                                               <  1:649101/250‑1 (MQ=255)                                                                                       aaGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAG‑TTTTCGGCATCGTTCAGGCa                                                                                                                                              <  1:816881/250‑1 (MQ=255)                                                                                              gAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGt                                                                                                                                                                  <  1:412189/224‑1 (MQ=255)                                                                                              gAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGtt                                                                                                                                       >  1:107636/1‑251 (MQ=255)                                                                                               aaTAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGttt                                                                                                                                      <  1:660267/251‑1 (MQ=255)                                                                                               aaTAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGttt                                                                                                                                      >  1:128110/1‑251 (MQ=255)                                                                                                   aaaTCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACg                                                                                                                                         <  1:118215/244‑1 (MQ=255)                                                                                                    aaTCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAg                                                                                                                                 <  1:447150/251‑1 (MQ=255)                                                                                                    aaTCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAg                                                                                                                                 <  1:551209/251‑1 (MQ=255)                                                                                                     aTCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGa                                                                                                                                <  1:466716/251‑1 (MQ=255)                                                                                                         gCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCa                                                                                                                                              >  1:568438/1‑233 (MQ=255)                                                                                                          cccGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTAc                                                                                                                           <  1:234232/251‑1 (MQ=255)                                                                                                          cccGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTAc                                                                                                                           <  1:70091/251‑1 (MQ=255)                                                                                                               aaCCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGt                                                                                                                      <  1:604487/251‑1 (MQ=255)                                                                                                               aaCCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGt                                                                                                                      <  1:230593/251‑1 (MQ=255)                                                                                                                       cTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCa                                                                                                                                              >  1:819667/1‑219 (MQ=255)                                                                                                                           aaCGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAG‑TTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATc                                                                                                          <  1:680977/250‑1 (MQ=255)                                                                                                                                                         gcACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAAc                                                                             <  1:377108/250‑1 (MQ=255)                                                                                                                                                          cACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCg                                                                           <  1:506599/251‑1 (MQ=255)                                                                                                                                                              gcgcTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCgtgt                                                                        >  1:322901/1‑250 (MQ=255)                                                                                                                                                                           aTGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGaaaaa                                                          <  1:333437/251‑1 (MQ=255)                                                                                                                                                                             gcgcAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACAt                                                                                                                         >  1:152765/1‑186 (MQ=255)                                                                                                                                                                                cAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAgtgc                                                     <  1:78601/251‑1 (MQ=255)                                                                                                                                                                                   cGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCgtgt                                                                        <  1:416990/229‑1 (MQ=255)                                                                                                                                                                                   cGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAgtgcggt                                                  >  1:647080/1‑251 (MQ=255)                                                                                                                                                                                       tgctgcAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAg                                                                                                                                                                                >  1:223475/1‑121 (MQ=255)                                                                                                                                                                                       tgctgcAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTgg                                                                                                               >  1:615080/1‑186 (MQ=255)                                                                                                                                                                                       tgctgcAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGaa                                                                                                 >  1:589422/1‑200 (MQ=255)                                                                                                                                                                                        gctgcAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGtt                                                                                                                                                                 >  1:876063/1‑135 (MQ=255)                                                                                                                                                                                                          aCTTGAACAAACACAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTAc                                                                                                                                                                                                        <  1:829599/78‑1 (MQ=255)                                                                                                                                                                                                               aacaaacaCAAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGaa                      >  1:162904/1‑251 (MQ=255)                                                                                                                                                                                                                     cacaAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGa                       <  1:16047/244‑1 (MQ=255)                                                                                                                                                                                                                     cacaAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGAAGTTATg                <  1:332136/251‑1 (MQ=255)                                                                                                                                                                                                                      acaAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGAAGTTATGt               <  1:697877/251‑1 (MQ=255)                                                                                                                                                                                                                       caAGAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAg                                                                                                                                 <  1:834495/136‑1 (MQ=255)                                                                                                                                                                                                                          gAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGAAGTTATGt               <  1:350348/247‑1 (MQ=255)                                                                                                                                                                                                                          gAAAAGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGAAGTTATGTTGAt           >  1:383899/1‑251 (MQ=255)                                                                                                                                                                                                                           aaaaGGCCTTTGGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTAAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAAc                                                                             <  1:96588/184‑1 (MQ=255)                                                                                                                                                                                                                           aaaaGGCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGAAGTTATGTTGATg          <  1:266022/251‑1 (MQ=255)                                                                                                                                                                                                                               ggCCTTTAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGAAGTTATGTTGATGTTTa      >  1:472092/1‑251 (MQ=255)                                                                                                                                                                                                                                   tttAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGAAGTTATGTTGATGTTTAAAGc  <  1:52216/251‑1 (MQ=255)                                                                                                                                                                                                                                   tttAGCATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGAAGTTATGTTGATGTTTAAAGc  <  1:669234/251‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      ATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGTTTATGCGCAACGGTTGCTGCAACTTTATGATTTACTTGAACAAACACAAGAAAAGGCCTTTAGAATTTTGAAAGACGAAAAGGTTCCGCTCAACAGTTTGCTTTCCGCTACCGCGCGGTTTTCGCTGCAAAGTAGCGCACGCAATCAGTTTTTCGGCATCGTTCAGGCATTACGTTTTGAAGATGTACATTGTTTTGTGGATATCGAGCTGGAAGGCTTTTCCGAGCCATTAACCGTGTCGATTACAGAAAAAAGTGCGGTGCGTTTAAAACTTGAATTAAATAAGGAAGTTATGTTGATGTTTAAAGC  >  1/1207296‑1207772

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