Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,589,465	A→C	S29A (TCA→GCA)	unknown ←	PTS IIA‑like nitrogen‑regulatory protein PtsN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,589,465	0	A	C	100.0%	182.3 / NA	63	S29A (TCA→GCA)	unknown	PTS IIA‑like nitrogen‑regulatory protein PtsN
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/36);  total (27/36)

CACATATTTTCCGGAATCAACACCGCAAAGATTAAATCCACTTCACGGTGATCCGCTGCCTCATAGTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGAGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTCACACACATCCAGTGTTTCGCGCACATTGTGGTCGGTAATT  >  1/1589218‑1589706                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  cacaTATTTTCCGGAATCAACACCGCAAAGATTAAATCCACTTCACGGTGATCCGCTGCCTCATAGTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGAt                                                                                                                                                                                                                                                <  1:77997/251‑1 (MQ=255)                 tCAACACCGCAAAGATTAAATCCACTTCACGGTGATCCGCTGCCTCATAGTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGc                                                                                                                                                                                                                                <  1:53687/251‑1 (MQ=255)                   aaCACCGCAAAGATTAAATCCACTTCACGGTGATCCGCTGCCTCATAGTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCtt                                                                                                                                                                                                                              <  1:440224/251‑1 (MQ=255)                                             ccGGTGATCCGCTGCCTCATAGTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGAc                                                                                                                                                                                                    <  1:741754/250‑1 (MQ=255)                                                gTGATCCGCTGCCTCATAGTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGc                                                                                                                                                                                                                                >  1:879907/1‑220 (MQ=255)                                                gTGATCCGCTGCCTCATAGTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGa                                                                                                                                                                                                 >  1:822486/1‑251 (MQ=255)                                                          gCCTCATAGTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCgg                                                                                                                                                                                       >  1:608594/1‑251 (MQ=255)                                                           ccTCATAGTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCgg                                                                                                                                                                                       <  1:234230/250‑1 (MQ=255)                                                                  gTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGAt                                                                                                                                                                                                                                                >  1:54763/1‑186 (MQ=255)                                                                  gTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTtaata                                                                                                                                                                               >  1:802026/1‑251 (MQ=255)                                                                  gTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTtaata                                                                                                                                                                               >  1:408872/1‑251 (MQ=255)                                                                    cAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGAt                                                                                                                                                                                                                                                <  1:395533/184‑1 (MQ=255)                                                                    cAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTaataat                                                                                                                                                                             <  1:601237/251‑1 (MQ=255)                                                                    cAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTaataat                                                                                                                                                                             >  1:630883/1‑251 (MQ=255)                                                                      aTAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGAt                                                                                                                                                                                                                                                >  1:14297/1‑182 (MQ=255)                                                                              gTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACttctgtctcttatacacactccgagc                                                                                                                                                                    >  1:546360/1‑226 (MQ=255)                                                                              gTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTtaat                                                                                                                                                                                >  1:208940/1‑238 (MQ=255)                                                                              gTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAAtt                                                                                                                                                                   >  1:463222/1‑251 (MQ=255)                                                                              gTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTACGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAAtt                                                                                                                                                                   >  1:450131/1‑251 (MQ=255)                                                                                           aatACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAg                                                                                                                                                      <  1:647710/251‑1 (MQ=255)                                                                                                 gCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGaa                                                                                                                                                >  1:165880/1‑251 (MQ=255)                                                                                                    aTAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAAt                                                                                                                                              <  1:240949/250‑1 (MQ=255)                                                                                                        gTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCgg                                                                                                                                           <  1:100567/249‑1 (MQ=255)                                                                                                                    ccGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACg                                                                                                                             <  1:612943/251‑1 (MQ=255)                                                                                                                               ccTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTc                                                                                                                                                                           <  1:409186/194‑1 (MQ=255)                                                                                                                                cTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTaa                                                                                                                                                                      >  1:55972/1‑198 (MQ=255)                                                                                                                                       cATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTATCTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATttgttg                                                                                                           <  1:200742/250‑1 (MQ=255)                                                                                                                                            ccACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCAt                                                                                                      <  1:810292/250‑1 (MQ=255)                                                                                                                                                      gccgAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCg                                                                                           >  1:872763/1‑251 (MQ=255)                                                                                                                                                       ccgAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCgg                                                                                          <  1:15367/251‑1 (MQ=255)                                                                                                                                                          agaCCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCgg                                                                                          >  1:352989/1‑248 (MQ=255)                                                                                                                                                          agaCCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGgtg                                                                                        <  1:405094/250‑1 (MQ=255)                                                                                                                                                                        cAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGGCAAAATCATGCCCTGACGGATACTTTCCGGTGTtaat                                                                                                                                                                                <  1:568114/148‑1 (MQ=255)                                                                                                                                                                             tttCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCAt                                                                                                      <  1:743457/217‑1 (MQ=255)                                                                                                                                                                                           aaacaCTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAat                                                      <  1:401549/251‑1 (MQ=255)                                                                                                                                                                                              cacTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGTCGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCgg                                                                                          >  1:557996/1‑212 (MQ=255)                                                                                                                                                                                              cacTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAAtata                                                   <  1:534650/251‑1 (MQ=255)                                                                                                                                                                                              cacTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAAtata                                                   <  1:609624/251‑1 (MQ=255)                                                                                                                                                                                                cTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAg                                                 <  1:236937/251‑1 (MQ=255)                                                                                                                                                                                                cTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAg                                                 <  1:828567/251‑1 (MQ=255)                                                                                                                                                                                                  cAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGc                                                                                                                                                                                                                         <  1:462981/81‑1 (MQ=255)                                                                                                                                                                                                  cAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAgcg                                               <  1:889970/251‑1 (MQ=255)                                                                                                                                                                                                    aaaCAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTtgttgt                                                                                                          >  1:668661/1‑190 (MQ=255)                                                                                                                                                                                                    aaaCAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCgg                                                                                  <  1:558033/214‑1 (MQ=255)                                                                                                                                                                                                    aaaCAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAgcgcg                                             >  1:631394/1‑251 (MQ=255)                                                                                                                                                                                                         gCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTc                                                                                                                                                                           >  1:517078/1‑120 (MQ=255)                                                                                                                                                                                                           aTTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTcacac                                      <  1:630595/251‑1 (MQ=255)                                                                                                                                                                                                           aTTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTcacac                                      <  1:699824/251‑1 (MQ=255)                                                                                                                                                                                                              gTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCgg                                                                                  >  1:232755/1‑204 (MQ=255)                                                                                                                                                                                                              gTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTcacacaca                                   >  1:105939/1‑251 (MQ=255)                                                                                                                                                                                                               ttCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTCACACACAt                                  <  1:664720/251‑1 (MQ=255)                                                                                                                                                                                                               ttCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTCACACACAt                                  <  1:266454/251‑1 (MQ=255)                                                                                                                                                                                                                      ttttCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTCACACACATCCAgtgt                           >  1:671705/1‑251 (MQ=255)                                                                                                                                                                                                                      ttttCCAGCTGCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTCACACACATCCAgtg                            <  1:686755/250‑1 (MQ=255)                                                                                                                                                                                                                                gCAATTTTTCAGCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTgttgtt                                                                                                         <  1:678167/163‑1 (MQ=255)                                                                                                                                                                                                                                           gCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATc                                                                          <  1:620770/183‑1 (MQ=255)                                                                                                                                                                                                                                           gCCACAATGTGTGCGATTAACTCAAACACACGCTTTTTGCCGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAgcg                                               <  1:449603/210‑1 (MQ=255)                                                                                                                                                                                                                                              acaATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTCACACACATCCAGTGTTTCGCGCACATTGTGGTCGGtaat   <  1:385637/251‑1 (MQ=255)                                                                                                                                                                                                                                               caATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTtaataa                                                                                                                                                                              >  1:239664/1‑79 (MQ=255)                                                                                                                                                                                                                                               caATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTcacac                                      <  1:247668/215‑1 (MQ=255)                                                                                                                                                                                                                                               caATGTGTGCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTCACACACATCCAGTGTTTCGCGCACATTGTGGTCGGtaatt  >  1:586248/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gtgtgCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTc                                                                                                                           >  1:724321/1‑126 (MQ=255)                                                                                                                                                                                                                                                   gtgtgCGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAATCGCGCTCACACACATCCAGTGTTTCGCGCACATTGTGGt          >  1:150244/1‑239 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  CACATATTTTCCGGAATCAACACCGCAAAGATTAAATCCACTTCACGGTGATCCGCTGCCTCATAGTCAATAGGCGAGTCTAATTGAATGAATACGGCGATAGGTTGGGCTTCTGCCGGCAAACGTCCTTTTGGCATTGCCACGCCGCCGCCGAGACCGGAGTTGCCCAGTTTTTCACGGCTGAACAAACACTCAAAACAGCATTGTTCACCGTTTTCCAGCTGCAATTTTTCAGCCACAATGTGTGAGATTAACTCAAACACACGCTTTTTGCTGGACAAAATCATGCCCTGACGGATACTTTCCGGTGTTAATAATTCAGTAAATTTCATAAGCTAAAGTTGGAAATCGGCGCCTAAATAAACGTCTTTTACCAATTTGTTGTTCATAATCTCTTGCGGTGTGCCGGTGGCAATCATTTTACCGGCGCTCACAATATAAGCGCGCTCACACACATCCAGTGTTTCGCGCACATTGTGGTCGGTAATT  >  1/1589218‑1589706

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