Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,735,484	T→C	V36A (GTC→GCC)	sodB →	superoxide dismutase, Fe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,735,484	0	T	C	100.0%	214.5 / NA	59	V36A (GTC→GCC)	sodB	superoxide dismutase, Fe
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/30);  total (29/30)

ATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGTCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGAACTGCCT  >  NC_000913/1735362‑1735619                                                                                                                           |                                                                                                                                        ataaAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                         <  1:947979/139‑1 (MQ=255)          agTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTaaa                                                                                                                >  1:646728/1‑139 (MQ=255)          agTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTaaa                                                                                                                >  2:604707/1‑139 (MQ=255)          agTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTaaa                                                                                                                >  1:374452/1‑139 (MQ=255)            tAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAgg                                                                                                              >  2:878177/1‑139 (MQ=255)            tAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAgg                                                                                                              <  1:878177/139‑1 (MQ=255)               cAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAc                                                                                                           <  1:604707/139‑1 (MQ=255)               cAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAc                                                                                                           >  2:721127/1‑139 (MQ=255)                 aTGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCTTCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACcg                                                                                                         <  2:646728/139‑1 (MQ=255)                    tCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                       <  2:1092705/106‑1 (MQ=255)                    tCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAc                                                                                                           <  2:878848/134‑1 (MQ=255)                    tCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAc                                                                                                           >  1:878848/1‑134 (MQ=255)                    tCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGt                                                                                                      >  1:67619/1‑139 (MQ=255)                            aTTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGt                                                                                              >  2:72646/1‑139 (MQ=255)                                       tACCATATGCTTAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGa                                                                                   <  1:721127/139‑1 (MQ=255)                                         ccATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTa                                                                                                                  >  1:477115/1‑106 (MQ=255)                                         ccATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTa                                                                                                                  <  2:477115/106‑1 (MQ=255)                                              tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGAt                                                                                                                    <  1:998981/99‑1 (MQ=255)                                              tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGAt                                                                                                                    >  2:998981/1‑99 (MQ=255)                                              tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGatt                                                                            >  1:269226/1‑139 (MQ=255)                                                   aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                         >  1:1037062/1‑89 (MQ=255)                                                   aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                         <  2:1037062/89‑1 (MQ=255)                                                   aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGAt                                                                                                                    >  2:97541/1‑94 (MQ=255)                                                   aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGAt                                                                                                                    <  1:97541/94‑1 (MQ=255)                                                        gCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGct                                                                  >  1:602691/1‑139 (MQ=255)                                                         ctctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACt                                                                                                                                     <  2:445016/71‑1 (MQ=255)                                                         ctctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACt                                                                                                                                     >  1:445016/1‑71 (MQ=255)                                                         ctctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGctc                                                                 <  2:602691/139‑1 (MQ=255)                                                          tctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGctct                                                                >  1:1077522/1‑139 (MQ=255)                                                                 ccGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTg                                                         <  1:871002/139‑1 (MQ=255)                                                                   gCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGc                                                       <  2:177097/139‑1 (MQ=255)                                                                    cacaTTTCTGCGGAAACCATCGAGTATCGCTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                         >  1:688995/1‑72 (MQ=255)                                                                    cacaTTTCTGCGGAAACCATCGAGTATCGCTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                                    >  1:688487/1‑93 (MQ=255)                                                                    cacaTTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                         <  2:688995/72‑1 (MQ=255)                                                                    cacaTTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAc                                                                                                           >  2:310259/1‑86 (MQ=255)                                                                    cacaTTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAc                                                                                                           <  1:310259/86‑1 (MQ=255)                                                                    cacaTTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                                    <  2:688487/93‑1 (MQ=255)                                                                        tttCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                       <  2:839521/54‑1 (MQ=255)                                                                        tttCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                       >  1:839521/1‑54 (MQ=255)                                                                            tGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTcaac                                              >  1:898081/1‑139 (MQ=255)                                                                                 aaaCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGc                                         >  1:1101733/1‑139 (MQ=255)                                                                                   aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                                     >  2:456828/1‑109 (MQ=255)                                                                                   aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                                     <  1:456828/109‑1 (MQ=255)                                                                                   aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTg                                                         <  2:563790/121‑1 (MQ=255)                                                                                   aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTg                                                         >  1:563790/1‑121 (MQ=255)                                                                                     cATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCt                                     <  1:278384/139‑1 (MQ=255)                                                                                           gTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTc                               <  2:67619/139‑1 (MQ=255)                                                                                           gTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTc                               <  1:72646/139‑1 (MQ=255)                                                                                                  tACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAAcc                        <  2:374452/139‑1 (MQ=255)                                                                                                         aGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                         <  2:322311/35‑1 (MQ=255)                                                                                                         aGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                         >  1:322311/1‑35 (MQ=255)                                                                                                         aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                       <  2:885905/85‑1 (MQ=255)                                                                                                         aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                       >  1:885905/1‑85 (MQ=255)                                                                                                         aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACttt                 >  1:636094/1‑139 (MQ=255)                                                                                                                tCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                                     <  1:310012/80‑1 (MQ=255)                                                                                                                tCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                                     >  2:310012/1‑80 (MQ=255)                                                                                                                  aGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGaa        <  2:636094/139‑1 (MQ=255)                                                                                                                  aGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGaa        <  2:898081/139‑1 (MQ=255)                                                                                                                        aTGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGAACTGCCt  >  1:183393/1‑139 (MQ=255)                                                                                                                           |                                                                                                                                        ATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGTCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGAACTGCCT  >  NC_000913/1735362‑1735619

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