Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	21,619	A→G	R71R (AGA→CGG)	ribF →	bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	21,619	0	A	G	100.0%	195.4 / NA	68	R71R (AGA→CGG)	ribF	bifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/29);  total (39/29)

ACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATAT  >  NC_000913/21375‑21858                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                tcTCGGAAGATTTGCGCTGTTATGAACCAGACATGAAGCGGATACGCGGCAGACATAAGCTCAGCCAGGCCCCGCAAGAAGGTTGTGTGCTGGCTATGGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGTGCTGTTACAGGGCTTGCAGGAAGAAGCGCGCACGCGCATCTGACCGGTGATGGTGATGCTTTGTGAACCGCAACCACTGGAACCGTTTGCGACCGATAAAGCCCCGGCACGGCTGAc                                                                                                                                                                                                                                            <  2:631940/249‑1 (MQ=255)      gACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCCACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGCCCCGGCt                                                                                                                                                                                                                                      >  2:644009/1‑251 (MQ=255)             ttCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGa                                                                                                                                                                                                                                <  1:290507/250‑1 (MQ=255)             ttCACTGTTTTGAGCCAGACATGAAGCTGATACGAGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGAGTGTGATGACTATTGGTAATTACGACGGCGTGCATCGCGGACATCGCGCGCTGTTAAAGGGCTTGCCGGAAGAAGGGCGCAAGCGCAACTTACCGGTGAAGGTGATGCTTTTTGAACCTCAACACCTGGAACTGTTTGCTACCGATAAAGCACCGGCACGGCTGACCCGACTGCggtg                                                                                                                                                                                                                                >  2:477049/1‑248 (MQ=255)                    ttttGAGCCCGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGAGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGGTGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGc                                                                                                                                                                                                                        >  2:330251/1‑251 (MQ=255)                    ttttGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGACGGGCGCAAGCGCAACTTACCGGTGATGGGGATGCTTTTTAAACCTCCACCACTGGAACTGTTTGCTCAATATAAAGCCCCGGACCGGCAGACCCGGCGGCGGGCAAAACTGc                                                                                                                                                                                                                        >  2:12132/1‑251 (MQ=255)                        gAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGAGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACAGGTGATGGTGATGCTTTTTGAACCTCAACACCTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCGGCGGGAAAAACCGCGta                                                                                                                                                                                                                     >  2:652962/1‑249 (MQ=255)                           ccAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATCAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCt                                                                                                                                                                                                                 >  1:676702/1‑251 (MQ=255)                            cAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCtt                                                                                                                                                                                                                <  1:417861/251‑1 (MQ=255)                              gACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGc                                                                                                                                                                                                              <  1:534156/251‑1 (MQ=255)                              gACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGc                                                                                                                                                                                                              >  2:5269/1‑251 (MQ=255)                                 aTGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGCTAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGc                                                                                                                                                                                                              >  2:199641/1‑248 (MQ=255)                                                  catacataATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACg                                                                                                                                                                                          >  1:448358/1‑251 (MQ=255)                                                           tctcAGCCAGTCCCCGCAAGCAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCACAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAAGTGCGTTAACTTGCAGAGGGTGGCGTTGAGTGAGTGCTgtgcg                                                                                                                                                                                 >  2:254638/1‑251 (MQ=255)                                                            ctcAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGAAGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGGGATGCTTTTTGAACCTCAAACACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGACACACTGCGTAACGTTGCAGCGGGTGGCGTTGGTTACGTGCTgtgcgt                                                                                                                                                                                >  2:628227/1‑251 (MQ=255)                                                             tcAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTgtgcgtg                                                                                                                                                                               >  1:570140/1‑251 (MQ=255)                                                              cAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTgtgcgtgc                                                                                                                                                                              >  1:158605/1‑251 (MQ=255)                                                               aGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGtgcgtgcg                                                                                                                                                                             >  1:261268/1‑251 (MQ=255)                                                                    gtcccctagataagGGTGTGTGCTGACTCTTGGTAATTTCGACGGCGCGCATCGCGGGCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTTGAACTGTTTGCTACCGATAAAGCCCTGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGATTGTGGCGTTTATTACGTGCGGTTCGTGCGTTTct                                                                                                                                                                        <  2:404067/240‑2 (MQ=255)                                                                      ccccGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGCCGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACGGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGATGACCCGGCTGCGGGCACAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTGAGTGCTGTGAGTGCGTTTCGAc                                                                                                                                                                      >  2:493111/1‑251 (MQ=255)                                                                           cAAGAAGGGTGAGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGCAGAAGGACGCAAGCGCAACTTACCGGAGATGTTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGACAACATGCGTTACCTTGCAGAGTGTGGAGGTGATTACGTGCTGGGCGTGCGTTTCGAcggtcg                                                                                                                                                                 >  2:522104/1‑246 (MQ=255)                                                                            aagaagGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATAGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCgg                                                                                                                                                                >  2:267595/1‑250 (MQ=255)                                                                             catgtacgatctggtGACTATTGGGAATTTAGACGGCGTGCAGCGCGGGCATCGCGCGCGGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGAGGGTGATGCTTTTTGAACCTCAACCACTGGAACTGGTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCggt                                                                                                                                                               <  2:298799/237‑3 (MQ=255)                                                                                agGGAGTGTGCTGACTATTGGTCATTTCGACGGCGTGCATCGCGGACATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCGGCGGGAAAAACTGCGTTCCCTTGCAGAGAGTGGCGTTGGTTACGTGATGTGCGTGCGATTCGACAGGCGTTTcg                                                                                                                                                            >  2:456149/1‑251 (MQ=255)                                                                                      tgtgCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGt                                                                                                                                                      >  1:17257/1‑251 (MQ=255)                                                                                         gCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTTTGCGTGCGTTTCGACAGGCGtt                                                                                                                                                               >  1:298799/1‑239 (MQ=255)                                                                                           tGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTCCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGCTACAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGCTTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAAcc                                                                                                                                                 >  1:630361/1‑251 (MQ=255)                                                                                            gactatgtgTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCCCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAAc                                                                                                                                                  <  1:608703/241‑1 (MQ=255)                                                                                                   ggTAACTTCGATGGTGTGCATCGCTGTCTTCTCGCGCTGTTACAGTGCTTTCAGGAAGTAGGGCGCTAACGCAACGTACCGGAGATGGTGATGCTTTTTGAACCTTAATCACTGGCACGGTTCGCTACCTATAAAGCCCCAGCACGGCTGTTCGGGTTGTGGGACACACTGCGTTACCTTGCAGAGTGTGGCTTTGATTACGTGCAGTGCGTGCGTTTCGACAGGCGTTTCGCGGCTTTAACCGCGCaaaa                                                                                                                                         <  2:151248/251‑1 (MQ=255)                                                                                                    gTAATTTCGCCGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAGAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGGGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTaa                                                                                                                                                   >  2:620319/1‑240 (MQ=255)                                                                                                    gTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAGAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTaa                                                                                                                                                   <  1:620319/240‑1 (MQ=255)                                                                                                     tAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGCTGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGCTTACGTGCTTTGCGTGCGTTTCGACAGGTGTTTCGCGGCGGTAACCGCGCAAAAtt                                                                                                                                       >  2:459223/1‑251 (MQ=255)                                                                                                      aaTTTCGTCGGCGTTCATCGCGGTCATCGCGCTCTGTTACAGGGCTTTCAGGAAGAAGGGCGCAAGCGCAACTTTCCGGTGATGTTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGCTAAAGCCCCGGCCCGGCTGACCCGGCTGCGGGAAAATCTGCGTTTCCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAAtt                                                                                                                                       <  2:557464/250‑1 (MQ=255)                                                                                                             aCGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGc                                                                                                                               <  1:256094/251‑1 (MQ=255)                                                                                                               ggCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCgg                                                                                                                             >  1:311823/1‑250 (MQ=255)                                                                                                                  gTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGAGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGAtct                                                                                                                          >  1:677287/1‑251 (MQ=255)                                                                                                                    gCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTTTATGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGAGTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGAtctt                                                                                                                         <  2:80415/250‑1 (MQ=255)                                                                                                                     cATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGAtcttct                                                                                                                       >  1:147778/1‑251 (MQ=255)                                                                                                                               aTCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTGGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGctt                                                                                                             <  2:56488/251‑3 (MQ=255)                                                                                                                                       cTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGt                                                                                                     >  1:130454/1‑251 (MQ=255)                                                                                                                                        tGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGt                                                                                                     >  1:525268/1‑250 (MQ=255)                                                                                                                                            tcAGGGCTGGCAGGAAGAAGGTCGCAAGCGCAATGTTCCGGTGCTGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCCACCGCTAAAGCCCCGGCACGGCTGACCCGTCTGCGGGAAAATCTGCGTTACCTTGCAGAGGGTGGCGTTGATTACGGGCTGTGCGTGCGTTTCGACAGGTTTGTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTGCTGGTGAAGCGTTTGCGTGTaaaa                                                                                                 <  2:406787/249‑1 (MQ=255)                                                                                                                                            aCAGGGCTTGCAGGAATAATGGCGCAAGAGCAACTTTCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTaaaa                                                                                                 <  2:590253/250‑1 (MQ=255)                                                                                                                                                    tGCAGGAAGAAAGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGcc                                                                                        >  1:152095/1‑251 (MQ=255)                                                                                                                                                         gaagaaGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAgg                                                                                   <  1:592882/251‑1 (MQ=255)                                                                                                                                                                 gcgcaagcgcaaCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGACGCTTTTGCGCGTAAATTTTCTTGCCGTGGGTGATGAtt                                                                           >  2:248967/1‑251 (MQ=255)                                                                                                                                                                     aagcgaaaCTGACCGTTGATGGTGATGCTTTCTGAACCACAACCACGGGAACTGTTTGCGACCGATCAAGACCCGGAACGGCTGACCGGGCGGCGGGAAAAACTGCGTTACCTTGTATAGTGTGGCGGTGATTACGTGCTGGTCGTGCGCTGCGATAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCACCATTGATCTTCTGGTGAAGCTTTTGCGCTTAAAATTTCTCGCCTTAGTTGATGAtttgcg                                                                       <  2:76007/251‑4 (MQ=255)                                                                                                                                                                        cgcaaCGTACCGGTTATGGTGATGCTGTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAAGTTCTTGCCGTAGGTGATGATTTCCGCtt                                                                    <  1:465965/251‑1 (MQ=255)                                                                                                                                                                         gcaaCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTAACTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGAGATGATTTCCGCttt                                                                   >  1:654642/1‑251 (MQ=255)                                                                                                                                                                                aCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTg                                                            <  1:562791/251‑1 (MQ=255)                                                                                                                                                                                         ggtgatgCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGaa                                                    <  1:5269/250‑1 (MQ=255)                                                                                                                                                                                            gatgCTTTGTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACCGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGc                                                 <  1:199722/250‑1 (MQ=255)                                                                                                                                                                                                   tttGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACTGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCt                                          <  1:307769/250‑1 (MQ=255)                                                                                                                                                                                                       aaCCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCACGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGAGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGCTCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGTTGATGATTTCCTCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGtt                                      >  2:254063/1‑250 (MQ=255)                                                                                                                                                                                                             aaCCACTGGACCTGTTTGCTACCGATAAAGCCCCGGCCCGGCTGACCCGGTTGTGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAg                               <  2:448358/251‑1 (MQ=255)                                                                                                                                                                                                                    ggAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAACCTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGAGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGAAAAATTTCTTGCCGTAGGGGATGATTTCCGCTTTGGAGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTg                        >  1:617233/1‑251 (MQ=255)                                                                                                                                                                                                                              gCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGCTCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTACGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGctgggaatc              >  2:177250/1‑243 (MQ=255)                                                                                                                                                                                                                               cTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATAc             >  2:133264/1‑251 (MQ=255)                                                                                                                                                                                                                                 aCCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTTCCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATAc             <  2:263403/249‑1 (MQ=255)                                                                                                                                                                                                                                 aCCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTGCAGCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACgg           <  2:419024/251‑1 (MQ=255)                                                                                                                                                                                                                                  ccGATAAAGCCCCGGCACGGCTGGCCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGc          <  2:17257/251‑1 (MQ=255)                                                                                                                                                                                                                                   cGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCt         >  1:142951/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGCTCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAACTGGCATGGAATACGGCttgg      >  2:207326/1‑247 (MQ=255)                                                                                                                                                                                                                                         aGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGata   <  2:261268/251‑1 (MQ=255)                                                                                                                                                                                                                                         aGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGGTGATTACGGGCTGTGCGTGCGGTTCGACAGGCGTTTCGAGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGata   <  1:611378/251‑1 (MQ=255)                                                                                                                                                                                                                                          gCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGatat  >  1:366516/1‑251 (MQ=255)                                                                                                                                                                                                                                          gCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGatat  <  1:644009/251‑1 (MQ=255)                                                                                                                                                                                                                                          gCCCCGGCACGGATGACCCGGCTGCGGGAAAAACTGAGTTACCTTGCAGAGTGTGGCGATGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCTCAAAATTTCATCAGCGATCTTATGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATAATTTCCGCTTTGGCGCTGTGCGTGAAGCCGATTTCTAGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGatct  >  2:533149/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                ACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATAT  >  NC_000913/21375‑21858

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