Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,573,001	G→A	L168L (CTG→TTG)	glgB ←	1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,573,001	0	G	A	100.0%	150.8 / NA	43	L168L (CTG→TTG)	glgB	1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/20);  total (23/20)

GAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAGGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACC  >  NC_000913/3572758‑3573235                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           gAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCAt                                                                                                                                                                                                                                      <  1:292528/250‑1 (MQ=255)  aaCTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTGTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCg                                                                                                                                                                                                                                    <  2:294930/251‑1 (MQ=255)     ttCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGt                                                                                                                                                                                                                                 >  1:400241/1‑251 (MQ=255)                  aTTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCAGCGGGTGACGGCGACCGTc                                                                                                                                                                                                                    <  1:333674/251‑1 (MQ=255)                                              tGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCga                                                                                                                                                                                         <  2:75450/250‑1 (MQ=255)                                                 gcTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGAcc                                                                                                                                                                                     >  1:215773/1‑251 (MQ=255)                                                         gTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGc                                                                                                                                                                             <  1:35049/251‑1 (MQ=255)                                                             gTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCTCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGttt                                                                                                                                                                         <  2:256417/251‑1 (MQ=255)                                                                   ccTTGTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGcc                                                                                                                                                                   <  2:274343/251‑1 (MQ=255)                                                                       ttCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGAAGTTGAATTGCCCAACCACCGAGCCCCGACGGGCGGTTGAAGccc                                                                                                                                                                  >  2:182220/1‑248 (MQ=255)                                                                             catcccGCAAAGAAGAGAGGCGGTTGCCGGGCGCAGTGGCGCTTCAAAGGCATAATGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCAGCTCGTATTTATAGAGCTTACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGTCGCTCTCTTTATGCAAGCTCATCGGTTGACGTCGACCGTTCCAGTAGTTGAATTTCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCTGACAgaga                                                                                                                                                         <  2:416990/248‑1 (MQ=255)                                                                                       aTAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTAcc                                                                                                                                               >  1:235738/1‑251 (MQ=255)                                                                                          agagCCGCTGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCtgt                                                                                                                                            <  2:36515/251‑1 (MQ=255)                                                                                              aCGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACg                                                                                                                                        <  1:218017/251‑1 (MQ=255)                                                                                               cgcgGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGc                                                                                                                                       >  1:156891/1‑251 (MQ=255)                                                                                                   ttttCTGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGCCGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTTTTTACGCAAGCGCATCGGGTGACGGCGCCCGTCCCAGTAGTTGCATTGCCCACCCCCCGAGACCCGACCGGCGTTTGGAGCCCAGGCAACGAAACGCGTACCttaattgccctc                                                                                                                                   >  2:28450/2‑240 (MQ=255)                                                                                                       ccGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCCCCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATa                                                                                                                               <  2:235738/251‑1 (MQ=255)                                                                                                       ccGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATa                                                                                                                               <  1:40139/251‑1 (MQ=255)                                                                                                              cATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCAGTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCGACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTg                                                                                                                        <  2:342790/251‑1 (MQ=255)                                                                                                                          aaaGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGCCAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGAGCGCCt                                                                                                             >  2:292590/1‑250 (MQ=255)                                                                                                                                        cGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACgg                                                                                               <  1:144374/250‑1 (MQ=255)                                                                                                                                         ggACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGgc                                                                                             <  1:365685/251‑1 (MQ=255)                                                                                                                                                 gACGCAAGTTGCCATTGGCATCAATCATCTCGTCTTTATCGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGAATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACAACCGAGACCCGACGGGCGTTTGGAGCCCAGAAAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGAATGCGCGCCTAAGGTTTCATACGGGCGCAGGTg                                                                                      >  1:353177/1‑250 (MQ=255)                                                                                                                                                   cGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTACGGTTTCATACGGGCGCAGGTGAGt                                                                                   >  2:208529/1‑251 (MQ=255)                                                                                                                                                    gCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGt                                                                                   >  2:160570/1‑250 (MQ=255)                                                                                                                                                                      caatcaTCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGcc                                                                  >  1:357179/1‑249 (MQ=255)                                                                                                                                                                        atcatcTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGc                                                              >  1:295149/1‑251 (MQ=255)                                                                                                                                                                         tcatcTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGc                                                              >  1:13032/1‑250 (MQ=255)                                                                                                                                                                                     tataGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCCTGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTgg                                                 >  2:341777/1‑251 (MQ=255)                                                                                                                                                                                       taGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGCCGGGCGTTTGGAGCCCAGAAAGAGAAACGCGTACATGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATCCGGGCGCAGGTGAGGACCTTCAGATAATAGCCAGGCATCCATTTcatngct                                               >  2:383261/1‑245 (MQ=255)                                                                                                                                                                                        agagCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGAtt                                              <  1:81499/251‑1 (MQ=255)                                                                                                                                                                                         gagCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGAtt                                              >  1:74205/1‑250 (MQ=255)                                                                                                                                                                                            cTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCg                                          >  1:82828/1‑251 (MQ=255)                                                                                                                                                                                               aCCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGAc                                       >  1:398477/1‑251 (MQ=255)                                                                                                                                                                                                  gTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCa                                     >  1:322896/1‑250 (MQ=255)                                                                                                                                                                                                    tATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCaaa                                   <  2:400241/250‑1 (MQ=255)                                                                                                                                                                                                      tGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACg                                >  1:332088/1‑251 (MQ=255)                                                                                                                                                                                                      tGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACg                                >  1:256330/1‑251 (MQ=255)                                                                                                                                                                                                        cgcCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCa                                     >  1:315851/1‑244 (MQ=255)                                                                                                                                                                                                        cgcCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCa                                     <  2:315851/244‑1 (MQ=255)                                                                                                                                                                                                                ggATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGatca                      <  1:160570/251‑1 (MQ=255)                                                                                                                                                                                                                           tCCCAGATGCCGCTCTCTTTACGCAAGCGCATCGTGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCGGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGt            <  2:405628/250‑1 (MQ=255)                                                                                                                                                                                                                                    ccGCTCTCTTTACGCGAGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGCGACCCGAACGGCGTTTGGCGCCCAGACAGAGAAACCTGTACCTATGACGCAATCAATAGTATCTGCCTGCGCGCCTAAGGCTTCTTCCGGGCGCAGGTGAGTACCTTCAGCTAATAACCAGGCATCCACTCCCTGGATTAGAGAACCAAAAAGGTAAGGATCATCAATCAGGTTATGCTgccc  >  2:37882/1‑248 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           GAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTACGCAGGCGCATCGGGTGACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACC  >  NC_000913/3572758‑3573235

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