Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,336,353	C→T	P175S (CCA→TCA)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,336,353	0	C	T	100.0%	161.4 / NA	49	P175S (CCA→TCA)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/29);  total (20/29)

CGGCGTACCTGCCGTTGTTGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGG  >  NC_000913/1336106‑1336597                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     cGGCGTACCTGCCGTTGTTGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTcag                                                                                                                                                                                                                                                   <  2:87311/251‑1 (MQ=255)     gaaCCTGCCGTTGTTGATCTGGCGGAAATGCGCGCAGCGGTTAAACGCATCGGCGGCGATACTGCAAAGGGTAACCCGCTCTCACAGGTCGACCTGGTCCTTG‑CCAACTCGGTGACCGACGATCGTTTTGGTGATGGTGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGTAAGCAAGAGTTAAGTCGGTTTAGCGTCGTGCCGTcaggcac                                                                                                                                                                                                                                               >  2:33799/3‑251 (MQ=255)           gCCGTTGTTGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTcaggcacaggc                                                                                                                                                                                                                                           <  1:95422/249‑1 (MQ=255)              gttgttGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTg                                                                                                                                                                                                                                      >  1:177763/1‑251 (MQ=255)                              cAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACTAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTc                                                                                                                                                                                                                      >  1:35575/1‑251 (MQ=255)                                               aaaCGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCGGGTCATTG‑ACCACTCGGTGACCGTCGAGCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCATCCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCAGCAGGTTAACCTCGAATATCTCGGCAAAg                                                                                                                                                                                                      <  2:338887/250‑1 (MQ=255)                                                   gCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAgtg                                                                                                                                                                                                 <  2:131847/251‑1 (MQ=255)                                                                        aGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGt                                                                                                                                                                            <  1:11018/251‑1 (MQ=255)                                                                        aGGTTAACCCGCTCTCACCGGTCGACCTGGTAATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGt                                                                                                                                                                            <  1:265002/251‑1 (MQ=255)                                                                             aaCCCGCTCTCACCGGTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAAt                                                                                                                                                                        <  1:113995/250‑1 (MQ=255)                                                                             aaCCCGCTCTATCCGGTCCACCTGCTAATTGCTCCAC‑CGGTGACCGATGATAGTTATGGTGATGATACCGCATTCCAAGAAAACGTCCGCCTAGAAATCGAGCGAAACCACGAACGTTATGTGCTACTGAAACGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGAATTTGCCATCAGGTTAACCTCGAATATCTCGGAAAAGCCGTGTGGAGTGAATTGCAGGACGGTGAATg                                                                                                                                                                       >  2:366335/1‑251 (MQ=255)                                                                                 cGCTCTCACCGGTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGAtt                                                                                                                                                                   <  2:264765/251‑1 (MQ=255)                                                                                      tcACCGGTCGACCTGGTCCTTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGCTGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAACGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCtt                                                                                                                                                               >  2:413545/1‑250 (MQ=255)                                                                                            gTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGACTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCgg                                                                                                                                                         >  1:544951/1‑250 (MQ=255)                                                                                                    ggTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTc                                                                                                                                                 <  1:549992/250‑1 (MQ=255)                                                                                                    ggTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCg                                                                                                                                                <  1:216631/251‑1 (MQ=255)                                                                                                           g‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTa                                                                                                                                          <  1:27685/250‑1 (MQ=255)                                                                                                                   cGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTc                                                                                                                                   <  2:342571/250‑1 (MQ=255)                                                                                                                   cGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTc                                                                                                                                   <  1:245374/250‑1 (MQ=255)                                                                                                                              gATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCAt                                                                                                                       >  2:467688/1‑251 (MQ=255)                                                                                                                                 cGTTTTGGTGATGATGAGGCATTTGAAGAAAACTTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTTTTCCTGAAATGGTGAAATCAAGCGTTCAGTCGGTTTAGCGTCGTTCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAATCAGTGTGTAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGCGGGTACCGACTCGCACACCACCATGAt                                                                                                                    <  2:449920/251‑1 (MQ=255)                                                                                                                                     ttGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAACTGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAAc                                                                                                                >  1:236685/1‑251 (MQ=255)                                                                                                                                      tGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAAc                                                                                                                <  1:413545/250‑1 (MQ=255)                                                                                                                                      tGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAAGGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATTGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCAATTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTTTGGAGTGAATTGCAGGACGTTGAATGGATTGCTGATCCGGATACTCTCGTTGGTACTGACTCGCACACCACCAAGATCAAc                                                                                                                <  2:38612/250‑1 (MQ=255)                                                                                                                                       ggTGATGATGAGGAATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACgg                                                                                                              <  2:229584/251‑1 (MQ=255)                                                                                                                                          gatgatgaGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGcc                                                                                                            <  1:462213/250‑1 (MQ=255)                                                                                                                                             gatgaGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTg                                                                                                         <  1:326124/250‑1 (MQ=255)                                                                                                                                                 agtcATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAATCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGTCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCCCCATGATCATCGGCCTTGGCGTg                                                                                                    <  2:422322/248‑1 (MQ=255)                                                                                                                                                 aGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGTCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCGTATCCGGATACACTCGTTGGTACTGACTCGCACACCTCTATGATCAACGGCCGTGGCGTg                                                                                                    <  2:578790/251‑1 (MQ=255)                                                                                                                                                  ggCATTTGAAGAAAACGTACGTCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGCCCTTGGCGTGc                                                                                                   >  1:291424/1‑251 (MQ=255)                                                                                                                                                  ggCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGc                                                                                                   <  1:353947/251‑1 (MQ=255)                                                                                                                                                     agttGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCtgg                                                                                                <  1:251307/249‑1 (MQ=255)                                                                                                                                                           agaaAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGa                                                                                                                                      >  1:372906/1‑207 (MQ=255)                                                                                                                                                           agaaAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGa                                                                                                                                      <  2:372906/207‑1 (MQ=255)                                                                                                                                                                    aCGCCTGGAAATGGCGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGCAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATCTCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGAAGGTGACTGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTtggtg                                                                                 >  1:66719/1‑251 (MQ=255)                                                                                                                                                                          ggAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTGTAGCGTCGTGCCGTCAGGCACAGGCAGTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCGGGGGTGGGGCGTTGGTGGGATc                                                                            <  2:258111/250‑1 (MQ=255)                                                                                                                                                                           gAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCg                                                                           >  1:456392/1‑250 (MQ=255)                                                                                                                                                                             aaTGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGa                                                                          >  1:371971/1‑249 (MQ=255)                                                                                                                                                                                                      gtgtTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGcc                                               <  2:287141/251‑1 (MQ=255)                                                                                                                                                                                                           ccTGAAATGGGGAAAGCCAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAAACTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGAACACCACCATGATCAACGGACTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGttt                                          >  2:588989/1‑251 (MQ=255)                                                                                                                                                                                                              gAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCa                                       <  1:273439/251‑1 (MQ=255)                                                                                                                                                                                                                 aTGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGc                                    >  1:234320/1‑251 (MQ=255)                                                                                                                                                                                                                   ggggAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCtt                                  >  1:170593/1‑251 (MQ=255)                                                                                                                                                                                                                    gggAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACTCTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTa                                 <  1:310055/251‑1 (MQ=255)                                                                                                                                                                                                                                     gTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGAAGTGGGc                >  2:56235/1‑251 (MQ=255)                                                                                                                                                                                                                                        ggTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGTGCTTc             <  2:35575/251‑1 (MQ=255)                                                                                                                                                                                                                                        ggTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTc             >  1:418786/1‑251 (MQ=255)                                                                                                                                                                                                                                               cgtcgtGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACtt       >  1:241775/1‑250 (MQ=255)                                                                                                                                                                                                                                                   gtGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCgg  >  1:542900/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     CGGCGTACCTGCCGTTGTTGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTG‑ACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCGG  >  NC_000913/1336106‑1336597

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