Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,150,418	A→G	intergenic (‑153/+29)	eptC ← / ← ppc	phosphoethanolamine transferase EptC/phosphoenolpyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,150,418	0	A	G	100.0%	157.7 / NA	45	intergenic (‑153/+29)	eptC/ppc	phosphoethanolamine transferase EptC/phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/16);  total (29/16)

AGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGATCGGCCATCAGATGGGAATCG  >  NC_000913/4150177‑4150657                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                agagTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACg                                                                                                                                                                                                                                          >  1:29582/1‑249 (MQ=255)       gagaAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTg                                                                                                                                                                                                                                  >  2:412545/1‑251 (MQ=255)             aaCCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGCAAGAACCTGTCTTGAAAAAATATCGACGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTTAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCaga                                                                                                                                                                                                                              >  1:324524/1‑249 (MQ=255)                 aaaaGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGCCTTACCCTGGACCTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTACTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTCGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAgg                                                                                                                                                                                                                         >  1:88991/1‑250 (MQ=255)                   aaGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAg                                                                                                                                                                                                                      <  2:184301/251‑1 (MQ=255)                     gTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTTAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATGTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGAt                                                                                                                                                                                                                    <  1:266106/251‑1 (MQ=255)                               gCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGCGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTAt                                                                                                                                                                                                          >  2:403404/1‑251 (MQ=255)                                  cAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTAc                                                                                                                                                                                                       <  2:228335/251‑1 (MQ=255)                                      ccAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCAt                                                                                                                                                                                                   <  1:575817/251‑1 (MQ=255)                                       cAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCAt                                                                                                                                                                                                   >  2:102770/1‑250 (MQ=255)                                                                      tGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGCAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCAtt                                                                                                                                                                    >  2:578882/1‑250 (MQ=255)                                                                         aCTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACg                                                                                                                                                                <  1:374607/251‑1 (MQ=255)                                                                          cTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAATTGCAGACAGAAATATATTGAAAACGAGGGAGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGCAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGc                                                                                                                                                               >  2:439502/1‑251 (MQ=255)                                                                           ttCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCt                                                                                                                                                              >  2:110310/1‑251 (MQ=255)                                                                                   aaTGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTGGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCGCGCGCAGAAGCACGAGGGTTGGCAGAAGAGGAAGATTAGCCGGTATTACTCATCCCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGAc                                                                                                                                                      <  2:465096/251‑1 (MQ=255)                                                                                           gTAAGAACCTGTTTTGAAAAAATATCGCCGAATGTAACGAAAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAAGGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGAt                                                                                                                                              <  2:550115/251‑1 (MQ=255)                                                                                           gTAAGAACCTGTCTTGAAAAAAGATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAATATTAGCCGGTATTACGCAGACCTGCCGCAATCCCGGCAATAGTGACCAGTAACGCTTGTTCGACGCGAGGAt                                                                                                                                              <  2:422792/251‑1 (MQ=255)                                                                                                     gTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTg                                                                                                                                    >  2:20929/1‑251 (MQ=255)                                                                                                            aaaaaTATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCcttct                                                                                                                             >  1:43515/1‑251 (MQ=255)                                                                                                             aaaaTATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCcttctt                                                                                                                            >  1:425685/1‑251 (MQ=255)                                                                                                                       cGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTg                                                                                                                  >  1:7672/1‑251 (MQ=255)                                                                                                                                  aCAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGt                                                                                                       <  1:110310/251‑1 (MQ=255)                                                                                                                                   cAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGAATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCATTCTTTTTCTGCCTGGCGGGAGAGGTg                                                                                                      >  2:402194/1‑251 (MQ=255)                                                                                                                                     aTTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTgca                                                                                                    <  2:351632/251‑1 (MQ=255)                                                                                                                                                       aaGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAAt                                                                                   >  1:250541/1‑250 (MQ=255)                                                                                                                                                          gTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGt                                                                               >  2:20957/1‑251 (MQ=255)                                                                                                                                                                  aGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCGTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCggg                                                                       <  1:75187/251‑1 (MQ=255)                                                                                                                                                                        tGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGgtg                                                                 >  2:38644/1‑251 (MQ=255)                                                                                                                                                                          aaTAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGgtgt                                                                >  1:5927/1‑250 (MQ=255)                                                                                                                                                                                   cagacagaAATATATTGAAAACGAGGGTGTTAGAACAGACGAATTTCAGACAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCg                                                      >  2:386923/1‑251 (MQ=255)                                                                                                                                                                                                  tGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAACCCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATCCCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCTAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGAct                                        <  1:361859/250‑1 (MQ=255)                                                                                                                                                                                                   gAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGActct                                      <  1:488793/251‑1 (MQ=255)                                                                                                                                                                                                        cGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTATGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAAt                                 >  1:448018/1‑251 (MQ=255)                                                                                                                                                                                                         gAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATc                                >  1:15367/1‑251 (MQ=255)                                                                                                                                                                                                         gAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCGCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGTTATTACGCATCCCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATc                                <  2:496488/251‑1 (MQ=255)                                                                                                                                                                                                          aGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATcc                               >  2:464311/1‑251 (MQ=255)                                                                                                                                                                                                            ggTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATcaac                             >  2:263578/1‑248 (MQ=255)                                                                                                                                                                                                             gtgtTATAACAGAATGATTTCAGAAACCCCTCGCGCAGACGCGCGAGGGTTCGCAAAAGAGGAAGATTAGCAGTTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCAGTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCATCAACTCGGCTTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTATCTGAATAGACTCTGCAATCCACg                            <  2:332329/251‑1 (MQ=255)                                                                                                                                                                                                               gtTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTCCGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGc                          >  2:292312/1‑251 (MQ=255)                                                                                                                                                                                                               gtTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGc                          >  2:75987/1‑251 (MQ=255)                                                                                                                                                                                                                    aaCAGAAGTATTTCAGAAAACCCTCGCGCAGCAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAGTCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGAGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGc                          >  1:225136/1‑246 (MQ=255)                                                                                                                                                                                                                     aaaGAATTATTTCAGAAAACCCTCGCGCAGAAGCACGGGGGTTGGCAGAAGAGGAAGATGAGCCGGTATCACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGATCg                    <  2:387419/249‑1 (MQ=255)                                                                                                                                                                                                                      cAGAAGTATTTCAGAAAACCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCACTCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGc                          >  2:375129/1‑244 (MQ=255)                                                                                                                                                                                                                                    aaaaCCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGATCGGCCATCAGATGGGa      >  2:101132/1‑250 (MQ=255)                                                                                                                                                                                                                                       aCCCTCGCGCAGAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGATCGGCCATCAGATGGGAATCg  >  1:142119/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                AGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACGCTTGTTCGACGCGAGGATCCGGTTCCTGGCCTTCTTTTTCTGCCTGGCGGGAGCGGTGCAGCAACTCGGCCTGCAATACGTTCAGCGGGTCGGTGTAAATATTCCGTAGCTGAATAGACTCTGCAATCCACGGCAGATCGGCCATCAGATGGGAATCG  >  NC_000913/4150177‑4150657

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