Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,978,375	T→C	P172P (CCA→CCG)	unknown ←	RTX toxin transporter, ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,978,375	0	T	C	100.0%	176.1 / NA	60	P172P (CCA→CCG)	unknown	RTX toxin transporter, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (38/22);  total (38/22)

CGAGCCACCGTATCCCCAACCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAATGACGAAAAAGCCTTGCTCCCAATTCGACATCAATACGACTTGTACTATGAGCAAAAACATAGGTTCTCAAGCCTCCCAAGATTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATTACATTCAATGTTGAAAATCCACGGTGAACCAGGACTTTATCCATCACGACTTGAAAAAAGAGTGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACGGCAGGAATAAACCAGGTAAAGTCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGATGCAATTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGAATTAGATCAAAAATGAGGTAACGATCGGTT  >  1/1978130‑1978615                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 cGAGCCACCGTATCCCCAACCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGgtgt                                                                                                                                                                                                                                             <  1:836711/251‑1 (MQ=255)     ccACCGTATCCCCAACCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAAt                                                                                                                                                                                                                                          >  1:56748/1‑250 (MQ=255)     ccACCGTATCCCCAACCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATc                                                                                                                                                                                                                                         <  1:490094/251‑1 (MQ=255)     ccACCGTATCCCCAACCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATc                                                                                                                                                                                                                                         <  1:67996/251‑1 (MQ=255)     ccACCGTATCCCCAACCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATc                                                                                                                                                                                                                                         <  1:555804/251‑1 (MQ=255)        ccGTATCCCCAACCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGt                                                                                                                                                                                                                                      >  1:609030/1‑251 (MQ=255)            aTCCCCAACCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATc                                                                                                                                                                                                                                         <  1:373170/244‑1 (MQ=255)               cccAACCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCaaaaa                                                                                                                                                                                                                               <  1:391527/251‑1 (MQ=255)                  aaCCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATc                                                                                                                                                                                                                                         <  1:375627/238‑1 (MQ=255)                                                    gCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGaaa                                                                                                                                                                                                             <  1:301673/232‑1 (MQ=255)                                                    gCTAATAAGTGACGAAAAAGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATaaa                                                                                                                                                                                          >  1:440328/1‑251 (MQ=255)                                                                  aaaaaGTCTTGCTCCCAGTTCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTAtt                                                                                                                                                                            >  1:180737/1‑251 (MQ=255)                                                                                 cagttCGACATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAg                                                                                                                                                                                                                <  1:115039/197‑1 (MQ=255)                                                                                         cATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATc                                                                                                                                                                                                                                         >  1:356743/1‑167 (MQ=255)                                                                                         cATCAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGt                                                                                                                                                                                                   >  1:491926/1‑205 (MQ=255)                                                                                            cAATACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACGACAATAGCAAGAGTCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTa                                                                                                                                                  <  1:455203/251‑1 (MQ=255)                                                                                             aaTACGACTTGTGCTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATc                                                                                                                                                                                                                                         >  1:108699/1‑163 (MQ=255)                                                                                                       gtgcTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATc                                                                                                                                                                                                                                         >  1:849596/4‑153 (MQ=255)                                                                                                        tgcTATGAGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAAttt                                                                                                                                       <  1:378813/248‑1 (MQ=255)                                                                                                               aGCAAACACATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAAcc                                                                                                                                                       <  1:179879/227‑1 (MQ=255)                                                                                                                    acaCATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTc                                                                                                                            >  1:95941/3‑249 (MQ=255)                                                                                                                     caCATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTa                                                                                                                         >  1:83531/2‑251 (MQ=255)                                                                                                                       cATAGGTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAAttt                                                                                                                                       >  1:703163/1‑235 (MQ=255)                                                                                                                           ggTTCTCAGCCCCCCCAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGAGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAg                                                                                                                   >  1:401112/1‑251 (MQ=255)                                                                                                                            gttctcagcccccccAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGaaa                                                                                                                                                                                                             <  1:758222/148‑1 (MQ=255)                                                                                                                            gttctcagcccccccAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAAt                                                                                                                     >  1:300907/13‑248 (MQ=255)                                                                                                                            gttctcagcccccccAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGa                                                                                                                  >  1:527488/13‑251 (MQ=255)                                                                                                                                       ccccAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCGATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGacg                                                                                                       >  1:62356/2‑249 (MQ=255)                                                                                                                                       ccccAAGACTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGacg                                                                                                       >  1:63219/2‑249 (MQ=255)                                                                                                                                                           aGACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCtt                                                                                   >  1:132698/1‑251 (MQ=255)                                                                                                                                                            gACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCttt                                                                                  >  1:227416/1‑251 (MQ=255)                                                                                                                                                            gACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCttt                                                                                  >  1:560011/1‑251 (MQ=255)                                                                                                                                                            gACTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGCGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTAc                                                                                                                        >  1:764145/1‑213 (MQ=255)                                                                                                                                                              cTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGaaa                                                                                                                                                                                                             >  1:284510/1‑126 (MQ=255)                                                                                                                                                              cTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATc                                                                                            >  1:829824/1‑239 (MQ=255)                                                                                                                                                              cTACAATAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGa                                                                                >  1:349983/1‑251 (MQ=255)                                                                                                                                                                    tAGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAAt                                                                                                                     >  1:876964/1‑208 (MQ=255)                                                                                                                                                                     aGCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAAt                                                                                                                                                                                                                              >  1:198157/1‑102 (MQ=255)                                                                                                                                                                      gCAAGAGCCACTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGtt                                                                         >  1:98813/1‑250 (MQ=255)                                                                                                                                                                              cACTGTAATCACGTTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTactc                                                                                                                                                                                                                                         <  1:785753/82‑4 (MQ=255)                                                                                                                                                                                cTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAAcc                                                                                                                                                       <  1:483790/162‑1 (MQ=255)                                                                                                                                                                                cTGTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCa                                                              >  1:702999/1‑251 (MQ=255)                                                                                                                                                                                aggTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCa                                                              >  1:47721/3‑251 (MQ=255)                                                                                                                                                                                  gTAATCACATTCAATGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATc                                                                                                                                                                                                                                         <  1:484897/78‑1 (MQ=255)                                                                                                                                                                                              aaTGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAAc                                                <  1:724083/251‑1 (MQ=255)                                                                                                                                                                                               aTGTTGAAAATCCACGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGc                                                                                                                                                                                  >  1:811098/1‑120 (MQ=255)                                                                                                                                                                                                             cGGTGAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAAc                                                >  1:598554/1‑236 (MQ=255)                                                                                                                                                                                                                 gAACCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCAtt                             >  1:333185/1‑251 (MQ=255)                                                                                                                                                                                                                  aaCCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATg                                                                                                                <  1:141377/167‑1 (MQ=255)                                                                                                                                                                                                                  aaCCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGt                                                                          >  1:683780/1‑205 (MQ=255)                                                                                                                                                                                                                   aCCAAGACCTTATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAg                           >  1:683547/1‑251 (MQ=255)                                                                                                                                                                                                                     caagaccttATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGAt                         <  1:314564/244‑1 (MQ=255)                                                                                                                                                                                                                        gaccttATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGaa                                                                                                                                                             <  1:446681/112‑1 (MQ=255)                                                                                                                                                                                                                        gaccttATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGcgcg                                                                                                            >  1:713415/5‑165 (MQ=255)                                                                                                                                                                                                                         accttATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTAc                                                                                                                        >  1:862197/4‑152 (MQ=255)                                                                                                                                                                                                                           cttATCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCaaaaa                   <  1:497234/250‑1 (MQ=255)                                                                                                                                                                                                                              aTCCATCACGACTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGaaa                                                                                                                                                                                                             >  1:100517/1‑62 (MQ=255)                                                                                                                                                                                                                                        aCTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGAt       <  1:72530/250‑1 (MQ=255)                                                                                                                                                                                                                                         cTTGAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCg     >  1:414826/1‑251 (MQ=255)                                                                                                                                                                                                                                            gAAAAAAGAGCGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACAGCAGGAATAAACCAGGTAAAATCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGtt  >  1:83480/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 CGAGCCACCGTATCCCCAACCCGACGAGCTTCAAAATAAGAAATCGGTAAAGCTAATAAATGACGAAAAAGCCTTGCTCCCAATTCGACATCAATACGACTTGTACTATGAGCAAAAACATAGGTTCTCAAGCCTCCCAAGATTATTTCAAATAAGACTACAATAGCAAGAGCCACTGTAATTACATTCAATGTTGAAAATCCACGGTGAACCAGGACTTTATCCATCACGACTTGAAAAAAGAGTGGTGTAATCAGTGCAAAAATTTGCAAAAATATAGAAACGATTAACGTTTCTATAAAAATTTTGCGGTATTTAATTACGGCAGGAATAAACCAGGTAAAGTCGAATTTTGCTAAATTTCCTACAATAGATGCGCGTGATGCAATTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGAATTAGATCAAAAATGAGGTAACGATCGGTT  >  1/1978130‑1978615

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