Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	42,219	G→A	P91L (CCA→CTA)	unknown ←	Cytochrome c heme lyase subunit CcmL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	42,219	0	G	A	100.0%	218.0 / NA	75	P91L (CCA→CTA)	unknown	Cytochrome c heme lyase subunit CcmL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/39);  total (36/39)

AAAAACAAATTAACGCGGCAATCAGCGTTAACACAAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATGGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTGTTATTCTGACATTGCGGACAGCGCAGTTCCCGCGTTAAGGCGTGATAGTCGGCTTCCTGCTGCGGCGAGCTGAAATTCAGCACGTCGATGGCGGCAAACACATT  >  1/41981‑42463                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     aaaaaCAAATTAACGCGGCAATCAGCGTTAACACAAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGGTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAg                                                                                                                                                                                                                                          <  1:85029/251‑1 (MQ=255)             aCGCGGCAATCAGCGTTAACACAAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGc                                                                                                                                                                                                                              <  1:644667/251‑1 (MQ=255)                  gCAATCAGCGTTAACACAAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAg                                                                                                                                                                                                                         >  1:231190/1‑251 (MQ=255)                     aTCAGCGTTAACACAAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGa                                                                                                                                                                                                                      >  1:312500/1‑251 (MQ=255)                          cGTTAACACAAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAACATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAA‑CCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATccac                                                                                                                                                                                                                 <  1:382813/250‑1 (MQ=255)                               acacAAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGt                                                                                                                                                                                                            >  1:398780/1‑251 (MQ=255)                               acacAAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGt                                                                                                                                                                                                            <  1:878692/251‑1 (MQ=255)                               acacAAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGt                                                                                                                                                                                                            <  1:230718/251‑1 (MQ=255)                                 acaAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTAt                                                                                                                                                                                                          <  1:162394/251‑1 (MQ=255)                                 acaAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTAt                                                                                                                                                                                                          <  1:208292/251‑1 (MQ=255)                                  caAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTAtt                                                                                                                                                                                                         >  1:711770/1‑251 (MQ=255)                                           ttGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACAt                                                                                                                                                                                                <  1:63320/251‑1 (MQ=255)                                               ccACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTg                                                                                                                                                                                            <  1:621588/251‑1 (MQ=255)                                                cacaAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATc                                                                                                                                                                                               >  1:571082/1‑247 (MQ=255)                                                       tCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTAt                                                                                                                                                                                                          <  1:860196/229‑1 (MQ=255)                                                        cATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAg                                                                                                                                                                                                                                          <  1:661403/196‑1 (MQ=255)                                                                 ttttCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGcac                                                                                                                                                                                                     >  1:532644/1‑224 (MQ=255)                                                                    tCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCAT‑AAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAAc                                                                                                                                                                      <  1:357824/251‑1 (MQ=255)                                                                     ccTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCACCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAAc                                                                                                                                                                      <  1:451850/251‑1 (MQ=255)                                                                              aCAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCaaa                                                                                                                                                              >  1:63781/1‑250 (MQ=255)                                                                              aCAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGCGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCaaa                                                                                                                                                              >  1:534443/1‑250 (MQ=255)                                                                              aCAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTCCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAac                                                                                                                                                             <  1:455835/251‑1 (MQ=255)                                                                                aaaGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAacac                                                                                                                                                           <  1:680247/251‑1 (MQ=255)                                                                                   gCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGAtt                                                                                                                                                                                    <  1:809684/223‑1 (MQ=255)                                                                                   gCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTg                                                                                                                                                                          >  1:132971/1‑233 (MQ=255)                                                                                   gCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACttt                                                                                                                                                        >  1:82863/1‑251 (MQ=255)                                                                                   gCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACttt                                                                                                                                                        >  1:145033/1‑251 (MQ=255)                                                                                      ttAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGt                                                                                                                                                     >  1:99989/1‑251 (MQ=255)                                                                                               tgctgTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATc                                                                                                                                            <  1:263827/251‑1 (MQ=255)                                                                                                 ctgTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGAtt                                                                                                                                                                                    <  1:625824/209‑1 (MQ=255)                                                                                                 ctgTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGCGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCca                                                                                                                                          >  1:265956/1‑251 (MQ=255)                                                                                                   gttgtGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCcaca                                                                                                                                        >  1:36253/1‑251 (MQ=255)                                                                                                     tgtgCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGc                                                                                                                                      <  1:688845/251‑1 (MQ=255)                                                                                                     tgtgCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGc                                                                                                                                      <  1:632275/251‑1 (MQ=255)                                                                                                      gtgCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGt                                                                                                                                                                                                            <  1:363646/180‑1 (MQ=255)                                                                                                      gtgCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTAt                                                                                                                                                                                                          >  1:862878/1‑182 (MQ=255)                                                                                                       tgCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCaa                                                                                                                                    <  1:754135/251‑1 (MQ=255)                                                                                                          ttCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAg                                                                                                                                 <  1:691186/251‑1 (MQ=255)                                                                                                             aCTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTg                                                                                                                                                                          <  1:257842/207‑1 (MQ=255)                                                                                                             aCTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAgcgg                                                                                                                              >  1:274532/1‑251 (MQ=255)                                                                                                              cTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCTTGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAgcgg                                                                                                                              >  1:567640/1‑250 (MQ=255)                                                                                                              cTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACaa                                                                                                                                                                    <  1:561520/212‑1 (MQ=255)                                                                                                              cTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAg                                                                                                                                 <  1:394218/247‑1 (MQ=255)                                                                                                              cTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTCTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAgcgg                                                                                                                              >  1:205978/1‑250 (MQ=255)                                                                                                                   ggCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAg                                                                                                                        >  1:489797/1‑251 (MQ=255)                                                                                                                     ccTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGaa                                                                                                                      <  1:882119/251‑1 (MQ=255)                                                                                                                     ccTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGaa                                                                                                                      <  1:199056/251‑1 (MQ=255)                                                                                                                           gcgcGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAg                                                                                                                                                                                                                                          >  1:160271/1‑129 (MQ=255)                                                                                                                            cgcgGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCg                                                                                                                <  1:596407/250‑1 (MQ=255)                                                                                                                                gCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTAt                                                                                                                                                                                                          <  1:204301/156‑1 (MQ=255)                                                                                                                                gCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATAtt                                                                                                           >  1:220409/1‑251 (MQ=255)                                                                                                                                 cATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTg                                                                                                          <  1:283877/251‑1 (MQ=255)                                                                                                                                 cATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTg                                                                                                          <  1:263685/251‑1 (MQ=255)                                                                                                                                            ccACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGcac                                                                                                                                                                                                     <  1:204778/149‑1 (MQ=255)                                                                                                                                            ccACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTg                                                                                                                                                    >  1:601915/1‑198 (MQ=255)                                                                                                                                              aCTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGc                                                                                                                                                                                                       >  1:276659/1‑145 (MQ=255)                                                                                                                                                    cATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTGTTATTCTGACATTGCGGAc                                                                                       <  1:552964/251‑1 (MQ=255)                                                                                                                                                                          aGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAAc                                                                                                                                                                      >  1:828342/1‑150 (MQ=255)                                                                                                                                                                                 ccAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTGTTATTCTGACATTGCGGACAGCGCAGTTCCCGCGTTAAGGCGTGATAg                                                          <  1:873309/251‑1 (MQ=255)                                                                                                                                                                                         cATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGcac                                                                                                                                                                                                     >  1:673072/1‑104 (MQ=255)                                                                                                                                                                                         cATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCcac                                                                                                                                         <  1:214533/164‑1 (MQ=255)                                                                                                                                                                                            aaGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTgc                                                                                                                                                   >  1:80661/1‑151 (MQ=255)                                                                                                                                                                                              gATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGcac                                                                                                                                                                                                     <  1:139098/99‑1 (MQ=255)                                                                                                                                                                                                    aGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTAtt                                                                                                                                                                                                         >  1:342661/1‑89 (MQ=255)                                                                                                                                                                                                             cGTAACCCATAAAACAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTGTTATTCTGACATTGCGGACAGCGCAGTTCCCGCGTTAAGGCGTGATAGTCGGCTTCCTGCTGCGGCGAGCTGaaa                               <  1:70504/250‑1 (MQ=255)                                                                                                                                                                                                                        aaaCAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAac                                                                                                                                                             >  1:767226/1‑113 (MQ=255)                                                                                                                                                                                                                        aaaCAAGGGTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTGTTATTCTGACATTGCGGACAGCGCAGTTCCCGCGTTAAGGCGTGATAGTCGGCTTCCTGCTGCGGCGAGCTGAAATTCAGCACGTCg                   >  1:376585/1‑251 (MQ=255)                                                                                                                                                                                                                                gTGGCAACGGTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTGTTATTCTGACATTGCGGACAGCGCAGTTCCCGCGTTAAGGCGTGATAGTCGGCTTCCTGCTGCGGCGAGCTGAAATTCAGCACGTCGATggcgg            <  1:222547/250‑1 (MQ=255)                                                                                                                                                                                                                                        ggTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATcc                                                                                                                                           >  1:524076/1‑115 (MQ=255)                                                                                                                                                                                                                                        ggTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTGTTATTc                                                                                                    >  1:596064/1‑154 (MQ=255)                                                                                                                                                                                                                                        ggTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTGTTATTCTGACATTGCGGACAGCGCAGTTCCCGCGTTAAGGCGTGATAGTCGGCTTCCTGCTGCGGCGAGCTGAAATTCAGCACGTCGATGGCGGCAAACACAt   >  1:369819/1‑251 (MQ=255)                                                                                                                                                                                                                                         gTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATccacca                                                                                                                                                                                                               >  1:133910/1‑46 (MQ=255)                                                                                                                                                                                                                                         gTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAacac                                                                                                                                                           >  1:189627/1‑98 (MQ=255)                                                                                                                                                                                                                                         gTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCcac                                                                                                                                         >  1:53488/1‑116 (MQ=255)                                                                                                                                                                                                                                         gTCAATAGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTGTTATTCTGACATTGCGGACAGCGCAGTTCCCGCGTTAAGGCGTGATAGTCGGCTTCCTGCTGCGGCGAGCTGAAATTCAGCACGTCGATGGCGGCAAACACAtt  >  1:825497/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     AAAAACAAATTAACGCGGCAATCAGCGTTAACACAAAAATGCTTGTCCACAAATTCATTAGTCATTTTCCTTATTTAACAAAGCATTAAGGCGCTGCTGTTGTGCTTCACTTAAGGCCTGTGGCGCGGCATTATCCGCCCCACTCGGCATCGTTCTTTTCCGCCCGAAAAGCCATACCAGCCCGCATAAGATCAAAGCAAAAGGCGTAACCCATAAAACAAGGGTGGCAACGGTCAATGGCGGGTCATAAGTAACAAAATTGCCATAGCGATCCACCATGTATTGCACCACATCCTGTTTGGATTTGCCTTCCTGTAACAATTCAAACACTTTAGTGCGCATATCCACAGCAATAGCGGCGTTAGAATCGGCGATATTGTTATTCTGACATTGCGGACAGCGCAGTTCCCGCGTTAAGGCGTGATAGTCGGCTTCCTGCTGCGGCGAGCTGAAATTCAGCACGTCGATGGCGGCAAACACATT  >  1/41981‑42463

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