Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	790,832	Δ1 bp	coding (198/1047 nt)	galT ←	galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	790,832	0	G	.	100.0%	202.9 / NA	44	coding (198/1047 nt)	galT	galactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (26/18);  total (26/18)

CCCATGGGTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCC  >  NC_000913/790588‑791073                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  cccAGGGGTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTAt                                                                                                                                                                                                                                            <  1:300920/251‑1 (MQ=255)         tACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTc                                                                                                                                                                                                                                    >  1:111946/1‑251 (MQ=255)           cGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGAGTTTTACTGTGATCCGGCGAAAAGGAGATCACCCGGCTGGTGCCGCGAGTGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCTCGTCCGCCATCACAGCCGACACGTCCTTAGTGAAAACGTAAGTACAGGTGAAATCCGG‑TTTTTATCGCCTATCAc                                                                                                                                                                                                                                  >  2:357057/1‑251 (MQ=255)                 ccccAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGca                                                                                                                                                                                                                             <  2:176301/250‑1 (MQ=255)                   ccAGTTCTGCGGTTTGCTCCTGCCAGGATTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCAGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACAt                                                                                                                                                                                                                          >  1:344016/1‑251 (MQ=255)                         cTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTg                                                                                                                                                                                                                    >  1:457826/1‑251 (MQ=255)                                        gCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCa                                                                                                                                                                                                       <  1:236070/249‑1 (MQ=255)                                           aGGTTTTGCCGCTTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGCTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCCAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTg                                                                                                                                                                                                  >  2:224926/1‑251 (MQ=255)                                              ttttGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGAt                                                                                                                                                                                               >  2:305606/1‑251 (MQ=255)                                                 tGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGt                                                                                                                                                                                            >  1:187731/1‑251 (MQ=255)                                                 tGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAAGCTGGATCg                                                                                                                                                                                             <  2:171626/250‑1 (MQ=255)                                                               aaTGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATACGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGACAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TGTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACAcc                                                                                                                                                                              <  2:315558/251‑1 (MQ=255)                                                                 tgctgcAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTg                                                                                                                                                                            >  1:152010/1‑251 (MQ=255)                                                                                    ggCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTgc                                                                                                                                                         <  1:353877/251‑1 (MQ=255)                                                                                     gCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTgcg                                                                                                                                                        <  2:48710/251‑1 (MQ=255)                                                                                          gTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGAGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGccccc                                                                                                                                                   >  1:439603/1‑251 (MQ=255)                                                                                             ttACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGc                                                                                                                                                >  1:151762/1‑251 (MQ=255)                                                                                             ttACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGAGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGAGTAATCGGG‑TTTTTATAGCCTGTCACCCGCACCTTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGc                                                                                                                                                >  1:455690/1‑251 (MQ=255)                                                                                                  gtgATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAggg                                                                                                                                           <  2:51492/251‑1 (MQ=255)                                                                                                        cGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCtt                                                                                                                                     >  2:169003/1‑251 (MQ=255)                                                                                                                  cAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCgcggtgc                                                                                                                           <  1:305606/251‑1 (MQ=255)                                                                                                                            ggCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCACGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAAcc                                                                                                                  >  1:260840/1‑250 (MQ=255)                                                                                                                               tGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGa                                                                                                               >  1:54136/1‑250 (MQ=255)                                                                                                                                                      cGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTg                                                                                       <  1:23783/251‑1 (MQ=255)                                                                                                                                                          tCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAgc                                                                                   >  1:268340/1‑251 (MQ=255)                                                                                                                                                                           tCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATc                                                                   <  2:268340/250‑1 (MQ=255)                                                                                                                                                                             tGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTACCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAAc                                                                >  2:217933/1‑251 (MQ=255)                                                                                                                                                                                cgcATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACggg                                                             >  1:69081/1‑251 (MQ=255)                                                                                                                                                                                       ggCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCTCGGTGTAATCGGG‑TTGTTATCGCCTTTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGTGGGCTTAGCGCGGTGCGGTGAAACCAGAATCCAGTGCCCGGTGAGCGGGTTGTAGCAGCGATGTGGATGATCAGCGGGATTaa                                                        <  2:206925/249‑1 (MQ=255)                                                                                                                                                                                          gtgtCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGAGGCGATGTGGATGATCAACGGGATTAAATTGc                                                   <  2:88020/251‑1 (MQ=255)                                                                                                                                                                                            gtCAGACATCAAAGCCGCAAAGTCAGTAGTGAAAACGTAAGTCCCGGGGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACGCCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGGGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGt                                                 <  2:221755/251‑1 (MQ=255)                                                                                                                                                                                               aGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCAt                                              <  2:111946/251‑1 (MQ=255)                                                                                                                                                                                                      aaaGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGtt                                       <  2:444940/251‑1 (MQ=255)                                                                                                                                                                                                           cGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCtt                                   >  1:385790/1‑250 (MQ=255)                                                                                                                                                                                                             cAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAAt                                >  1:62667/1‑251 (MQ=255)                                                                                                                                                                                                                 gTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCggg                            >  1:348045/1‑251 (MQ=255)                                                                                                                                                                                                                 gTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCggg                            >  1:160858/1‑251 (MQ=255)                                                                                                                                                                                                                          gAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCtg                   <  1:243545/251‑1 (MQ=255)                                                                                                                                                                                                                                      cccGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTg        >  1:15306/1‑250 (MQ=255)                                                                                                                                                                                                                                       ccgttgtAATCGGT‑TTTTTATCGCCTGTCACCCGCACATTACCTCCGCAGAGGAAGCAAGCTGGATCGTTCGCGGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCATGGGCGCTTAGCGCGGTGCGGTGAAATTAGAATCCATTGCCCGGTGAGCGGGTTGTAGCTGCGATGTGGATGATCTACGGGATTAAATTGCGTCGTGGTTGTTCCTGATTCGGGagttcccttatgagggcgtgcc      <  2:269209/247‑22 (MQ=255)                                                                                                                                                                                                                                        cGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCGGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGTTGGCGTGACt     >  2:394851/1‑251 (MQ=255)                                                                                                                                                                                                                                         ggTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAACTTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtg    >  2:528/1‑251 (MQ=255)                                                                                                                                                                                                                                          gtgtAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgc   >  1:103605/1‑251 (MQ=255)                                                                                                                                                                                                                                           tgtAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgcc  >  2:45181/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  CCCATGGGTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCC  >  NC_000913/790588‑791073

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