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%	215.6 / 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)

TATGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGTCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTT  >  NC_000913/1735348‑1735604                                                                                                                                         |                                                                                                                         tATGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                        >  2:189805/1‑139 (MQ=255)   tGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACt                                                                                                                      <  1:1050370/139‑1 (MQ=255)          taataataaAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGa                                                                                                               >  1:165224/1‑139 (MQ=255)          taataataaAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGa                                                                                                               <  1:382349/139‑1 (MQ=255)               ataaAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                          <  2:69072/139‑1 (MQ=255)               ataaAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                          <  1:189805/139‑1 (MQ=255)                  aaGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTg                                                                                                       <  2:166018/139‑1 (MQ=255)                          tAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAgg                                                                                               >  1:163977/1‑139 (MQ=255)                             cAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAc                                                                                            >  2:155919/1‑139 (MQ=255)                               aTGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTa                                                                                                   >  2:179830/1‑130 (MQ=255)                               aTGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTa                                                                                                   <  1:179830/130‑1 (MQ=255)                               aTGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACcg                                                                                          >  2:1134004/1‑139 (MQ=255)                               aTGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACcg                                                                                          <  2:165224/139‑1 (MQ=255)                                    aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTa                                                                                                                     >  1:686991/1‑107 (MQ=255)                                    aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTa                                                                                                                     <  2:686991/107‑1 (MQ=255)                                    aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                     >  2:491181/1‑139 (MQ=255)                                    aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                     <  1:1134004/139‑1 (MQ=255)                                       cGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGaa                                                                                  <  1:264710/139‑1 (MQ=255)                                            tACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTaa                                                                             <  1:1099315/139‑1 (MQ=255)                                                      aCCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGaa                                                                   <  1:491181/139‑1 (MQ=255)                                                            tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAAcc                                                                                                         <  2:1127227/95‑1 (MQ=255)                                                            tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAAcc                                                                                                         >  1:1127227/1‑95 (MQ=255)                                                            tGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGatt                                                             >  1:296692/1‑139 (MQ=255)                                                            tGCTAAAGATGATCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGatt                                                             >  1:77584/1‑139 (MQ=255)                                                               tAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGattatt                                                          >  1:1039411/1‑139 (MQ=255)                                                                      gCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                          >  1:94728/1‑84 (MQ=255)                                                                      gCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                          <  2:94728/84‑1 (MQ=255)                                                                      gCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGTAGct                                                   >  1:748890/1‑139 (MQ=255)                                                                      gCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGct                                                   >  1:299488/1‑139 (MQ=255)                                                                        tctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                     >  1:58051/1‑103 (MQ=255)                                                                        tctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                     <  2:58051/103‑1 (MQ=255)                                                                        tctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                      <  2:504568/134‑1 (MQ=255)                                                                        tctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                      >  1:504568/1‑134 (MQ=255)                                                                          tGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGa                                               <  1:151617/139‑1 (MQ=255)                                                                           ggCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATTCCACTAACCTGAACAACCTg                                                                                                       <  2:667915/82‑1 (MQ=255)                                                                           ggCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTg                                                                                                       >  1:667915/1‑82 (MQ=255)                                                                            gCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAg                                             <  2:1039411/139‑1 (MQ=255)                                                                            gCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAg                                             <  2:299488/139‑1 (MQ=255)                                                                                                 aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAAcc                                                                                                         <  1:814974/58‑1 (MQ=255)                                                                                                 aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAAcc                                                                                                         >  2:814974/1‑58 (MQ=255)                                                                                                      cGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAg                   >  1:519744/1‑139 (MQ=255)                                                                                                         gTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTc                <  1:155919/139‑1 (MQ=255)                                                                                                                       aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAAcc         <  2:27965/132‑1 (MQ=255)                                                                                                                       aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAAcc         >  1:27965/1‑132 (MQ=255)                                                                                                                       aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACttt  >  2:978444/1‑139 (MQ=255)                                                                                                                            cATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGa                                               >  1:90806/1‑89 (MQ=255)                                                                                                                            cATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGa                                               <  2:90806/89‑1 (MQ=255)                                                                                                                             aTCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTa                                     <  1:825866/98‑1 (MQ=255)                                                                                                                             aTCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTa                                     >  2:825866/1‑98 (MQ=255)                                                                                                                              tCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                        >  2:611948/1‑78 (MQ=255)                                                                                                                              tCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                        <  1:611948/78‑1 (MQ=255)                                                                                                                              tCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                      >  2:567329/1‑80 (MQ=255)                                                                                                                              tCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                      <  1:567329/80‑1 (MQ=255)                                                                                                                              tCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACg                           >  2:1024293/1‑107 (MQ=255)                                                                                                                              tCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACg                           <  1:1024293/107‑1 (MQ=255)                                                                                                                               cAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAgg                                                                                >  1:942712/1‑53 (MQ=255)                                                                                                                               cAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAgg                                                                                <  2:942712/53‑1 (MQ=255)                                                                                                                                  aCTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAAt                                                                           >  2:338022/1‑55 (MQ=255)                                                                                                                                  aCTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAAt                                                                           <  1:338022/55‑1 (MQ=255)                                                                                                                                         |                                                                                                                         TATGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGTCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTT  >  NC_000913/1735348‑1735604

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