Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,786,064	G→A	E269K (GAA→AAA)	gpmM →	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,786,064	0	G	A	100.0%	120.2 / NA	36	E269K (GAA→AAA)	gpmM	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/19);  total (17/19)

CGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAATATGCCCACGTTACTTTCTTCTTC  >  NC_000913/3785826‑3786291                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    cgagATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTg                                                                                                                                                                                                                         <  1:95534/248‑1 (MQ=255)    gATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCtt                                                                                                                                                                                                                      >  1:268760/1‑251 (MQ=255)       ggACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGGAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCg                                                                                                                                                                                                                   <  1:127035/251‑1 (MQ=255)              gATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGc                                                                                                                                                                                                            >  2:199834/1‑251 (MQ=255)                      ttGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCg                                                                                                                                                                                                    >  2:128960/1‑251 (MQ=255)                       tGGGATCGCGAAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCg                                                                                                                                                                                                    <  2:153600/250‑1 (MQ=255)                               cgTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTcg                                                                                                                                                                                           >  1:323618/1‑251 (MQ=255)                                       aaGCTTATTACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCCCTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAg                                                                                                                                                                                   <  1:216099/251‑1 (MQ=255)                                             aTGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGt                                                                                                                                                                              <  1:199834/250‑1 (MQ=255)                                                aCCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGtt                                                                                                                                                                          >  1:76771/1‑251 (MQ=255)                                                 cctgcggcctctGGCGCAGGGAGATTTCCAGGCCGATCCCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTa                                                                                                                                                                         <  2:366262/243‑1 (MQ=255)                                                    gctgACCCGGGCGCAGGGCGAGTTCCAGGCCGATCCCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTaa                                                                                                                                                                        <  2:231752/249‑1 (MQ=255)                                                       gACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCg                                                                                                                                                                   >  1:84420/1‑251 (MQ=255)                                                                  aGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATg                                                                                                                                                        >  2:98676/1‑251 (MQ=255)                                                                       gAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTACAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGAc                                                                                                                                                   >  1:278901/1‑251 (MQ=255)                                                                                gCCGATACCGCCGTTGCTGGTGTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATNcgc                                                                                                                                          >  2:369664/1‑251 (MQ=255)                                                                                  cggtACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATTACTTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGATGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGGCCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATAcgccg                                                                                                                                        <  2:157070/248‑1 (MQ=255)                                                                                  cGATACCGGCGTTTCTGGTTTGCAGGCTGCTGATGGTGGCGACGAAGATGAGGAATTCGTGAAAGCGACCTTTATCCGTGCTGAAGGGCAGACAGACGCGGCAATTGAAGACGGCGATGCACTGAGTTGCATGAACGTCCGTGTTTGCCGCACGCGGAAAAGCACGCGTGCGTGTGTGAGCGCTGATTGCGATGGGTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGcggaatgaggacgccg                                                                                                                                        <  2:255626/251‑16 (MQ=255)                                                                                        cgccgTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGAc                                                                                                                                   >  2:66371/1‑250 (MQ=255)                                                                                                              gctTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTac                                                                                                              <  1:253606/249‑1 (MQ=255)                                                                                                                               aaaTGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTa                                                                                           <  1:185160/251‑1 (MQ=255)                                                                                                                                     tgaATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCt                                                                                     >  1:359808/1‑251 (MQ=255)                                                                                                                                               aaaGCGCCCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTg                                                                           >  1:154845/1‑251 (MQ=255)                                                                                                                                                 aGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGa                                                                         <  2:359808/251‑1 (MQ=255)                                                                                                                                                          tATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGATGGCGaaa                                                                 <  2:278901/250‑1 (MQ=255)                                                                                                                                                           aTCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGATGGCGaaaa                                                                <  2:302113/250‑1 (MQ=255)                                                                                                                                                                    gAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTGCGGCGAGTGGATGGCGAAAAACGACaaaa                                                       <  2:71237/250‑1 (MQ=255)                                                                                                                                                                      aGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTGTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATAAATCGTGCTTCCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGACAGTGGTTAACGGCGATTTCGTGAAGCTGACCGAATAAGCAGCTGCCATCAAAACTGCGGTTGCTTACCCAACCGCATCCCTGGTTAACACATTCGGCGGGTGGATGGCAAAAACCGACacaccgc                                                    >  2:216808/1‑245 (MQ=255)                                                                                                                                                                       ggTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGGGTGGATGGCGAAAAACGACAAAACTCa                                                   <  1:211563/251‑1 (MQ=255)                                                                                                                                                                       ggTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTc                                                    >  1:231287/1‑250 (MQ=255)                                                                                                                                                                        gTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTc                                                    >  2:304840/1‑249 (MQ=255)                                                                                                                                                                                gACGCGGCAATGGAAGCCGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCTCCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTCAGTTGCGTAt                                          <  1:101184/251‑1 (MQ=255)                                                                                                                                                                                       cAATGGAAGACGGCGTTGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAGGATCACTCGTGCTTGCGTGAACGCTGATTTCGATGGCTTCGCACGTAAGAAAGTGGATAGCGGCGCTTTCGGGAGGCTGACCGAATACGCCTCTGACATCACAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTCATTTGCGTATTGccga                                    >  2:80/1‑250 (MQ=255)                                                                                                                                                                                                  ggCGATGCGCTGATTTTCATGAACTGTCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTGCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAata                        <  2:236203/251‑1 (MQ=255)                                                                                                                                                                                                  ggCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTAAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAat                         >  1:31797/1‑250 (MQ=255)                                                                                                                                                                                                                        aCTTCCGTTGTGTCCGCGCGCGTAAAATCACTCGTTCTTTCGTGAACGCTGATTTCCATGGCTTCGCGCGGAAGCAAGTGGTTAACGTCGATTTCGTGATTCTGACCGAATCCGCCGCTGACATCAAAACTGCGGTTTCTTACCCACCCGCTTCCCTGGTTAACCCCTTCGGCGAGTGTATGGCGAAAAACGACAAAACTCAGTTGCGTATTTCCTAACCCGAAAAATATGCCCACGTTCCTttcttcttc  <  2:295978/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    CGCGATGGACCGCGATAACCGTTGGGATCGCGTAGAAAAAGCTTATGACCTGCTGACTCTGGCGCAGGGCGAGTTCCAGGCCGATACCGCCGTTGCTGGTTTGCAGGCTGCTTATGCTCGCGACGAAAATGATGAATTCGTGAAAGCGACCGTTATCCGTGCTGAAGGTCAGCCAGACGCGGCAATGGAAGACGGCGATGCACTGATTTTCATGAACTTCCGTGCTGACCGCGCGCGTGAAATCACTCGTGCTTTCGTGAACGCTGATTTCGATGGCTTCGCGCGTAAGAAAGTGGTTAACGTCGATTTCGTGATGCTGACCGAATACGCCGCTGACATCAAAACTGCGGTTGCTTACCCACCCGCTTCCCTGGTTAACACCTTCGGCGAGTGGATGGCGAAAAACGACAAAACTCAGTTGCGTATTTCCGAAACCGAAAAATATGCCCACGTTACTTTCTTCTTC  >  NC_000913/3785826‑3786291

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