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

CCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTC  >  NC_000913/1337151‑1337635                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  ccAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAttt                                                                                                                                                                                                                                            <  1:364210/251‑1 (MQ=255)    aGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCCCCAGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGGAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAAc                                                                                                                                                                                                                                         >  2:180538/1‑251 (MQ=255)            gaatgCCGCTGGCTTGCTGACGAAAAAAGCCGTAACTCTGGGCCTCAAGTGGCGACCATGGTTCAAAGCGTCGCTGGCACCGCGTTTGTAAGTCGTTTCTGTTTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGTGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGGTCCTATCGATACGGCAATCAAAAAAGGCGATTTAACCATcgg                                                                                                                                                                                                                                   <  2:21703/245‑1 (MQ=255)                               gCGAAAAAGGCCGTAACACTGGTCCTCAAGCGGCAACCATGGGTCAAAGCGTCTCTGGCACCGGGTTCGAAAGTCATTTCTGATTATCTGGCAAAAGCGAAACTGACACCGGATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTAGTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCTGTGCTGTCCGGCAACCg                                                                                                                                                                                                              <  2:111062/251‑1 (MQ=255)                                 aaaaaaaGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTa                                                                                                                                                                                                            <  2:169023/250‑1 (MQ=255)                                                 ctGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTat                                                                                                                                                                                            >  2:160142/1‑251 (MQ=255)                                                   gggCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatc                                                                                                                                                                                           <  1:269020/250‑1 (MQ=255)                                                     gCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccat                                                                                                                                                                                        <  2:147429/251‑1 (MQ=255)                                                     gCCTAAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccat                                                                                                                                                                                        <  2:147424/251‑1 (MQ=255)                                                        tCAAGCGGCAACCATGGGTCACAGCGTCGCTGGCACCGGGTTCGACAGACGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAAGTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTCAGTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAGCCGTCGGTGCGGTGCTGGCCGGCAACCGTAACTTTGAAGGCCGTatcctgccg                                                                                                                                                                                     >  2:150953/1‑246 (MQ=255)                                                                 aaCCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGCTCCTATCGCAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCGGGTTaa                                                                                                                                                                             >  2:9530/1‑250 (MQ=255)                                                                 aaCCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAgggaacggcaaagatatcgatacc                                                                                                                                                                            >  1:150016/1‑229 (MQ=255)                                                                    cATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaact                                                                                                                                                                         <  1:129764/251‑1 (MQ=255)                                                                                                ttCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtgg                                                                                                                                             >  1:449663/1‑251 (MQ=255)                                                                                                  cGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTggtggt                                                                                                                                            <  1:374176/250‑1 (MQ=255)                                                                                                    aaaGTCGTTTCTGATTATCTGGCAATAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGc                                                                                                                                         >  1:441550/1‑251 (MQ=255)                                                                                                      aGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCGATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCt                                                                                                                                       <  1:186539/251‑1 (MQ=255)                                                                                                       gTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTAGACGCACTGGGGTTTACCCCTGTGGGATACGGTTGTACCACCTGTATTGGAAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCtc                                                                                                                                      >  1:40554/1‑250 (MQ=255)                                                                                                         cGTTTCTGATTATCTGTCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGCCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCGTCCGTTTGGTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTAt                                                                                                                                     <  1:465617/250‑1 (MQ=255)                                                                                                               tGATTATCTGGAAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTgg                                                                                                                              >  2:197965/1‑251 (MQ=255)                                                                                                                   taacTGGCAAAAGCGAACCGGACACCGTTGCTCTTCGAACGTGGGTTGACCCTTGTGAGATACGGTTGTACCAGCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCAATCGAAACGGCAATCAAAAAAGTCGATTTAACCGTCGGGGCGGTGCTTTCCGGCAACCGTAACTTTGAAGGCCGTTTCCATCCTCTGGTTAAACCTAACTGGCTTGCCTCCCCGCCGCTGGTGGTTTCCTATGCCCTGGCgg                                                                                                                           <  2:286735/247‑1 (MQ=255)                                                                                                                        ggCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGaaa                                                                                                                       <  2:118598/249‑1 (MQ=255)                                                                                                                           aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCCGGCGGGAAATATg                                                                                                                   >  1:103233/1‑250 (MQ=255)                                                                                                                              aGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAat                                                                                                                >  1:89910/1‑250 (MQ=255)                                                                                                                                    aCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAAcc                                                                                                         <  1:234883/251‑1 (MQ=255)                                                                                                                                     cTGACACCGTATATCGACGAACTGGGGTTTCACCTTGTGGGCTACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAACCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAGCCGTAACTTTGAAGGCCGTATCCATCAGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGGAATATGAATATCAAccc                                                                                                        >  2:217374/1‑250 (MQ=255)                                                                                                                                       gACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCt                                                                                                        <  1:57019/249‑1 (MQ=255)                                                                                                                                         cacCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCt                                                                                                    >  1:131272/1‑251 (MQ=255)                                                                                                                                           ccGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCTATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAGCCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGTCCTCGCCGCCGCGGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTc                                                                                                  <  1:217374/251‑1 (MQ=255)                                                                                                                                               aTCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTAt                                                                                                                                     >  1:117704/1‑212 (MQ=255)                                                                                                                                               aTCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTAt                                                                                                                                     <  2:117704/212‑1 (MQ=255)                                                                                                                                                   cgacgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAg                                                                                              <  2:131272/247‑1 (MQ=255)                                                                                                                                                   cgacgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGACGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAAGATCAACCTGGCTTCTGAGCCt                                                                                           <  1:86761/250‑1 (MQ=255)                                                                                                                                                    gacgcACTGGGGCTTAACCTTGTGGGCTACGGTTGTACCACCTGTCTTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAACAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTCTCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCAGGCGGGGAAAATGGATTTCTACCTGGCTTCTGAGCCTAt                                                                                         >  2:336455/1‑251 (MQ=255)                                                                                                                                                         aCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGcc                                                                                    >  1:346224/1‑251 (MQ=255)                                                                                                                                                                 ttAACCTTGTGGGATACGGTTGTACCACCTGTATGGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCg                                                                             <  1:178736/250‑1 (MQ=255)                                                                                                                                                                   aaCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCaa                                                                          <  1:177959/251‑1 (MQ=255)                                                                                                                                                                    aCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGACGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTAACGGCCATGAtagcca                                                                          >  2:261319/1‑245 (MQ=255)                                                                                                                                                                                 cGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCCGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGttt                                                            <  2:103233/251‑1 (MQ=255)                                                                                                                                                                                                 aTTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCt                                             <  1:180538/250‑1 (MQ=255)                                                                                                                                                                                                     gTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCa                                        <  1:160142/251‑1 (MQ=255)                                                                                                                                                                                                           ctGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAACCTATCGGCCATGATCGCAAAGGCGAACCGGGTTATCTGAAAGATATCTGGCCATccgcc                                  >  2:396357/1‑247 (MQ=255)                                                                                                                                                                                                                 cGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCACTGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATTTGGCCATCGGCACAAGa                             >  1:416871/1‑250 (MQ=255)                                                                                                                                                                                                                 cGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGACGATTTACCCCTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCCCAAAGGCGATCCGGTTTCTCTGAAAGATTACTGGCCATCGGCACAAGa                             >  1:38317/1‑250 (MQ=255)                                                                                                                                                                                                                  gCTGCCCGCCCCTATCGAAACGGAAATCAAAAACGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCACCCGTAAATTTGAAGGCCATCTCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAAGCTGACAATCCACCTGGCTTTTGAGCCAAACGGACCTGATGGCAAAGGCAGTCCAGTTTTTCTGTAAGAAATCTGGCCATCgcacgcagga                            >  2:266730/1‑241 (MQ=255)                                                                                                                                                                                                                                 cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATCTCAACCTGGCTTCTCAGCCTAACGGCCATGATCGCAAAGGCGAACCGGTTTATCTGAAAGATATCTGGCCCTCGGGACAAGAAATTGCCCGTGCGGtcg            >  1:340915/1‑249 (MQ=255)                                                                                                                                                                                                                                 cGAAACGGCAATCAAAAAAGGCGATGTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCAATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGGTGATCTGAAAGATATCTGGCTATCGGCACAAGAAATTGCCCGTGCGGTa             <  2:235059/250‑1 (MQ=255)                                                                                                                                                                                                                                    aaCGGCAATCAAAAAAGGTGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGGAACTTGGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTTGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGAACAAGAAATTGCCCGTGCGTTAGAAc         <  2:299531/251‑1 (MQ=255)                                                                                                                                                                                                                                      cGGCAATCAAAAAAGGCTATTTAATCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCTTCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAAGATCAACCTGGCTTCTGAGCCTATCGGCCTTGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTTGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACaa       <  1:266730/251‑1 (MQ=255)                                                                                                                                                                                                                                      cGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCAGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCACTATCGCAAAGGCGATCCGGTTTATCTGGAAGCTATCTGGCCATCGGCCAAAGGAAATGCCCGTgtcggagaaca        >  2:270950/1‑240 (MQ=255)                                                                                                                                                                                                                                          aaTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGGCCGGCAACCGTAGCTTTGAAGGGCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCGCGCCGCCGCTGGTGGTTGCCTATGCGCGGGCGGGAAATAGGAATATCACCCTGGCTTCTGGGCCTAGTGGCCATGAGCGCAAAGGCGAGCCTGTTTATCTGAAAGATATCTTGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGtct   <  1:214023/251‑1 (MQ=255)                                                                                                                                                                                                                                           aTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGtctc  >  1:173372/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATAAACCTGGGTTCTGAGCCTAACGGCCATGGTCGCAAAGGCGGACCGGTTTATCTGAAAGAAATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGGAGAACAAGtctc  >  2:120574/1‑251 (MQ=255)                                                                                                                                                                                                                                            tCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTg                                            <  1:426100/208‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  CCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTC  >  NC_000913/1337151‑1337635

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: ‑