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

TATGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGTCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGAA  >  NC_000913/1735348‑1735613                                                                                                                                         |                                                                                                                                  tATGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCa                                                                                                                                 <  1:289599/139‑1 (MQ=255)   tGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACt                                                                                                                               <  2:369165/139‑1 (MQ=255)    gCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTa                                                                                                                              >  2:946392/1‑139 (MQ=255)     cAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTaa                                                                                                                             <  1:2290/139‑1 (MQ=255)       aaTTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAAcc                                                                                                                           >  1:34691/1‑139 (MQ=255)               ataaAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                   <  1:819017/139‑1 (MQ=255)               ataaAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                   <  1:580361/139‑1 (MQ=255)                    ggAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGAt                                                                                                              >  2:189386/1‑139 (MQ=255)                      agagTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTa                                                                                                            <  2:34691/139‑1 (MQ=255)                              aaTGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTAcc                                                                                                    <  2:241366/139‑1 (MQ=255)                                    aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                              <  1:189386/139‑1 (MQ=255)                                    aTTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGttt                                                                                              >  1:847247/1‑139 (MQ=255)                                     ttCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGG‑ACCGCGTTTGa                                                                                            >  2:765775/1‑139 (MQ=255)                                                                 aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                   <  2:82159/89‑1 (MQ=255)                                                                 aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGaacaac                                                                                                                   >  1:82159/1‑89 (MQ=255)                                                                 aaGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                 >  1:882580/1‑139 (MQ=255)                                                                  aGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGc                                                                <  1:309877/139‑1 (MQ=255)                                                                      gCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                 <  2:229894/134‑1 (MQ=255)                                                                      gCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                 >  1:229894/1‑134 (MQ=255)                                                                       ctctGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGG‑ACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGctct                                                          <  1:765775/139‑1 (MQ=255)                                                                              aCCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACt                                                                                <  1:152202/111‑1 (MQ=255)                                                                              aCCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACt                                                                                >  2:152202/1‑111 (MQ=255)                                                                               ccGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTg                                                   <  1:92592/139‑1 (MQ=255)                                                                               ccGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTg                                                   <  2:882580/139‑1 (MQ=255)                                                                                 gCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTg                                                                                             <  1:804994/95‑1 (MQ=255)                                                                                 gCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTg                                                                                             >  2:804994/1‑95 (MQ=255)                                                                                  cacaTTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAgaga                                                                        <  2:162313/115‑1 (MQ=255)                                                                                  cacaTTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAgaga                                                                        >  1:162313/1‑115 (MQ=255)                                                                                                 aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACg                                    <  1:727999/136‑1 (MQ=255)                                                                                                 aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACg                                    >  2:727999/1‑136 (MQ=255)                                                                                                 aCCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAg                                 >  2:446874/1‑139 (MQ=255)                                                                                                      cGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGCTTa                                                                                                            <  2:691654/59‑1 (MQ=255)                                                                                                      cGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTa                                                                                                            >  1:691654/1‑59 (MQ=255)                                                                                                      cGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGc                                                 <  2:167026/118‑1 (MQ=255)                                                                                                      cGAGTATCACTACGGCAAGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGc                                                 >  1:167026/1‑118 (MQ=255)                                                                                                                       aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAGATCAc                                                                                 <  2:526664/69‑1 (MQ=255)                                                                                                                       aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAGATCAc                                                                                 >  1:526664/1‑69 (MQ=255)                                                                                                                       aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACg                                    <  2:829722/114‑1 (MQ=255)                                                                                                                       aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACg                                    >  1:829722/1‑114 (MQ=255)                                                                                                                       aGCACCATCAGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACttt           >  2:292519/1‑139 (MQ=255)                                                                                                                                aGACTTATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGaa  <  1:617185/139‑1 (MQ=255)                                                                                                                                    ttATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                               <  2:406901/74‑1 (MQ=255)                                                                                                                                    ttATGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCa                                                               >  1:406901/1‑74 (MQ=255)                                                                                                                                      aTGCCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCg                                                                 <  2:53915/70‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  TATGCAAATTAATAATAAAGGAGAGTAGCAATGTCATTCGAATTACCTGCACTACCATATGCTAAAGATGCTCTGGCACCGCACATTTCTGCGGAAACCATCGAGTATCACTACGGCAAGCACCATCAGACTTATGTCACTAACCTGAACAACCTGATTAAAGGTACCGCGTTTGAAGGTAAATCACTGGAAGAGATTATTCGCAGCTCTGAAGGTGGCGTATTCAACAACGCAGCTCAGGTCTGGAACCATACTTTCTACTGGAA  >  NC_000913/1735348‑1735613

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: ‑