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%	245.9 / NA	74	P175S (CCA→TCA)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/36);  total (38/36)

GGCGTACCTGCCGTTGTTGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCG  >  NC_000913/1336107‑1336596                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ggCGTACCTGCCGTTGTTGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTAGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTcagg                                                                                                                                                                                                                                                 >  2:161101/1‑251 (MQ=255)            cGTTGTTGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCAttt                                                                                                                                                                                                                                      <  1:339918/251‑1 (MQ=255)               tgttgATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATCCCGCAAAGGTTAACCCGCTATCACCGGTCGAACTGGTCATTGACCAATCGGTGACCGTCGCTCATTTTGGTGATGAAGAGGCATCTGAAGAAAAAGTACGCCTGGAACTGTAGCGCAACCACGAACGTTATGTGATCCTCAAATGGGGACAGCAAGCGTTCAGTCGGTTTAGCGTCGGGCAGTCAAGCACAGGTATTTGcc                                                                                                                                                                                                                                   >  1:252749/1‑251 (MQ=255)                     tCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAgg                                                                                                                                                                                                                             >  2:163075/1‑251 (MQ=255)                     tCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAgg                                                                                                                                                                                                                             >  2:339834/1‑251 (MQ=255)                            gcAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCt                                                                                                                                                                                                                      <  1:531810/251‑1 (MQ=255)                                     gAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATAtct                                                                                                                                                                                                             <  2:456895/251‑1 (MQ=255)                                         cGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGc                                                                                                                                                                                                         <  1:1945/251‑1 (MQ=255)                                          ggTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCa                                                                                                                                                                                                        <  1:398670/251‑1 (MQ=255)                                                   ccTCGGCGGCGATACTGCAAATGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAgtgt                                                                                                                                                                                               <  2:299837/251‑1 (MQ=255)                                                    cTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGtgtg                                                                                                                                                                                              >  1:129890/1‑251 (MQ=255)                                                      cggcggCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGtgtg                                                                                                                                                                                              <  1:30011/249‑1 (MQ=255)                                                      cggcggCGATACTCCAAAGGTTAACCCGCTCTCACCGGTCGGCCTGGTCATTGACCACGCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTGCCTGAAATGGGGAAAGCAAGCGTGCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCGCGGCAAAGCAGTGTGGa                                                                                                                                                                                            <  2:547558/251‑1 (MQ=255)                                                       ggcggcGCTACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGCGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTgg                                                                                                                                                                                             >  2:352464/1‑249 (MQ=255)                                                            cgATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGCTGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTgcat                                                                                                                                                                                      >  2:27232/1‑248 (MQ=255)                                                                 cTGCAAAGGATAACAAGCTCTCATCGGTCGACCTGGTCATTGACCACCCGGTGACCGTCGATCGTTTTGGTGATCATGAGGCATTTGAAGCAACCGTACGCCTGGAAATGGAGCGCAACCACGAATGTTATGTGTTCCTGAAATGCGGAACGCAAGCGTTCAGTCGGTTTCGCGTCGTGCCGTAAGGCACAGGCAGTTGCCATCCGGTTTAACTCGAATATCTCGGACAAGCAGTGGGGAGAGCATTGc                                                                                                                                                                                   >  2:84176/1‑249 (MQ=255)                                                                        ggTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTg                                                                                                                                                                          <  1:299264/251‑1 (MQ=255)                                                                        ggTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTAGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTg                                                                                                                                                                          >  2:480318/1‑251 (MQ=255)                                                                            aaCCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGGACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATg                                                                                                                                                                      >  1:563216/1‑251 (MQ=255)                                                                               ccgctactaCCGGTCGACCTGGTCATTGACCAGTCGGTGACAGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTggatggat                                                                                                                                                                   <  2:478930/243‑1 (MQ=255)                                                                                  ctctCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGc                                                                                                                                                                >  2:216446/1‑251 (MQ=255)                                                                                    ctcACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCtt                                                                                                                                                              >  2:320984/1‑251 (MQ=255)                                                                                     tcACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCtt                                                                                                                                                              >  1:136457/1‑250 (MQ=255)                                                                                         aggTCGACATGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGAGATGATGAGGCATTTGAAGCAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACCGGCATTTGCCATCCGGTTAACCTCGAATATCTCGGCAAAGCAGTGGGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCg                                                                                                                                                         >  2:6759/2‑251 (MQ=255)                                                                                           gTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCa                                                                                                                                                                                                        >  1:151328/1‑202 (MQ=255)                                                                                           gTCGACCTGGTCAGTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCa                                                                                                                                                                                                        <  2:151328/202‑1 (MQ=255)                                                                                             cGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGAt                                                                                                                                                      >  1:405180/1‑250 (MQ=255)                                                                                             cGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGAt                                                                                                                                                      <  1:120110/250‑1 (MQ=255)                                                                                                 cTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGTAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTGAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACt                                                                                                                                                 >  1:542347/1‑251 (MQ=255)                                                                                                    gTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGt                                                                                                                                              <  2:133358/251‑1 (MQ=255)                                                                                                      cATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGtt                                                                                                                                             <  1:399703/250‑1 (MQ=255)                                                                                                              aCTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGAc                                                                                                                                    <  1:327682/251‑1 (MQ=255)                                                                                                               cTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGCTTGCTTATCCGGATACACTCGTTGGTACTGAcc                                                                                                                                   >  2:182949/1‑250 (MQ=255)                                                                                                                       ccGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCAcacc                                                                                                                           <  1:179431/251‑1 (MQ=255)                                                                                                                       ccGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCAcacc                                                                                                                           <  1:226294/251‑1 (MQ=255)                                                                                                                          gcGATCGTTTTGGTGAGGATGAGGCATTGGAAGAAAACGTACGGCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACaccacc                                                                                                                        <  1:341246/250‑1 (MQ=255)                                                                                                                            gATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCAt                                                                                                                      <  1:49018/251‑1 (MQ=255)                                                                                                                                           gatgaGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCACCCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTgg                                                                                                       >  1:512385/1‑251 (MQ=255)                                                                                                                                             tgaGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCGTGATCAACGACCTTGGCg                                                                                                     >  2:464378/1‑251 (MQ=255)                                                                                                                                               aGGCATTGGAAGAAAAAGTACGCCCAGGAATGTAGCGAAAACACGAACGTTATGTGTTCATGAAACGGGAACAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGAATTTGCCATCAGGTTAACCTGGAATATCTCGGCAAAGCAGTGTGGTGTGAATTGCAGGACGGTGAATGGATTGCTCATCAGGATACACTCGTAGGTACTGAATCGCACACCACCATGATAGACGGCCTTGGCGTg                                                                                                   >  2:31440/1‑251 (MQ=255)                                                                                                                                                     ttGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCtgggg                                                                                             >  1:71300/1‑251 (MQ=255)                                                                                                                                                        aagaaAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCtggggtgg                                                                                          >  1:500681/1‑251 (MQ=255)                                                                                                                                                                 tACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTtggt                                                                                 >  1:468405/1‑251 (MQ=255)                                                                                                                                                                  aCGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGACATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTCCTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTtggtg                                                                                >  2:116882/1‑251 (MQ=255)                                                                                                                                                                        ggAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTTGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATc                                                                           <  2:275904/250‑1 (MQ=255)                                                                                                                                                                        ggAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACt                                                                                                                                                 <  1:121546/180‑1 (MQ=255)                                                                                                                                                                        ggAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACt                                                                                                                                                 >  2:121546/1‑180 (MQ=255)                                                                                                                                                                        ggAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCg                                                                          <  1:272545/251‑1 (MQ=255)                                                                                                                                                                          aaaTGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGaa                                                                        <  1:193029/251‑1 (MQ=255)                                                                                                                                                                          aaaTGGAGAGCAACAAAGAACGTTATGTGTTCCTGACATGCGGAAAGCAAGCGTTCAGTCGGATTAGCGTCGTGACGTCAGGCCCAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGCTTGCTTATCCGGATACACTCGTTGGTACTGACTAGCACACCACCGTAATCAACGGCCTTCGCGCGCTGGGGGGGGGCGGTGGTGGGATCgc                                                                         >  2:361229/1‑249 (MQ=255)                                                                                                                                                                            aTGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGAACGAAGc                                                                      >  2:286345/1‑251 (MQ=255)                                                                                                                                                                             tGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGc                                                                      >  1:554248/1‑250 (MQ=255)                                                                                                                                                                              ggAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCa                                                                     <  1:339834/250‑1 (MQ=255)                                                                                                                                                                               gAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGa                                                                   <  1:163075/251‑1 (MQ=255)                                                                                                                                                                                 gcgcAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAg                                                                 <  1:489945/251‑1 (MQ=255)                                                                                                                                                                                        cACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAAt                                                          <  2:468405/251‑1 (MQ=255)                                                                                                                                                                                          cGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGt                                                        <  1:161101/251‑1 (MQ=255)                                                                                                                                                                                              cGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAgg                                                    >  1:317394/1‑251 (MQ=255)                                                                                                                                                                                                ttATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGcc                                                  >  1:406768/1‑251 (MQ=255)                                                                                                                                                                                                    gtgtTCCTGACATGGGGAAAGCAAGCGTCCAGTCGGCTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGcc                                              >  2:103212/1‑251 (MQ=255)                                                                                                                                                                                                    gtgtTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGcc                                              <  1:352464/251‑1 (MQ=255)                                                                                                                                                                                                           tGAAATGGGGCAAGCAAGCGTTCCGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGAGGAGTGCATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTcc                                       >  2:247539/1‑251 (MQ=255)                                                                                                                                                                                                            gAAATGGGGAAAGCAAGCGCTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGCCTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGACGCCGCAATGTTAGGCCAGCCGGTTTcc                                       >  2:365539/1‑250 (MQ=255)                                                                                                                                                                                                                 ggggAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCtt                                 <  1:27232/251‑1 (MQ=255)                                                                                                                                                                                                                      aagcaagcGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATccc                            <  2:129890/251‑1 (MQ=255)                                                                                                                                                                                                                       agcaagcGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCCGCCGGTTTCCATGCTTATCCCg                           >  1:131106/1‑251 (MQ=255)                                                                                                                                                                                                                       agcaagcGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCg                           <  1:320984/251‑1 (MQ=255)                                                                                                                                                                                                                             cGTTAAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTa                     <  2:110218/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGc               >  1:463723/1‑251 (MQ=255)                                                                                                                                                                                                                                           aGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAAc        <  1:75573/250‑1 (MQ=255)                                                                                                                                                                                                                                             cgtcgtGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGCATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCCCCATGATCAACGGCCTTGGCGTGCTGGGGAGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGATAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTCGTGGGCTTCAAACtt      >  1:535677/1‑250 (MQ=255)                                                                                                                                                                                                                                              gtcgtGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTAc    >  1:43581/1‑251 (MQ=255)                                                                                                                                                                                                                                               tcgtGCCGTCAGGCACAGGCCTTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGATGGAATTGCAGGACGTTGAATGGATTGCTTATCAGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCATAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCGTATCCCGGATGTAGTGGGCTTCAAGCTTa     <  2:444980/249‑1 (MQ=255)                                                                                                                                                                                                                                                cgtGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCg  <  2:405180/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    GGCGTACCTGCCGTTGTTGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTACCG  >  NC_000913/1336107‑1336596

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