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	98.2%	182.1 / ‑6.6	57	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (26/30);  minor base T (0/1);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCACAG  >  NC_000913/1337150‑1337640                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       cccAAGTGGGCTGCTGGCCGCAGGCTGGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCGGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAtt                                                                                                                                                                                                                                                  <  2:28804/251‑1 (MQ=255)  ccAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAttt                                                                                                                                                                                                                                                 >  1:172976/1‑251 (MQ=255)         ggCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                                                                                                                                                                                                          <  1:322740/250‑1 (MQ=255)              aTGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgc                                                                                                                                                                                                                                     >  1:378379/1‑251 (MQ=255)                  ccGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGATGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAANCAAAAAAGGCGATTTAACCGTcggtgcggtg                                                                                                                                                                                                                                 >  2:382236/1‑251 (MQ=255)                    gCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAACCTGACACCGTATCTCGACGAACTGGGGGTTAACCTTGTGGGATACGGTTGTACCCCCTGTATTGGTAACTCTGGGCCGCGGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCt                                                                                                                                                                                                                               <  2:745292/251‑1 (MQ=255)                        gCTTGCTGGCGAAAAAAGCCGTCACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGAGGCGCTGGC‑ACCGGGATCGAAAGCCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGGGGGATACGGTTGTACCACCTGCATTGGTAACTCTGGGCCGCTGCCCGATCCTATAGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTcc                                                                                                                                                                                                                           >  2:55808/1‑251 (MQ=255)                           tGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCgg                                                                                                                                                                                                                         <  1:11865/250‑1 (MQ=255)                              gggCGAAAAAAGCCGGAACTCGGGGCTTCAAGGGGCAACCAGGGGTCAAAGCGGCGCTGGC‑ACCGGGTTCGAAAGTCGGTTCTGATTATCTGTCAAAAGCGAAACTGACTCCGTATCTCGACAAACTGGGGTTTAACCTTGTGTGATACGGTTGTACCACCTGTATTGGTAACTCTGTGCCGCGGCCCGATCCTATCGTAACGGCAATCAAAAAAGGCGATTTGACCGTCGGTGCGGTTCGGTCCGTCAAc                                                                                                                                                                                                                     <  2:737748/250‑1 (MQ=255)                                   aaaaaaGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCGGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACTGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTaa                                                                                                                                                                                                                <  2:663122/251‑1 (MQ=255)                                     aaaaGCCGAAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGTACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCACGCAAAAAAGGCGATTTAACCGGCGGTGCGGGGCTGGCCGGCAAACGTAACt                                                                                                                                                                                                              >  2:35455/1‑251 (MQ=255)                                       aagacGTAACTCTTGGCCACAAGCGTCAACCATGGGTCAAAGCGGCGCTGGC‑ACCGGGTTCGAAAGCCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGACCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAGCCGTAACttt                                                                                                                                                                                                            <  2:546176/247‑1 (MQ=255)                                           cGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGa                                                                                                                                                                                                          >  1:693182/1‑249 (MQ=255)                                                ctctGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑AACGGGTTCGAAAGTCGTTTCTGATTAACAGGCAAAAGAGAAACTGACACCGTGTCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCCATCAAAAAAGGCGATTTAACAGCCGGTGCGGTGCTGTCCGGCAACCGTAACTTTgccgaccgt                                                                                                                                                                                                   >  2:453856/1‑243 (MQ=255)                                                          caatagtcAACCATGTGTCAAAGCGTCGCTGGC‑ACATGAGTCGAAAGTCGTTTATGATTATGTGGCAAAAGCGAAATTTACACCGTATCTCGTCGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTTGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGTCGATTTAACCGTCGCTGCGGTGCTGTCCGGAAACCGTAACTTTGAAGGCCGTatccagcc                                                                                                                                                                                           <  2:118412/242‑4 (MQ=255)                                                              gcgtacccATGGGTGAAAGCGTCGCTGTC‑ACCGGGTTCGAAAGTCGTTTCATATTATCTGGCAAAAGCGAAACTGACACCGTAGCTCGACGTACTGGGGTTTAACCGTGTGGGGTACGGTTGTACCACCTGTATTTGTCACTCTGCGCCGCTGCCCGCTCCTATCGTAACGGCAATCAAAAAAGGCGATTGAACCGTCGGTGCGGTGGTGTCCGTCAACCGTAGCTTTGAAGCCCGTATCCATCCTCTGGt                                                                                                                                                                                     <  2:631574/245‑1 (MQ=255)                                                               ggCAACCATGGGATAACGCGTCGCTGGCCACGGGGCT‑GAAAGTCGTTTGTGATTATCTGGCAAAAGCTAAACTGACACCGTATCTCGACGAACTGGGGTTTATCCGTGTGGGATACGGTTGTACCACCTGTATTCGGAACGCTGGGCCGGTGCCCGATCCTATCGAAACTGCAATCAAAAAAGGCGATGTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACGTTGAAGGCCTTATCCATCCGCTGGtg                                                                                                                                                                                    <  2:710699/251‑2 (MQ=255)                                                                 gaaCCAGGGGACAAAGCGGCGCTGGC‑ACCGGGTGCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTAGCTCGACGAACTGGGGTTTAGCCTTGTGGGATACGGTGGTACCACCTGTATTGGTAACTCTTGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTTCGGCAGCCGTAACTTTGAAGGCCGTATCCCTCCGCTGGTTaa                                                                                                                                                                                  <  2:621419/250‑1 (MQ=255)                                                                        gggTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaactaac                                                                                                                                                                           >  1:44380/1‑251 (MQ=255)                                                                          gTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGACGGCCGTATCCATCCGCTGGTTAAAACTAActg                                                                                                                                                                         >  1:106594/1‑251 (MQ=255)                                                                                     gCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgc                                                                                                                                                              <  1:72987/251‑1 (MQ=255)                                                                                               gggTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTAACCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtgg                                                                                                                                                     >  2:561615/1‑251 (MQ=255)                                                                                                ggTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggg                                                                                                                                                    >  1:184739/1‑250 (MQ=255)                                                                                                 gTTCGAAAGTCGTTTCTGATTATCTGGCATAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTACCCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTCACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGc                                                                                                                                                                                                      >  2:499149/1‑200 (MQ=255)                                                                                                 gTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCGCCTGTATGGTTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGc                                                                                                                                                                                                      <  1:499149/200‑1 (MQ=255)                                                                                                 gTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtg                                                                                                                                                   >  1:478754/1‑251 (MQ=255)                                                                                                     gAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTg                                                                                                                                               <  1:316823/251‑1 (MQ=255)                                                                                                                          ggCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGaaaa                                                                                                                           >  2:230523/1‑249 (MQ=255)                                                                                                                          ggCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGaa                                                                                                                             <  1:403812/248‑1 (MQ=255)                                                                                                                             aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGGAATATg                                                                                                                        >  1:678520/1‑250 (MQ=255)                                                                                                                             aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCGCGCCGCCGCTNGTGGTTGCCTATGCGCTGGCGGGAAATATg                                                                                                                        >  2:375146/1‑250 (MQ=255)                                                                                                                                aGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTGGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCTGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAata                                                                                                                    <  1:55808/251‑1 (MQ=255)                                                                                                                                   gAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATc                                                                                                                  >  1:298236/1‑250 (MQ=255)                                                                                                                                              cGTATCTCGAAGAACTGGGGTTTAACCTTGGGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCTCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCt                                                                                                      <  1:35455/251‑1 (MQ=255)                                                                                                                                                 aTCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTCCCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGa                                                                                                    <  2:89058/250‑1 (MQ=255)                                                                                                                                                    tcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATTACTATCAACCCTGCTTCTTAGCCt                                                                                                >  2:7925/1‑251 (MQ=255)                                                                                                                                                    tcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCt                                                                                                >  1:634427/1‑251 (MQ=255)                                                                                                                                                       acgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCGGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATc                                                                                             <  2:30853/251‑1 (MQ=255)                                                                                                                                                             tGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCAt                                                                                       >  1:763237/1‑251 (MQ=255)                                                                                                                                                              ggggTTTATCCTTGTGGGGTACGGTTGTACCACCTGTATTGGTAACTCGGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCAt                                                                                       <  2:678520/250‑1 (MQ=255)                                                                                                                                                                     aaCCTTGTGGGATACGGTTGTACCACCGGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGAAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCACCCTGGCTTCTTAGCCTATCGGCCATGCTCGCa                                                                                >  2:654685/1‑250 (MQ=255)                                                                                                                                                                          tgtgGGATGCGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAgg                                                                            <  2:763237/249‑1 (MQ=255)                                                                                                                                                                          taggggATACTGTTGTACCACCGGTATTGGNAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATACCAACCTGGCTTCTGAGCCTATCGGCCCTGATCGCAAAGGc                                                                           <  2:378379/247‑1 (MQ=255)                                                                                                                                                                            ggggtagctggtGTACCGCCTGGATTGGTATCTCTGGGACGTTGCCCGAACCTACTGAAACGGAAAGCAAAAAATGCGATTTAACAGTCGGTGCGGTGCGTTGCGGCAACCGAACCTGTGAAGGCCGGATCCAGCCGTTGGTTAAAACTATCTGGCTGGCCTCGCCGCCGCTGGTGGTGGCAGAGGCGCTGGCGGGAAAAATGAATAGCAACCTGGCTTCTGAGCCTATCGGCCATGAtttgaaagtcgat                                                                        <  2:282941/240‑13 (MQ=255)                                                                                                                                                                             gggATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTTCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGGATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATc                                                                       <  1:734621/251‑1 (MQ=255)                                                                                                                                                                                  aCGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGtt                                                                  >  1:614750/1‑251 (MQ=255)                                                                                                                                                                                                 gTATTTGTAAATCTGGGCCGCTCCCCGATCCTATCGAAACGGCAATTAAAAAAGGCGATTGAACCGTCGGTTCGGTGCGGTCTGTGAGCCGTAACGTTGAATGCCGTATCCATCCGCTGGTTAAATCTTACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCTTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTTTGCGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGatat                                                    <  1:450722/250‑1 (MQ=255)                                                                                                                                                                                                   aTTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCtt                                                                                                        <  1:339098/196‑1 (MQ=255)                                                                                                                                                                                                   aTTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCCGGCtt                                                                                                        >  2:339098/1‑196 (MQ=255)                                                                                                                                                                                                     tatgaaccccgggCCGCTGCCCGATCCTTTCGAAACGGCAATAAAAAAAGGCGATTTAACCGTCGGTGCGGTGTTGTCCGGCATCCTTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTgg                                                <  1:339004/240‑1 (MQ=255)                                                                                                                                                                                                      ggTAACTCGGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGcc                                              <  1:62359/251‑1 (MQ=255)                                                                                                                                                                                                       gTAACTCTGGGCCGCTTCCCGACGCTATCCAAACGGCAACCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCGATGCGCTGGCGGGAAATATGAATATCATCCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCa                                             <  2:478754/251‑1 (MQ=255)                                                                                                                                                                                                                       gCCCGATCCTATCGAACCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAAt                               <  2:634427/249‑1 (MQ=255)                                                                                                                                                                                                                              ccTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCATCTGAACCCATCGGCCATTAAAGCAAAGGCGATCCGGCTTATCTGAAAGCAATCTGGCCATCGGCCCAAGGAATTGCCCGTgg                      >  2:736995/1‑250 (MQ=255)                                                                                                                                                                                                                                         ggCAATCAAAAAAGGCGATTTACCCGTCGGTGCGGTGATGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACa             >  1:581379/1‑249 (MQ=255)                                                                                                                                                                                                                                               cAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGGGATCCGGTTTATCTGAAAGATATTTGGCCATCGGCACANAAAATTGCCCGTGCGGGAGAACAAGtccccc     >  2:373837/1‑247 (MQ=255)                                                                                                                                                                                                                                               cAAAAAAGGCGATTTAACCGGCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTcc      <  1:382236/250‑1 (MQ=255)                                                                                                                                                                                                                                                aaaaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCcac    <  1:637238/251‑1 (MQ=255)                                                                                                                                                                                                                                                  aaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCACAg  <  1:230523/251‑1 (MQ=255)                                                                                                                                                                                                                                                  aaaagGCGAGTTAACGGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCACAg  <  2:44380/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       CCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGC‑ACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCACAG  >  NC_000913/1337150‑1337640

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