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

TGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGTCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGAACTGCC  >  NC_000913/1735350‑1735618                                                                                                                                       |                                                                                                                                       tGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                      >  1:949086/1‑137 (MQ=255)  gCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTa                                                                                                                                   >  1:178956/1‑139 (MQ=255)                aaGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                        >  2:400407/1‑136 (MQ=255)                aaGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                        <  1:400407/136‑1 (MQ=255)                      agTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTa                                                                                                                                   <  2:724394/119‑1 (MQ=255)                      agTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTa                                                                                                                                   >  1:724394/1‑119 (MQ=255)                      agTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTaaa                                                                                                               >  2:1139142/1‑139 (MQ=255)                      agTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCAGCAGACTTATGCCACTAACCTGAACAACCTGATTaaa                                                                                                               >  1:844385/1‑139 (MQ=255)                        tAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAgg                                                                                                             >  2:1139188/1‑139 (MQ=255)                        tAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAgg                                                                                                             >  1:764743/1‑139 (MQ=255)                                  aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaac                                                                                                                           <  1:707096/115‑1 (MQ=255)                                  aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaac                                                                                                                           >  2:707096/1‑115 (MQ=255)                                        aTTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGt                                                                                             <  2:224648/139‑1 (MQ=255)                                                 aCTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAc                                                                                                          >  1:305866/1‑117 (MQ=255)                                                 aCTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAc                                                                                                          <  2:305866/117‑1 (MQ=255)                                                          tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                      <  2:491525/80‑1 (MQ=255)                                                          tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                      >  1:491525/1‑80 (MQ=255)                                                          tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAAAAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGatt                                                                           >  2:961696/1‑139 (MQ=255)                                                           gCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTaa                                                                                           <  2:988185/122‑1 (MQ=255)                                                           gCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTaa                                                                                           >  1:988185/1‑122 (MQ=255)                                                             tAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGattatt                                                                        <  1:961696/139‑1 (MQ=255)                                                               aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                      <  2:764743/139‑1 (MQ=255)                                                               aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                      >  1:615901/1‑139 (MQ=255)                                                                    gCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGct                                                                 >  2:727592/1‑139 (MQ=255)                                                                    gCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGct                                                                 >  1:622050/1‑139 (MQ=255)                                                                     ctctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGc                                                                     >  2:588476/1‑134 (MQ=255)                                                                     ctctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGc                                                                     <  1:588476/134‑1 (MQ=255)                                                                     ctctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGACTAAAGGTACcg                                                                                                        <  1:873079/99‑1 (MQ=255)                                                                     ctctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGACTAAAGGTACcg                                                                                                        >  2:873079/1‑99 (MQ=255)                                                                      tctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                                   >  1:326417/1‑103 (MQ=255)                                                                      tctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                                   <  2:326417/103‑1 (MQ=255)                                                                          gCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAg                                                           <  2:844385/139‑1 (MQ=255)                                                                           cACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAgg                                                          <  1:1294403/139‑1 (MQ=255)                                                                               gCACATTTCTGCGGAAACCAGCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGc                                                      >  2:791955/1‑139 (MQ=255)                                                                                       cTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTcaa                                              <  1:1139142/139‑1 (MQ=255)                                                                                               aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAAcc                                                                                                                       >  1:1057773/1‑58 (MQ=255)                                                                                               aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAAcc                                                                                                                       <  2:1057773/58‑1 (MQ=255)                                                                                               aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                                    <  2:1267429/109‑1 (MQ=255)                                                                                               aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                                    >  1:1267429/1‑109 (MQ=255)                                                                                                       gTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAagag                                                                              >  1:692691/1‑91 (MQ=255)                                                                                                       gTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAagag                                                                              <  2:692691/91‑1 (MQ=255)                                                                                                       gTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTc                              <  2:178956/139‑1 (MQ=255)                                                                                                         aTCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGtttaccttc                                                                                             <  2:1174168/74‑7 (MQ=255)                                                                                                         aTCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGtttaccttc                                                                                             >  1:1174168/1‑68 (MQ=255)                                                                                                           cACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGa                          >  1:785029/1‑139 (MQ=255)                                                                                                             cTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                                    >  2:1112332/1‑95 (MQ=255)                                                                                                             cTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                                    <  1:1112332/95‑1 (MQ=255)                                                                                                              tACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAAcc                       <  1:914221/139‑1 (MQ=255)                                                                                                              tACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAAcc                       >  2:914221/1‑139 (MQ=255)                                                                                                                  gCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGctc                                                                <  2:1163251/94‑1 (MQ=255)                                                                                                                  gCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGctc                                                                >  1:1163251/1‑94 (MQ=255)                                                                                                                     aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                                   <  2:61398/56‑1 (MQ=255)                                                                                                                     aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                                   >  1:61398/1‑56 (MQ=255)                                                                                                                     aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACttt                >  2:1226283/1‑139 (MQ=255)                                                                                                                     aGCACCATCAGACTTATGCCACTAACCTGAACAACCTAATTa                                                                                                                 >  1:900024/1‑42 (MQ=255)                                                                                                                     aGCACCATCAGACTTATGCCACTAACATGAACAACCTGATTa                                                                                                                 <  2:900024/42‑1 (MQ=255)                                                                                                                         ccATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTAc            <  1:727592/139‑1 (MQ=255)                                                                                                                           aTCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTg          >  2:638736/1‑139 (MQ=255)                                                                                                                            tCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAg                                                           <  1:861121/89‑1 (MQ=255)                                                                                                                            tCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAg                                                           >  2:861121/1‑89 (MQ=255)                                                                                                                              aGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGaa       <  2:615901/139‑1 (MQ=255)                                                                                                                              aGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTAGaa       <  2:622050/139‑1 (MQ=255)                                                                                                                                   tATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTTGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGAACTGcc  <  2:785029/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       TGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGTCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGAACTGCC  >  NC_000913/1735350‑1735618

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