Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,337,135	T→C	M110V (ATG→GTG)	unknown ←	3‑deoxy‑manno‑octulosonate cytidylyltransferase (EC 2.7.7.38)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,337,135	0	T	C	100.0%	198.0 / NA	66	M110V (ATG→GTG)	unknown	3‑deoxy‑manno‑octulosonate cytidylyltransferase (EC 2.7.7.38)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/33);  total (33/33)

GCGTAAATGCCGATGTGACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCATGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAA  >  1/1336891‑1337371                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             gCGTAAATGCCGATGTGACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATaa                                                                                                                                                                                                                                         >  1:68524/1‑250 (MQ=255) gCGTAAATGCCGATGTGACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATaa                                                                                                                                                                                                                                         >  1:26891/1‑250 (MQ=255)  cGTAAATGCCGATGTGACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGAGCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAAc                                                                                                                                                                                                                                        <  1:455331/250‑1 (MQ=255)        tGCCGATGTGACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGa                                                                                                                                                                                                                                 >  1:564111/1‑251 (MQ=255)                tgACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCgg                                                                                                                                                                                                                         <  1:201377/251‑1 (MQ=255)                 gACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCg                                                                                                                                                                                                                          >  1:825075/1‑249 (MQ=255)                      aaaTATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTtcgtc                                                                                                                                                                                                                   <  1:542001/251‑1 (MQ=255)                        atatTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTtcgtcg                                                                                                                                                                                                                  >  1:314633/1‑250 (MQ=255)                        atatTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGc                                                                                                                                                                                                                 <  1:201993/251‑1 (MQ=255)                         tatTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGcc                                                                                                                                                                                                                <  1:458575/251‑1 (MQ=255)                            tCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTg                                                                                                                                                                                                             <  1:107544/251‑1 (MQ=255)                               gTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAAt                                                                                                                                                                                                          >  1:62773/1‑251 (MQ=255)                                 aTTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGt                                                                                                                                                                                                        <  1:61940/251‑1 (MQ=255)                                     aGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGa                                                                                                                                                                                                     >  1:72584/1‑250 (MQ=255)                                      gCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGt                                                                                                                                                                                                        >  1:84992/1‑246 (MQ=255)                                                  ccGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCAtc                                                                                                                                                                                        <  1:797493/250‑1 (MQ=255)                                                         gTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAAttt                                                                                                                                                                                            >  1:174263/1‑239 (MQ=255)                                                                      aaTTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGttttt                                                                                                                                                                    >  1:222238/1‑250 (MQ=255)                                                                      aaTTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGttttt                                                                                                                                                                    >  1:633870/1‑250 (MQ=255)                                                                       aTTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTcc                                                                                                                                                                  <  1:388352/251‑1 (MQ=255)                                                                         tGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTccac                                                                                                                                                                <  1:312701/251‑1 (MQ=255)                                                                            tCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTccac                                                                                                                                                                <  1:122648/248‑1 (MQ=255)                                                                                         gAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAAATTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGt                                                                                                                                                >  1:358297/1‑251 (MQ=255)                                                                                                 gcgcGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAg                                                                                                                                                                         >  1:502934/1‑218 (MQ=255)                                                                                                          aaaTAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATtgtg                                                                                                                               >  1:365693/1‑251 (MQ=255)                                                                                                           aaTAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTgg                                                                                                                              >  1:317196/1‑251 (MQ=255)                                                                                                            aTAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGt                                                                                                                             >  1:133523/1‑251 (MQ=255)                                                                                                            aTAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGt                                                                                                                             >  1:466936/1‑251 (MQ=255)                                                                                                                 gCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAAGTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGaa                                                                                                                                    <  1:130103/239‑1 (MQ=255)                                                                                                                         ccGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCa                                                                                                                 <  1:235030/250‑1 (MQ=255)                                                                                                                             ttttGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGcac                                                                                                            <  1:721049/251‑1 (MQ=255)                                                                                                                                          cACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGTGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGcac                                                                                                            <  1:879398/238‑1 (MQ=255)                                                                                                                                          cACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGccc                                                                                                  <  1:200924/248‑1 (MQ=255)                                                                                                                                           aCTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACCTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGaa                                                                                              >  1:106069/1‑251 (MQ=255)                                                                                                                                            cTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGaaa                                                                                             >  1:490006/1‑251 (MQ=255)                                                                                                                                              tttACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAAtt                                                                                           >  1:290586/1‑251 (MQ=255)                                                                                                                                                                     ttCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAAtt                                                                    <  1:538875/251‑1 (MQ=255)                                                                                                                                                                      tCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTTCCACCaaat                                                                    >  1:255042/1‑248 (MQ=255)                                                                                                                                                                      tCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGGCACCCAtt                                                                    >  1:55726/1‑250 (MQ=255)                                                                                                                                                                          tCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAATGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTa                                                               <  1:445441/251‑1 (MQ=255)                                                                                                                                                                          tCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTa                                                               <  1:350340/251‑1 (MQ=255)                                                                                                                                                                           cgcgTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCacac                                                                                                          <  1:468809/207‑1 (MQ=255)                                                                                                                                                                           cgcgTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTAt                                                              <  1:895635/251‑1 (MQ=255)                                                                                                                                                                                   caATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGGTTCCACCAATTCGTTTTCGGTGGCg                                                      >  1:351418/1‑251 (MQ=255)                                                                                                                                                                                      tttttACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCCTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAGTTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGccccccactccttttcgggggcggta                                                   >  1:443080/1‑227 (MQ=255)                                                                                                                                                                                            ccGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCCATTCGttttc                                                             >  1:242670/1‑232 (MQ=255)                                                                                                                                                                                            ccGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGAcg                                             <  1:139232/251‑1 (MQ=255)                                                                                                                                                                                            ccGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGAcg                                             <  1:149433/251‑1 (MQ=255)                                                                                                                                                                                                         cATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCgg                                <  1:199111/251‑1 (MQ=255)                                                                                                                                                                                                           tGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTc                                                                   <  1:827147/214‑1 (MQ=255)                                                                                                                                                                                                           tGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGGCGCGGTCGGCACCGGAt                              <  1:45459/251‑1 (MQ=255)                                                                                                                                                                                                           tGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGAt                              <  1:603107/251‑1 (MQ=255)                                                                                                                                                                                                           tGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATACACACTTCCGCCCCGAAATTTCTTGCGGCGGTTTCCACCACTTCGTTTTCGGGGGCGGTAATGACGCCGTCGGcccccggt                              >  1:448830/1‑244 (MQ=255)                                                                                                                                                                                                           tGTTCACAGGGTATTTTGCAAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGAt                              <  1:805891/251‑1 (MQ=255)                                                                                                                                                                                                            gTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTTTCGGTGGCGGTAATGACGCGGGCGGCACCGGAtt                             >  1:821645/1‑251 (MQ=255)                                                                                                                                                                                                            gTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAAATCCTTTTTGGTGGGGATAATGAAGCGGGCGGCACCGGAtt                             >  1:868895/1‑251 (MQ=255)                                                                                                                                                                                                                   gggTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCCATTCGTTATCGGGGGCGATAATGACGCGGGCGGCACCGGGTTGttttgg                      >  1:343817/1‑248 (MQ=255)                                                                                                                                                                                                                    ggTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGc                      <  1:823464/250‑1 (MQ=255)                                                                                                                                                                                                                     gTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCCCCAATTCGTTATCGGGGGGGGTAATTACGCGGGCGGCACCGGGTTGTTTgggtt                    >  1:806276/1‑248 (MQ=255)                                                                                                                                                                                                                     gTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCCATTCGTTATCGGGGGGGGTAATGACGCGGGCGGCACCGGGTTGttttggtt                    >  1:781112/1‑246 (MQ=255)                                                                                                                                                                                                                                 gTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCCCCAATTCGTTATCGGGGGGGGTAATGACGCGGTCGGCACCCGGTTGTTTTGCTTGTTTTAAAACgg        >  1:759217/1‑250 (MQ=255)                                                                                                                                                                                                                                 gTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCCATTCGTTATCGGTGGCGATAATGACGCGGGCGGCACCGGATTGTTTGGCTTGTTCaaaaaac         >  1:90981/1‑248 (MQ=255)                                                                                                                                                                                                                                 gTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCaaaaa           <  1:603961/248‑1 (MQ=255)                                                                                                                                                                                                                                 gTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAAc          <  1:415414/249‑1 (MQ=255)                                                                                                                                                                                                                                       cGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGaa  <  1:501495/251‑1 (MQ=255)                                                                                                                                                                                                                                           aCTTGTTGCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGtt     <  1:868176/244‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             GCGTAAATGCCGATGTGACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCATGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAA  >  1/1336891‑1337371

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