Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,304,469	T→C	intergenic (+76/+639)	eco → / ← mqo	serine protease inhibitor ecotin/malate:quinone oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,304,469	0	T	C	100.0%	121.9 / NA	35	intergenic (+76/+639)	eco/mqo	serine protease inhibitor ecotin/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/23);  total (12/23)

TGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCC  >  NC_000913/2304336‑2304586                                                                                                                                      |                                                                                                                      tgtgAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTg                                                                                                                 >  2:966770/1‑140 (MQ=255)  gtgAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGa                                                                                                                >  2:678874/1‑140 (MQ=255)  gtgAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGa                                                                                                                >  2:534914/1‑140 (MQ=255)   tgAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGAt                                                                                                               >  1:232558/1‑140 (MQ=255)     aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGc                                                                                                             >  2:448015/1‑140 (MQ=255)         aCCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACg                                                                                                         <  1:330264/140‑1 (MQ=255)         aCCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACg                                                                                                         >  1:929010/1‑140 (MQ=255)             cgTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTgg                                                                                                     <  1:158490/140‑1 (MQ=255)             cgTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTgg                                                                                                     <  2:458873/140‑1 (MQ=255)             cgTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTgg                                                                                                     <  2:631956/140‑1 (MQ=255)             cgTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTgg                                                                                                     >  2:312735/1‑140 (MQ=255)                  ggAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCGGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGt                                                                                                <  1:570211/140‑1 (MQ=255)                    aaGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTgg                                                                                                     >  1:615621/1‑133 (MQ=255)                    aaGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTgg                                                                                                     <  2:615621/133‑1 (MQ=255)                       gCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAACCTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                           <  2:857139/140‑1 (MQ=255)                       gCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                           <  2:141646/140‑1 (MQ=255)                       gCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                           <  1:678874/140‑1 (MQ=255)                       gCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                           <  2:953206/140‑1 (MQ=255)                       gCGGAAGAAAAAATTCAAAACGCGGTAGTTGGCTAAACTGCCGTGAAGGGCGGCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCGAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                           <  2:952384/140‑1 (MQ=17)                              agaAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCt                                                                                    <  2:190999/140‑1 (MQ=255)                                      gACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCt                                                                            <  1:312735/140‑1 (MQ=255)                                              ggTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTg                                                                    >  1:894433/1‑140 (MQ=255)                                              ggTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTg                                                                    <  2:929010/140‑1 (MQ=255)                                                tAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCg                                                                  <  1:534914/140‑1 (MQ=255)                                                  gTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGc                                                                <  1:966770/140‑1 (MQ=255)                                                   ttCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCa                                                               <  2:276714/140‑1 (MQ=255)                                                   ttCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCa                                                               <  2:962061/140‑1 (MQ=255)                                                       cTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCAcccc                                                           >  1:867972/1‑140 (MQ=255)                                                        tAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCg                                                          <  2:441445/140‑1 (MQ=255)                                                           aCTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAg                                                       <  2:867972/140‑1 (MQ=255)                                                               ccGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCa                                                   >  2:543329/1‑140 (MQ=37)                                                                 gTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGa                                                 >  1:557720/1‑140 (MQ=37)                                                                           gCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTc                                       <  1:448015/140‑1 (MQ=255)                                                                                             aaGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATcc                     >  2:901659/1‑140 (MQ=255)                                                                                                     cACGCCGCATCCGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCg             <  2:329474/140‑1 (MQ=255)                                                                                                                cGACATCCAACGCCCGAGCCGGCTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGcc  <  1:901659/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                      TGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCC  >  NC_000913/2304336‑2304586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑