Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,337,394	A→G	S522G (AGC→GGC)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,337,394	0	A	G	100.0%	118.7 / NA	40	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/20);  total (20/20)

AGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAG‑AAATTGCCCGTGCGGTAGAACAAGTC  >  NC_000913/1337154‑1337633                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 aGTGTGCGGATGGCAGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCTTAACGGCGATCAAAAAAGGCGAGTTAAc                                                                                                                                                                                                                                        <  2:281285/251‑1 (MQ=255)    gtgctgagggctgcagcactGCTGGAGAACAAAGCCGTAACTATGGGCCGCCAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGGTCGAAAGGCGATTCTGAGTATCTGGCAAAAGCGAAACGGCCACCGTATCTCGACGAACGGGGGTTCAACCTTGTGTGAGACGGTTTTACCACATGTATTGGTCACTCTGGGCCGCGGCCCGCTCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCgc                                                                                                                                                                                                                                     <  2:107933/232‑2 (MQ=255)    gtgCGGATGGCCGCAGGCTTGCTGGCGAAAAAGGCCGTAACTCTGGGCCACAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGGTGCTGATAATCTGGCAAAAGCGAAGCTGACACCGTATCTCGACGAACTGGGTTTTAACCTTGTGGGATGCGGGTGTACCACCTGTATTGGTAACTTTTGTCCGCTTCCCGAGCCTATCGATACTGCAATCAAAAAAGGCGATTGAACCGt                                                                                                                                                                                                                                     <  2:431716/251‑1 (MQ=255)      gCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcg                                                                                                                                                                                                                                   <  2:22812/251‑1 (MQ=255)       cTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAAATCTGGGCCTCACGCGGCAACCATGGAACAAAGCGTTGCTGGCACCGGATTCGAACGTCGCTTATGCTTATCTGGCAACCGCGAAAATGACAACGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcgg                                                                                                                                                                                                                                  >  1:38512/1‑251 (MQ=255)          aTGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAACAAGGCGATTTAACCGt                                                                                                                                                                                                                                     >  2:118063/1‑245 (MQ=255)          aTGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCTATTTAACCGt                                                                                                                                                                                                                                     <  1:118063/245‑1 (MQ=255)                  aGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCTGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGt                                                                                                                                                                                                                       >  1:253914/1‑251 (MQ=255)                  aGGCTTGCTGGCGAAAAAAGCCGTAACGCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGt                                                                                                                                                                                                                       <  2:95191/251‑1 (MQ=255)                    gCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTTGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTcc                                                                                                                                                                                                                     <  2:158024/251‑1 (MQ=255)                        gCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCa                                                                                                                                                                                                                 <  2:263633/251‑1 (MQ=255)                             ccaaaaaaGCCGTAACTCGTGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTGGTGGGATACGGTGGTACCACCTGTAGTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGt                                                                                                                                                                                                            <  2:307311/249‑1 (MQ=255)                                      ccGTAACACTGGGCCTCAAGCGGAAACCATGGGTCAAAGCGTCGCTGCCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTACCCTTGTGGGATACGGTTGTACCACCTGTCTTGGTAACTCAGGGCCGCAGCCCGATCCTATCGCAACGGCAATCACAAAAGGCGATTTTACCGTCGGTGCGGTGCTGTCCGGCAACCGTCAATTTGa                                                                                                                                                                                                    >  2:395350/1‑250 (MQ=255)                                           aCTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGGATTGGTAACTCTGGTCCGCTGCCCGATCCTAGCGAAACGGCAATCAAAAAAGGCGAGTTAACCGTCGGTGCGGTTCTGTCCGGCAACCGTAACTTTGAAGGCCg                                                                                                                                                                                              <  2:232151/251‑1 (MQ=255)                                                gggCCTCAAACGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatc                                                                                                                                                                                          <  1:251228/250‑1 (MQ=255)                                                        aGCGGCCACCATGGACCAAAGCGTCGCTGGCACCGGGTTCGAAAGCCGTTTCTGATTATCTGGCAACAGCGAAAATGAAACCCTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGACACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTg                                                                                                                                                                                 >  2:411725/1‑251 (MQ=255)                                                                         aaaGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActggct                                                                                                                                                                >  2:93371/1‑251 (MQ=255)                                                                          aaGCGTCGCTGGCACCGGCTACGAAAATCGTTTATGATTATCTGGCCACAGCGAAACTGACCACGTATCTCGCCGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGCAACTCTGGCCCGCTGCCCGACCATATCGAAACGGCAATCAAAAAAGGCGATTTCACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActggctg                                                                                                                                                               >  2:448141/1‑251 (MQ=255)                                                                                cGCTGGCGACGGGTACGAAAGTCGTTTCTGATTATTTTGCAAAAGCGAAACTGACATCGCATTTCGTCGAACTGGGGTTTAACCTTGTGGTATACGGTTTGACCACCTGTATTTGTAACTCTGGGCCGCTGCTCGATGCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTTCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCAGCCGCTGGTTAAATCTTACTGGCTGGCCTc                                                                                                                                                          <  2:326644/250‑1 (MQ=255)                                                                                       aCCGGGATCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGACACTGACACCGTATCGCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGCAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGACGGCCGTATCCATCCGCTGGTTACAACTAACTGGCTGGCCTcgccgccgc                                                                                                                                                  >  2:214072/1‑251 (MQ=255)                                                                                                                         aaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGaa                                                                                                                <  1:17713/251‑1 (MQ=255)                                                                                                                           aGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAGCTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGGAAAATGAat                                                                                                               >  2:216555/1‑250 (MQ=255)                                                                                                                            gCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGgaatatgaatat                                                                                                             >  2:337215/1‑251 (MQ=255)                                                                                                                                    gACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTgg                                                                                                     <  1:344945/251‑1 (MQ=255)                                                                                                                                          gTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTg                                                                                               >  1:122697/1‑251 (MQ=255)                                                                                                                                          gTATCTCGACGAACAGGGGTTTAACCATGTGGAATACGGTTGTAACACCTGTATTGGTACCTCTGGGCAGCTGCACGATCCTATCGAAACGGAAATCACAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTACGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAAAATGAATATCAACCCGGCttttg                                                                                               >  2:298833/1‑248 (MQ=255)                                                                                                                                                cgacgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCt                                                                                          <  1:337215/250‑1 (MQ=255)                                                                                                                                                   cgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCg                                                                                      >  1:139464/1‑251 (MQ=255)                                                                                                                                                                         ggATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGGGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATcc                                                                >  1:256504/1‑251 (MQ=255)                                                                                                                                                                                   gTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTg                                                                                               >  1:178096/1‑210 (MQ=255)                                                                                                                                                                                   gTACCACCTGTATTGGTAACTCGGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTg                                                                                               <  2:178096/210‑1 (MQ=255)                                                                                                                                                                                                                gtGCCCGGTCCTATCAAAACGACATACAAAAAAGGCGATTTAACCGGCGGTGCGGTGCTGTTCGGCAACCTGAGCTTTGATGGCCGTATCCAGCCGCTGGTTAAACCTAACGGGCTGGCCTCGCCGCCGCGGGTGGTTGCCTATGCGCTGTCGGGAAATATGAGTATCATCTTGGCTTCTGAGCCGATCGGCCATGATCGCAAAGGCGATCCGGTTTAGCTTAAAGATATCTGGCAATCTGCTCAAG‑AAAt                        <  2:284267/250‑1 (MQ=255)                                                                                                                                                                                                                    ccGATCCTATCGAAACGGCAATCAACAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGACATATGAATATCCACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGAGAACCGGATTATCTGTAAGATTTCTGGCCCTCGGCCCAAA‑AAATTGcc                    >  2:207027/1‑251 (MQ=255)                                                                                                                                                                                                                       aTCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGTGTCTAGCGGCCATGATCGCAAAGGCGATCCGGTGTATCTGAAAGATATCTGGCCATCGGCACAAG‑AAATTGCCCGt                 <  1:298833/251‑1 (MQ=255)                                                                                                                                                                                                                               gAAACGGCAATCAAAAAAGGCCATTTAACCGTCGGTGCGATCCTGTACGGCAACCGTAACTTTGAAGGCCGTATCAATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCAGGCGGGAAATATGAATACCCACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCCTCGGCACCAG‑AAATTGCCCGTGCGGt            >  1:422698/1‑248 (MQ=255)                                                                                                                                                                                                                                aaaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAG‑AAATTGCCCGTGCGGTAGa         >  1:341022/1‑250 (MQ=255)                                                                                                                                                                                                                                 aaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAata                                                                                                              <  2:28476/149‑1 (MQ=255)                                                                                                                                                                                                                                 aaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAata                                                                                                              >  1:28476/1‑149 (MQ=255)                                                                                                                                                                                                                                    ggCAATCAAAACAGGCGATTTACCCGTCGGTGCGGTGCTGTCCGGCACCCGTACCTTTGCAGGCCGTATCCATCCGCTGGTAAAAACTAACTGGCAGGCCTCGCCGCAGCTGGAGGTTGCCTCTGCGCTGGCGGGAAAAATGAATATCAACCTGGCTTCTGAGCCTATCGGACATGATCGCAAAGACCATCCGGTTTTTCTCACAGATAACCGGCCCTCCGCCCAAGGAAATCGCC‑GGGCGGTAGAACaat    >  2:417266/1‑250 (MQ=255)                                                                                                                                                                                                                                      cAATCAAAAAAGGCGATTTAACCTTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCGGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTTGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAG‑AAATTGCCCGTGCGGTAGAACAAGtc  <  1:249960/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 AGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAG‑AAATTGCCCGTGCGGTAGAACAAGTC  >  NC_000913/1337154‑1337633

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