Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,955,096	T→C	E303G (GAA→GGA)	recB ←	exodeoxyribonuclease V subunit RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,955,096	0	T	C	100.0%	149.1 / NA	45	E303G (GAA→GGA)	recB	exodeoxyribonuclease V subunit RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/18);  total (27/18)

CACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑TCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACG  >  NC_000913/2954965‑2955229                                                                                                                                      |                                                                                                                                      cACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCg                                                                                                                               >  2:519553/1‑140 (MQ=255) cACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCg                                                                                                                               >  2:40157/1‑140 (MQ=255)            ttttCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCg                                                                                                                    <  1:519553/140‑1 (MQ=255)              ttCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCg                                                                                                                    >  2:786185/1‑138 (MQ=255)              ttCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCg                                                                                                                    <  1:786185/138‑1 (MQ=255)                        aCTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGGTGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGaa                                                                                                        >  1:435793/1‑140 (MQ=255)                          tGTTTCGCGGATCTCTGCCAATCCGCGGGGGATCACCCGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAAc                                                                                                      <  1:382378/140‑1 (MQ=255)                              tCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCGTCTTGCGATTTTCTAAGAAACGCTg                                                                                                  >  1:377961/1‑140 (MQ=255)                                 cgGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTCCCCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTgg                                                                                                 >  1:10589/1‑140 (MQ=255)                                 cgGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCTCCCCGGCGTTCGTGGGATCTTGTAAGAAACGCTGGga                                                                                               >  1:222518/1‑140 (MQ=255)                                 cgGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGga                                                                                               >  1:747922/1‑140 (MQ=255)                                              aTGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGACCTTCGTGCGATCCTCTAAGAAACGCTGGGAGAATTTTTCCAGc                                                                                  >  1:537470/1‑140 (MQ=255)                                                gcgcGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCGGGCCTTCGGGCGGTCTGCTAAGAAACGCTGGGAGAATTTTTCCAGCGa                                                                                >  2:321041/1‑140 (MQ=255)                                                 cgcgGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGAc                                                                               <  1:690616/140‑1 (MQ=255)                                                   cgGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTc                                                                             >  1:138992/1‑140 (MQ=255)                                                      gTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGACTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCgg                                                                          <  2:138992/140‑1 (MQ=255)                                                       tGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGc                                                                         <  1:575298/140‑1 (MQ=255)                                                                gATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATaa                                                                >  2:78468/1‑140 (MQ=255)                                                                    gcgGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGt                                                            >  1:456085/1‑140 (MQ=255)                                                                     cgGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTTGGAAAATGTTTCCAGCGACTCGGGCAACTGATAACTGtt                                                           >  2:320595/1‑140 (MQ=255)                                                                     cgGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGAt                                                                  <  1:192778/133‑1 (MQ=255)                                                                     cgGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGAt                                                                  >  2:192778/1‑133 (MQ=255)                                                                        aTCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGt                                                        >  2:498829/1‑140 (MQ=255)                                                                          cGACAATGGTTCTGCCAGCAGTTGAACGAGCGCCTCATACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtct                                                      <  1:243656/140‑1 (MQ=255)                                                                           gACAATGGTTCTGCACGCAGTTAATTGATCGCCTCAAGCAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtctc                                                     <  1:320595/140‑1 (MQ=255)                                                                           gACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACACTGGGAGAATTTGTCCGGCGACTCCAGCAAATGATAACTGTTTGtctc                                                     >  2:243656/1‑140 (MQ=255)                                                                                tGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCtcttctt                                                >  1:217230/1‑140 (MQ=255)                                                                                     cTGCAAGCAGTTGATCGATCGCATCAAACAGTGTAGGTCTCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGcc                                           <  2:377961/140‑1 (MQ=255)                                                                                      ggCAAGCAGTTGATCCATCGCCCCAAACACTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGccc                                          <  1:277099/139‑1 (MQ=255)                                                                                       gCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCa                                         <  1:78468/140‑1 (MQ=255)                                                                                        cAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGttt                                                          <  2:212335/122‑1 (MQ=255)                                                                                        cAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGttt                                                          >  1:212335/1‑122 (MQ=255)                                                                                           gCAGTTGATCGATCGCCTCAAATAGTGGACGTCGCGGTGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGgcg                                     <  1:309109/140‑1 (MQ=255)                                                                                                  atcgatcgCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACACCGGCAACAGATCACTATTTGTCTCGTCTTGTGCCCAGGCGCTGATcc                              >  2:277099/1‑139 (MQ=255)                                                                                                       tcgCCTCAAACAGTGGATGTCGCGGGGTTC‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTc                          >  2:380277/1‑140 (MQ=255)                                                                                                       tcgCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCg                         >  1:499266/1‑140 (MQ=255)                                                                                                        cgCCTCAAACAGTGGCTGTCGCCGGTTT‑‑CGCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGa                        <  1:321041/140‑1 (MQ=255)                                                                                                        cgCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGa                        <  2:782847/140‑1 (MQ=255)                                                                                                        cgCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGa                        <  1:498829/140‑1 (MQ=255)                                                                                                         gCCTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGAt                       >  2:496404/1‑140 (MQ=255)                                                                                                           cTCAAACAGTGGATGTCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATcc                     >  2:194866/1‑140 (MQ=255)                                                                                                                          tCGCGGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGc      <  1:496404/140‑1 (MQ=255)                                                                                                                            gcgGGGTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCACTTAGCCTGATTGCTa    >  1:505412/1‑140 (MQ=255)                                                                                                                              ggggTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTGTTCTGCCCCGGCGCTGATCTTGTCGATCCACTTAGCCTGATTGCTACg  >  2:373043/1‑140 (MQ=255)                                                                                                                              ggggTT‑‑CCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGGCGTTCCATTGAGCCTGATTGCTACg  >  1:341591/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                      CACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTT‑‑TCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACG  >  NC_000913/2954965‑2955229

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