Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,486,205	C→A	T29K (ACG→AAG)	crp →	DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,486,205	0	C	A	100.0%	171.4 / NA	49	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (30/19);  total (30/19)

TACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACC  >  NC_000913/3485958‑3486450                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         tACAGTAGTACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGAAGTACAGTTGATAGCCCCTTCCAAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGAGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACACGACTCTCGAATGGTTCTTGACTAATTGCCACATTCATAAGTACCCATCCAACAGCAAGCt                                                                                                                                                                                                                                                       >  2:135554/1‑251 (MQ=255)    aGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGCCTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGCGCAAGCTTAt                                                                                                                                                                                                                                                    >  2:375994/1‑251 (MQ=255)    aGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGAAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACATACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTAt                                                                                                                                                                                                                                                    <  1:139165/251‑1 (MQ=255)     gTAATACATTGCTGTAATGCATGAATGCAAAGGACGTCAAATTACCGTGCAGTAAAGTTGATAGCCCCTTCCCAGGTAGCAGGAAGCATATTTCGGCAATCCAGAGACACCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATATGTACCCATCCAAGAGCAAGCTTAtt                                                                                                                                                                                                                                                   >  2:141377/1‑251 (MQ=255)              tGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGt                                                                                                                                                                                                                                          >  1:79834/1‑251 (MQ=255)                aTGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTATCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTg                                                                                                                                                                                                                                         >  2:353435/1‑250 (MQ=255)                 tGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATACAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGAAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCCACAGGGTGa                                                                                                                                                                                                                                        >  2:81021/1‑250 (MQ=255)                        cATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGTGCAAGCTTATTCACCAGGGTGAAAAAGCgg                                                                                                                                                                                                                                >  2:420834/1‑251 (MQ=255)                             aTGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAc                                                                                                                                                                                                                            <  1:445446/250‑1 (MQ=255)                                 aaaGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGt                                                                                                                                                                                                                       >  1:31127/1‑251 (MQ=255)                                    gggcGTCACATTACCGTGCAGTAAAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTCCGGAAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAACCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAACCCGCAACCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAACCGCTGtact                                                                                                                                                                                                                    <  2:407143/248‑1 (MQ=255)                                    ggACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGCAAAAGCGGCAACGCTGgact                                                                                                                                                                                                                    >  2:61930/1‑251 (MQ=255)                                       cGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGAACCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGATGtactac                                                                                                                                                                                                                  >  1:277563/1‑250 (MQ=255)                                        gTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCCGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGGGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAt                                                                                                                                                                                                                >  2:225128/1‑251 (MQ=255)                                                 ccGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATAGTTAAA‑gg                                                                                                                                                                                                      >  2:317398/1‑251 (MQ=255)                                                      cAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTACA‑GGCTCtg                                                                                                                                                                                                 >  2:40082/1‑251 (MQ=255)                                                        gTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCtgtg                                                                                                                                                                                               <  1:9454/251‑1 (MQ=255)                                                                gATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCt                                                                                                                                                                                       <  1:143230/251‑1 (MQ=255)                                                                                gTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAagag                                                                                                                                                                       <  1:326013/251‑1 (MQ=255)                                                                                  aGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAggg                                                                                                                                                                     >  2:211244/1‑251 (MQ=255)                                                                                     gggAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTa                                                                                                                                                                   >  1:142984/1‑250 (MQ=255)                                                                                       gaggCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAg                                                                                                                                                                <  2:167562/248‑1 (MQ=255)                                                                                                  aggaaaTCCAGAGACAGCGGCGTTATCTGGCTGTGGAGAAAGCTTAGAACAGAGGATAACCGCGCATGGTGCTTGGCAAAACGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACAGTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGGACTACAGCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCct                                                                                                                                                      <  2:103370/246‑1 (MQ=255)                                                                                                   ggCAAGCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCctct                                                                                                                                                    <  2:79834/251‑1 (MQ=255)                                                                                                       aTCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCt                                                                                                                                                 >  1:292515/1‑250 (MQ=255)                                                                                                          cAGAGACAGCGGCTTTATCTGGCTCTGGAGAAAGCTTAGAACAGAGGATAACCGCGCATGGTGCGTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCt                                                                                                                                             <  2:51104/251‑1 (MQ=255)                                                                                                                aagaggcgTTATATGGCTCCGGCCAAAGCTTATCAAAGCGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACCGACCCGACTCTCGACTGGGTCTTGTCTCATTGCCACATTCATAAGTACCCATCCCACAGCCAGCTTATTCATCTGGGTGAAAAAGCGGCAACGCAGTACACCATAGTTTAA‑GGCTCTGTGGCAGCGCTGATCAAAGACCAAGAGTGCAAAGAAAGGATCCTCTCCTATCTGAATcc                                                                                                                                       >  2:298760/5‑250 (MQ=255)                                                                                                                   cggcgTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAggg                                                                                                                                    >  1:138663/1‑251 (MQ=255)                                                                                                                                            cTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CAGGGCCTg                                                                                                          >  2:310777/1‑251 (MQ=255)                                                                                                                                                 aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTg                                                                                                      >  1:252825/1‑250 (MQ=255)                                                                                                                                                                           tGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATg                                                                            <  2:142984/249‑1 (MQ=255)                                                                                                                                                                           tGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATg                                                                            <  1:420834/249‑1 (MQ=255)                                                                                                                                                                             gCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGt                                                                         >  1:198094/1‑250 (MQ=255)                                                                                                                                                                              cAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTAc                                                                       <  1:198951/251‑1 (MQ=255)                                                                                                                                                                                aaCCGCAAACAGACCCGACTCTCGAATGGGTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGt                                                                     <  2:252825/251‑1 (MQ=255)                                                                                                                                                                                 aCCGCAAACAGACCCGACTCTCGAATGGATCTTGTCTCATTGCCACATTAATAAGTACCAATACAAGCGCAAGCTTATTCACCTGGGTGAAAAAGCGGAAACGCTGAACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATAAAAGACGAAGAGGGTAAAGAAATGATCCTATCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAAGGTAGCGCATGGGAACGTg                                                                    >  1:107300/1‑251 (MQ=255)                                                                                                                                                                                       aaCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGTGCGaaaa                                                              <  1:211244/251‑1 (MQ=255)                                                                                                                                                                                          aGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCg                                                           >  1:260443/1‑251 (MQ=255)                                                                                                                                                                                                        aaTGGTTCTTGTCTCATTGCCACATTAATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCt                                             >  2:362092/1‑251 (MQ=255)                                                                                                                                                                                                        aaTGGCTCTAGTCTCATTGCCACATTCATAAGTACCCATCCACGAGCAAACTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGAGATGATCAAAGACGAAGCGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATATTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTATCGCATGGGTACGTGCGAACACCGCCTGTGAAGGGGCt                                             >  2:362098/1‑251 (MQ=255)                                                                                                                                                                                                               cTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAAttt                                      >  1:254408/1‑251 (MQ=255)                                                                                                                                                                                                                      catggCCACATTCAAATGTACCCACCCAAGAGCAAGCTTATTCACAAGGGTGAAAAAGCGGAACCGCTGTACTACAGCGTAAAA‑GGCTTTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGACCCTCTCCAACCTGAATCAGGGTGATTTTATTGGCGAC‑CGGGGCTTTTTTGAAG‑AGGGCAAGAAACGAGGCGTATGGGTACGTGGGAAATCCGCCTGTGAAGTGGCTGAAATTTCGTACaa                               <  1:31559/247‑1 (MQ=255)                                                                                                                                                                                                                                     aaGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAAtt                >  1:142062/1‑251 (MQ=255)                                                                                                                                                                                                                                       gtccccATCCAAGAGCAACCTTATTCCCCAGGGTGAAAAAGGGGAAACGCGGAACTAAACCGTAAAA‑GGCCCTGTGGCAGGGCTGAACAAAGACGAAGGGGGAAAAAAAAGAACCCTCCCCTATAAGAATCGGGGGAATTTTTTTGG‑GAAGCTGGGCTTGTT‑GAAGGAGGGCCGGAAACGTACCGTATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGa              <  1:173751/248‑1 (MQ=255)                                                                                                                                                                                                                                         aCCCATCCAAGTGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCGGGAACGTAGCGCATGGGTACGTGCGAAAATCGCCTGTGAAGTGGCTGAAATGTCGTAAAAAAAGTTTCGCCAATTGAtt            >  2:273762/1‑251 (MQ=255)                                                                                                                                                                                                                                           caaTCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAT‑GTTAAAAGGCTCTGTGGCGGTGTTGAACAAGGACGAAGAGGTTAAAGAAATGGTCCTCTCCTAGCGGAAGCAGGGTGTTTTTATTGGCGAC‑CGGGGCCTGTTTGAGG‑AGGGCCAGGAACGAAGCGCATGGGTACGTGCTAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCa          <  2:87297/249‑1 (MQ=255)                                                                                                                                                                                                                                               ccAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGt       >  1:127857/1‑250 (MQ=255)                                                                                                                                                                                                                                                   gagCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAcc  >  2:156959/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gagCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAcc  <  1:273130/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         TACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAA‑GGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAA‑CTGGGCCTGTTTGAAG‑AGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACC  >  NC_000913/3485958‑3486450

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