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%	151.5 / NA	48	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/30);  total (18/30)

TAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCACA  >  NC_000913/1337147‑1337639                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      tAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCg                                                                                                                                                                                                                                                    >  1:337268/1‑251 (MQ=255)             gCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                                                                                                                                                                                                         <  1:452026/250‑1 (MQ=255)                   ggCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAAGCAAAAAAGGCGATTTAACCGTcggtgcgg                                                                                                                                                                                                                                  <  1:157238/251‑1 (MQ=255)                       gCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCt                                                                                                                                                                                                                              <  1:129262/251‑1 (MQ=255)                        cAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTACCCGTCGGTGCGGTGCTg                                                                                                                                                                                                                             >  2:312098/1‑251 (MQ=255)                              tgcgggCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGGTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGGATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCgg                                                                                                                                                                                                                        <  2:445728/246‑1 (MQ=255)                                  ggCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAc                                                                                                                                                                                                                    <  1:424421/250‑1 (MQ=255)                                                                        cATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaact                                                                                                                                                                             <  1:251817/251‑1 (MQ=255)                                                                                aaagcgtcgacaggaCCGGGTTCGAAC‑GTCGCTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACGCTGTGGGATACAGTTGTACCACCTGTATTGGAAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATTAAAAAAGGCGATTTACCCGTCGGCGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCTCTGGTTAAAACTAActggct                                                                                                                                                                     >  2:516526/15‑251 (MQ=255)                                                                                    cGTCGCTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGTTGGcc                                                                                                                                                                 <  2:566375/251‑1 (MQ=255)                                                                                       cGCTGGCACCGGGTTCGAAAAGTCGTT‑CTGATTAACTTGCAAAAGCGGAACTGGCCCCGGATCCCCGCCGA‑CTTGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcg                                                                                                                                                              <  1:470317/251‑1 (MQ=255)                                                                                         cTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgcc                                                                                                                                                            >  2:53546/1‑251 (MQ=255)                                                                                                      cGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGtt                                                                                                                                               >  1:1961/1‑251 (MQ=255)                                                                                                              cgtttgttaTTATCTGGCAAAAGCTAAACGGACACCGTATCTC‑GACGAACTGGGTTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGTGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATATATGCGTGGGTGCGGTGCTGTCCTGCAACCGCAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAATTAACTTGCTGGCCTCGTCGCCGCTGGTGGTTGCCTAt                                                                                                                                         <  2:568987/242‑1 (MQ=255)                                                                                                              cGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATg                                                                                                                                        <  1:267704/251‑1 (MQ=255)                                                                                                               gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgc                                                                                                                                       <  2:171187/251‑1 (MQ=255)                                                                                                                     gATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCt                                                                                                                                    <  1:312098/248‑1 (MQ=255)                                                                                                                     gATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGGAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTgg                                                                                                                                  <  2:487033/250‑1 (MQ=255)                                                                                                                          tctggcaacagccacccTGACACCGTATATA‑GACAAACCGGGGTCTAGCCTTCTGAGATACGGCTGAACCCCCTGTCTTGGTAACTCTGGGCCGCTGCAAGATCCTAACGAAACGGCAATCAAAAAAGGCTATTTAACCGTCGGAGCGGGGCTGTCCGGAAACCGTCATTTTGAAGGCCGTATCCATCCGATGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCATATGCGCTGGCGGGa                                                                                                                             >  2:67627/17‑250 (MQ=255)                                                                                                                                 aaaGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGTCAATCAAAAAAGGCGATTTAACCGTCGGTTCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGaa                                                                                                                     <  1:527066/251‑1 (MQ=255)                                                                                                                                  aaGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGAATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGCGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAAATAACTGACTGGCCTCGCCGACGCTGGCGGTTGCCTATGCGCTGGCGGGAAATATGAat                                                                                                                    >  1:33694/1‑251 (MQ=255)                                                                                                                                   aGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAat                                                                                                                    >  1:11723/1‑250 (MQ=255)                                                                                                                                    gCGAATCTGACACCGTGTCTC‑GCCGAACTGTGGTTTAACCTTGTGGGATACGTTTGTACCCCCTGTATTGGTAACTCTGTGCCGCTGCCCGATCCTATCGAAACGGCATTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGGAACTTTGAAGGCCGTTGCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCTCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAatat                                                                                                                  <  1:70745/251‑1 (MQ=255)                                                                                                                                    gCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAatat                                                                                                                  >  1:100437/1‑251 (MQ=255)                                                                                                                                         aCTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTAGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAAcc                                                                                                             <  1:310316/251‑1 (MQ=255)                                                                                                                                             acacCGTATCTC‑GACGAACTGGGGTTTACCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGc                                                                                                         <  2:337268/251‑1 (MQ=255)                                                                                                                                                  gTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAACGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCTTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATTACTATCAACCTGGCTTCtt                                                                                                    >  1:537/1‑250 (MQ=255)                                                                                                                                                      ccc‑gacgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGTCGGGAAATATGAATATCAACCTGGCTTCTGAGcc                                                                                                <  1:448985/249‑1 (MQ=255)                                                                                                                                                            cgcaCTGGGGTTTAACCTTGTGGGATACGTTTGTACCACCGGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTTGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTTATCCTATCg                                                                                           <  2:68320/248‑1 (MQ=255)                                                                                                                                                            cgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCg                                                                                           <  1:247649/251‑1 (MQ=255)                                                                                                                                                              aaCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGc                                                                                         >  1:329341/1‑251 (MQ=255)                                                                                                                                                                 tGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCAt                                                                                      >  2:441331/1‑251 (MQ=255)                                                                                                                                                                       ttAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATAGCAACCTGGCTTCTGAGCCTATCGGCCATGATCg                                                                                 <  2:287083/250‑1 (MQ=255)                                                                                                                                                                          actttGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCAAGATCGCaa                                                                              >  2:421842/4‑250 (MQ=255)                                                                                                                                                                                    aTACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAACATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCgg                                                                   <  1:142901/251‑1 (MQ=255)                                                                                                                                                                                         gTTGTACAACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTAt                                                              <  1:284459/251‑1 (MQ=255)                                                                                                                                                                                                   cTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATTGCAAAGGCGATCCGGTTTATCTGAAAGat                                                     >  2:452918/1‑250 (MQ=255)                                                                                                                                                                                                     gTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGatat                                                   <  1:321888/250‑1 (MQ=255)                                                                                                                                                                                                     gTATTGGTAACGCTGGGCCGCTGCCCGTTCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATc                                                  <  2:100437/251‑1 (MQ=255)                                                                                                                                                                                                       aTTGGTAAATCTGGGCAGCTGCCTGATCCTTTCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAGCTTTGAAGGCTGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAAGATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTGTATCTGAAAGATATcg                                                 <  2:331914/250‑2 (MQ=255)                                                                                                                                                                                                             aaCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGGCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGTCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCAt                                           <  2:247613/250‑1 (MQ=255)                                                                                                                                                                                                                   gggCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGc                                       <  1:53546/248‑1 (MQ=255)                                                                                                                                                                                                                              cGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGccc                         <  1:172902/251‑1 (MQ=255)                                                                                                                                                                                                                                    tATCGAAACGGCAACCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCCTCGGCACAAGAAATTGCCCGTGCgg                   >  1:190931/1‑251 (MQ=255)                                                                                                                                                                                                                                     aTCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCgg                   <  2:170743/250‑1 (MQ=255)                                                                                                                                                                                                                                             ggCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAAGTGCCCGTGCGGTAGAACa            >  2:128753/1‑249 (MQ=255)                                                                                                                                                                                                                                               cAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGtc        >  1:174058/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCAACGGCACAAGAAATTGCCCGTGCGGTAGAACAAGtctac     >  1:414542/1‑248 (MQ=255)                                                                                                                                                                                                                                                     aaaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCcaca  <  1:351932/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      TAACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAA‑GTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTC‑GACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCACA  >  NC_000913/1337147‑1337639

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