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	96.1%	164.6 / NA	51	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (19/30);  minor base T (0/2);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAG  >  NC_000913/1337148‑1337625                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        aaCCCAAGGGTGCTTATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTCTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCg                                                                                                                                                                                                                                      <  2:275902/250‑1 (MQ=255)   cccAAGTGTGCTGATGGCCTCTGGCTTGCTGGGGAAAAAAGCCGTAACTCTGGGCCTCAATCGTCAACCAGGGGTCAAAGCGT‑CGCTGGCACCGGGTGCGAAAGTCGGTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAtt                                                                                                                                                                                                                                   <  2:290089/251‑1 (MQ=255)     cAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTa                                                                                                                                                                                                                                 <  1:550515/251‑1 (MQ=255)          gtgCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGt                                                                                                                                                                                                                            >  2:133346/1‑251 (MQ=255)                    ccGCAGGCTTTCGGGCGAACAAAGCCGTAACTCTGGGCCGCAAGCGGCAGCCATGGGTCAAAGCTG‑CGCTGGCAACGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCTATTTAACCGTcggtgcggtg                                                                                                                                                                                                                  <  2:359473/251‑1 (MQ=255)                           cTTGCTGGCGAAAAAAGCCGTAAATCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCCATCATAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCg                                                                                                                                                                                                           >  2:438844/1‑251 (MQ=255)                                                    ctGGGCCTCAAGCGGCAACCATGGGTCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGt                                                                                                                                                                                    <  1:118147/249‑1 (MQ=255)                                                      gggCCTCATGCGGCAAACATGGGTCAATG‑GTACGCTGGCGCCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAAGCAAAAAAGGCGATTTATCCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatcc                                                                                                                                                                                <  2:396299/251‑1 (MQ=255)                                                         ccTCAAGCGGCAACCATGGGTCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAATGCCGTatccatc                                                                                                                                                                             <  1:412329/251‑1 (MQ=255)                                                              aGCGGCAACCATGGGTCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTg                                                                                                                                                                        <  1:10632/251‑1 (MQ=255)                                                                 ggCATCCATTGGTCGAAGCTG‑CGCTGGCCCCGGGTTCGAAAGTCGTTTCGGTTTATCTGGCAAACGCGAAACTGACGCCATATCTCGACGAACTGGGGTTTAACCTTGGGGGATCCGGTTGTACCACCTGTATTGGTAACTCGGGGCCGCTGCCCGGTACTATCGAAACTGCAATCAAAAAAGGCGATTTAACCTTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTCTCCATCCGCTGGtt                                                                                                                                                                     <  2:3154/251‑1 (MQ=255)                                                                    aaCCATGGTTCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCGCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTaaa                                                                                                                                                                  <  1:550357/251‑1 (MQ=255)                                                                             tCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGAATCGGTTGTACCACCTGTATTGGTAACTCTGAGCCGCTGCCCGATCCTATCGAAACGGCGATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActgg                                                                                                                                                         >  2:11445/1‑251 (MQ=255)                                                                                  gCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACtggctggc                                                                                                                                                    <  1:239279/251‑1 (MQ=255)                                                                                        gCTGGCACCGGGTTCGAAAGTCGTTTCTGATGATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGTTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgc                                                                                                                                               <  1:481802/251‑1 (MQ=255)                                                                                           ggCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgcc                                                                                                                                              <  1:185667/249‑1 (MQ=255)                                                                                                 gggTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtgg                                                                                                                                      >  1:110471/1‑251 (MQ=255)                                                                                                  ggTTCGAAAGTCGTTTCTGATTCTCCGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGCTTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatacat                                                                                                                                                                              >  2:232404/1‑210 (MQ=255)                                                                                                  ggTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCGATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccat                                                                                                                                                                              <  1:232404/210‑1 (MQ=255)                                                                                                   gTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtg                                                                                                                                    >  1:65454/1‑251 (MQ=255)                                                                                                      cGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGtt                                                                                                                                 >  1:522575/1‑251 (MQ=255)                                                                                                           gTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTa                                                                                                                            <  1:470850/251‑1 (MQ=255)                                                                                                             cGTTTCTGATTATCTGGCAAAAGAGAAACTGACACCTTATCTAGACGAATGGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTTCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAATCTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATg                                                                                                                          <  2:331699/251‑1 (MQ=255)                                                                                                                   tGATTATCGGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAAACTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACAGTCGGTGCGGTGATGTCCGGCAACCGTAACTTTGAAGGCCGTATCAATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCTACGCCGGAGGTTGCCTATGCGCt                                                                                                                      >  2:63920/1‑249 (MQ=255)                                                                                                                    gattagcgggCAAAAGCGAAACTGACACCGTAACACGACGAACTGGGGTTTAACCTTGGGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACGGGCGGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGc                                                                                                                   <  2:169234/243‑1 (MQ=255)                                                                                                                    gATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCt                                                                                                                      <  1:221894/248‑1 (MQ=255)                                                                                                                    gATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCt                                                                                                                      <  1:498670/248‑1 (MQ=255)                                                                                                                    gATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTg                                                                                                                     <  1:133346/249‑1 (MQ=255)                                                                                                                        aTCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCggg                                                                                                                <  1:245263/250‑1 (MQ=255)                                                                                                                            ggCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGa                                                                                                               <  1:452330/247‑1 (MQ=255)                                                                                                                               aaaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGACATATg                                                                                                         >  1:154787/1‑250 (MQ=255)                                                                                                                               aaaaGCCAAACGGACCCCGTATCTCGACGAACTGGGGTTTAACCTTTTGGGATACGGTTGTACCACCTTTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAACCGGCAATCAAAAAATGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGATATATGa                                                                                                        <  2:12790/251‑1 (MQ=255)                                                                                                                                     gAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTA‑CTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCa                                                                                                  >  1:434603/1‑250 (MQ=255)                                                                                                                                       aaCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAAc                                                                                                >  1:320773/1‑251 (MQ=255)                                                                                                                                           gACACCGTATCTCGACGAACTGTGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGCGACGCTGCCCGCTCCTATCGACACGGCAATCAAAAAAGGCGCTTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGACGTATCCATCCGCTGGTTAAAACTAAATGGCTGGACTCGCCGCCGCTGGTGGTTGCCTATGATCTGGCGGAAAATATGAATATCAACCTgg                                                                                            >  1:517870/1‑251 (MQ=255)                                                                                                                                           gACACCGTATCTCGACGAACTGGGGGTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGAGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTgg                                                                                            >  1:488091/1‑251 (MQ=255)                                                                                                                                                                 ggttgtAACCTTGTGGGATACGGTTTTACCACCTTTATGGGGACCTCTGGGCCGCTGCCCGATCCAATCGATACGGCAATCAAAAAAGGCGATTTAACCTTCGGCGCGGTGCTGACCGGCAACCGTAACTTTGAAGGTCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGa                                                                      <  2:20872/246‑1 (MQ=255)                                                                                                                                                                  ggTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGAt                                                                     <  1:209796/251‑1 (MQ=255)                                                                                                                                                                       aaCCTTGTGGGATACGGTTGTACCACCTGTATTGGTATCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCa                                                                 >  1:472959/1‑250 (MQ=255)                                                                                                                                                                             gtgGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCg                                                           <  1:328598/250‑1 (MQ=255)                                                                                                                                                                              ggggATACTGTTTTACCACCTGTTGTGGTAACTCTGTGCCGCGGCCCGATCCGATCGAAACGGTAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGTCCGTATCCATCCGCTTGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGAt                                                         <  2:306600/250‑1 (MQ=255)                                                                                                                                                                                    aCGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGtt                                                   >  1:142875/1‑251 (MQ=255)                                                                                                                                                                                      ggtggtaccccttggaTTGGAAACTTTGGGCCGGTGCCCAATCCTATCGAAACGGCAATCAAAAAAGGAGAGTTAACCGCCGGGGCGGTGCTGTCCGGCAACAGTACCTTTGAAGGCCGTAGCCATCCGCGTGTTAAAACTAACTGGATGGCCTCGCCGCCGCTGGTGGGTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGCTCGCAAAGGCGATCCGGtgta                                                 <  2:218134/236‑4 (MQ=255)                                                                                                                                                                                               cccttgaTTGGTAACCCTGGGCCGTGGCCGGACCCTAACGAAACGGCAATCAAAAAATGCGATTTACCCGTCGGTGCGGTGCTGTCCGGCAACAGTATCTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATAGCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCAGGTTTATCTGaatga                                        <  2:156784/245‑4 (MQ=255)                                                                                                                                                                                                   gTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATc                                    <  1:494736/251‑1 (MQ=255)                                                                                                                                                                                                                            cGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGccc           >  1:192565/1‑251 (MQ=255)                                                                                                                                                                                                                                  tATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCgg     >  1:353929/1‑251 (MQ=255)                                                                                                                                                                                                                                     cgacacggtaagcAAAAAAGGTGATTTAACCGTTGGTGCGGGTCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAGCTAACGTGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCGGGCTTCTGAGCCTATCGGCCATGAGCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGTCCAGCGGGACAAGAAATTGCCCGTGCGGTAg  <  2:566587/239‑1 (MQ=255)                                                                                                                                                                                                                                     cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATc                                                        <  2:147841/197‑1 (MQ=255)                                                                                                                                                                                                                                     cGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATc                                                        >  1:147841/1‑197 (MQ=255)                                                                                                                                                                                                                                      gAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGt    >  1:154985/1‑248 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        AACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGT‑CGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAG  >  NC_000913/1337148‑1337625

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