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%	223.2 / NA	64	V36A (GTC→GCC)	sodB	superoxide dismutase, Fe
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/31);  total (33/31)

ATGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGTCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGAACTGCCT  >  NC_000913/1735349‑1735619                                                                                                                                        |                                                                                                                                        atacAAATTAAGAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTGCCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCAc                                                                                                                                      >  1:1089858/4‑139 (MQ=255)   gCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTa                                                                                                                                    <  2:584912/139‑1 (MQ=255)   gCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTa                                                                                                                                    >  1:584912/1‑139 (MQ=255)   gCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTa                                                                                                                                    >  2:1416864/1‑139 (MQ=255)    cAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTaa                                                                                                                                   <  2:1089858/139‑1 (MQ=255)                         tAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTAATTAAAga                                                                                                              >  1:554641/1‑138 (MQ=255)                             aaTGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAcc                                                                                                          <  1:1416864/139‑1 (MQ=255)                              aTGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACcg                                                                                                         >  1:1616988/1‑139 (MQ=255)                                 tCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGt                                                                                                      >  1:1808700/1‑139 (MQ=255)                                   aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGAt                                                                                                                    >  1:1044069/1‑123 (MQ=255)                                   aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGAt                                                                                                                    <  2:1044069/123‑1 (MQ=255)                                   aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCAGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAgg                                                                                                              <  1:541882/129‑1 (MQ=255)                                   aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCAGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAgg                                                                                                              >  2:541882/1‑129 (MQ=255)                                           tACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTaa                                                                                            <  1:1060411/139‑1 (MQ=255)                                             ccTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                       >  1:1636615/1‑94 (MQ=255)                                             ccTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                       <  2:1636615/94‑1 (MQ=255)                                             ccTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAAt                                                                                          <  2:1699940/139‑1 (MQ=255)                                               tGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCa                                                                                        >  2:983637/1‑139 (MQ=255)                                                 cACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACt                                                                                      >  2:1675763/1‑139 (MQ=255)                                                      ccATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAag                                                                                 <  1:1675763/139‑1 (MQ=255)                                                         tatGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAgaga                                                                              >  1:373300/1‑139 (MQ=255)                                                           tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaa                                                                                                                             >  1:364230/1‑90 (MQ=255)                                                           tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaa                                                                                                                             <  2:364230/90‑1 (MQ=255)                                                           tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGatt                                                                            <  2:1808700/139‑1 (MQ=255)                                                              tAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                       <  2:808309/77‑1 (MQ=255)                                                              tAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                       >  1:808309/1‑77 (MQ=255)                                                              tAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGattatt                                                                         <  1:983637/139‑1 (MQ=255)                                                                aaGATGCTCTGGCACCGCACATTTCTGTGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                       >  1:515321/1‑139 (MQ=255)                                                                aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGt                                                                                              <  1:446253/116‑1 (MQ=255)                                                                aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGt                                                                                              >  2:446253/1‑116 (MQ=255)                                                                aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                       <  2:554641/139‑1 (MQ=255)                                                                aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                       >  1:847740/1‑139 (MQ=255)                                                                 aGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAAt                                                                                          <  1:797505/119‑1 (MQ=255)                                                                 aGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAAt                                                                                          >  2:797505/1‑119 (MQ=255)                                                                     gCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGct                                                                  >  2:85218/1‑139 (MQ=255)                                                                      ctctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGctc                                                                 >  2:1508046/1‑139 (MQ=255)                                                                           gCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAg                                                            <  2:1616988/139‑1 (MQ=255)                                                                           gCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAg                                                            <  2:847740/139‑1 (MQ=255)                                                                              ccGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAAACTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTg                                                         >  1:291780/1‑139 (MQ=255)                                                                                gCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGc                                                       <  2:291780/139‑1 (MQ=255)                                                                                 cacaTTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTaa                                                                                                                                   >  2:930369/1‑62 (MQ=255)                                                                                 cacaTTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTaa                                                                                                                                   <  1:930369/62‑1 (MQ=255)                                                                                 cacaTTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCg                                                      >  1:1079221/1‑139 (MQ=255)                                                                                        cTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTcaa                                               <  1:1387309/139‑1 (MQ=255)                                                                                                aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGatt                                                                            >  2:1475226/1‑102 (MQ=255)                                                                                                aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGatt                                                                            <  1:1475226/102‑1 (MQ=255)                                                                                                aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAg                                       >  2:573578/1‑139 (MQ=255)                                                                                                aCCATCGAGTAACACTACGACAAGCACCATCAGACTTATCCCACTAACCTGAACAAACGGATTAAAGGTCCCGCGTTTGAAGGTAAGTCACTTGAAGAGATTATTCGCAGCCCAGAAAGTGGCGTATTCAACAACGCAg                                       >  1:936847/1‑139 (MQ=255)                                                                                                  cATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCt                                     <  2:515321/139‑1 (MQ=255)                                                                                                   aTCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTc                                    <  2:936847/139‑1 (MQ=255)                                                                                                        gTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTc                               <  1:573578/139‑1 (MQ=255)                                                                                                         tATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCt                              >  1:832991/1‑139 (MQ=255)                                                                                                             aCTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGaa                          <  1:1182779/139‑1 (MQ=255)                                                                                                             aCTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCGGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGaa                          <  1:85218/139‑1 (MQ=255)                                                                                                              cTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAAcc                                                                                                                        <  1:1529377/44‑1 (MQ=255)                                                                                                              cTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAAcc                                                                                                                        >  2:1529377/1‑44 (MQ=255)                                                                                                                      aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTc                                    <  2:1874149/120‑1 (MQ=255)                                                                                                                      aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTc                                    >  1:1874149/1‑120 (MQ=255)                                                                                                                      aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACttt                 >  1:906593/1‑139 (MQ=255)                                                                                                                            aTCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAag                                                                                 <  2:579730/69‑1 (MQ=255)                                                                                                                            aTCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAag                                                                                 >  1:579730/1‑69 (MQ=255)                                                                                                                              cAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGt                                                                                              <  2:529466/54‑1 (MQ=255)                                                                                                                              cAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGt                                                                                              >  1:529466/1‑54 (MQ=255)                                                                                                                                     aTGCCACAAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGGTTATTCGCAGCTCTGAAGGTGGCGCATTCAACAACGCAGCTCAGGTCTGGAACCATTCTTTCTACTGGAACTGcat  >  1:1229630/1‑137 (MQ=255)                                                                                                                                        |                                                                                                                                        ATGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGTCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGAACTGCCT  >  NC_000913/1735349‑1735619

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