Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,955,071	C→A	E311D (GAG→GAT)	recB ←	exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,955,071	0	C	A	100.0%	91.0 / NA	27	E311D (GAG→GAT)	recB	exonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/14);  total (13/14)

CCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCC  >  NC_000913/2954933‑2955210                                                                                                                                           |                                                                                                                                            ccGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTg                                                                                                                                   <  5:288/149‑1 (MQ=255)        aaCATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTcg                                                                                                                            >  5:86587/1‑149 (MQ=255)             gTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCgggg                                                                                                                       >  2:99977/1‑149 (MQ=255)                               gCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTc                                                                                                      <  8:20331/148‑1 (MQ=255)                                             tttCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTaa                                                                                        <  3:99799/148‑1 (MQ=255)                                                   gcgcTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGc                                                                                 <  6:86587/149‑1 (MQ=255)                                                        aCTGTTTAGCGGATCTCAGCCAATGCGCGGGAGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCAACAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGAGCGATCTTCTAAGAAACGCTgg                                                                              >  6:142725/1‑147 (MQ=255)                                                                gcgGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAAttttt                                                                    >  6:133434/1‑149 (MQ=255)                                                                       tcAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGa                                                             >  1:63255/1‑149 (MQ=255)                                                                       tcAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGa                                                             >  3:57185/1‑149 (MQ=255)                                                                          gCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAtcttc                                                                                           <  1:80082/116‑1 (MQ=255)                                                                          gCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAtcttc                                                                                           >  2:80082/1‑116 (MQ=255)                                                                          gCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGaa                                                                         >  2:133089/1‑134 (MQ=255)                                                                          gCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGaa                                                                         <  1:133089/134‑1 (MQ=255)                                                                                           cACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACt                                                  <  2:941/140‑1 (MQ=255)                                                                                           cACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACt                                                  >  1:941/1‑140 (MQ=255)                                                                                            aCCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTTATAACTGt                                         <  2:63255/148‑1 (MQ=255)                                                                                             ccAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGttt                                       <  2:116754/149‑1 (MQ=255)                                                                                                           gACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGc                                                               >  3:69949/1‑111 (MQ=255)                                                                                                           gACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGc                                                               <  4:69949/111‑1 (MQ=255)                                                                                                                  gTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGgcg                  >  1:44506/1‑149 (MQ=255)                                                                                                                   ttCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGgcgc                 <  1:99977/149‑1 (MQ=255)                                                                                                                       gCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGAt             >  1:155350/1‑149 (MQ=255)                                                                                                                         aaGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGGAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCt           >  2:64936/1‑149 (MQ=255)                                                                                                                            cAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGt        <  2:44506/149‑1 (MQ=255)                                                                                                                                  atcgatcgCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATcc  <  4:57185/149‑1 (MQ=255)                                                                                                                                  atcgatcgCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATc   <  5:133434/148‑1 (MQ=255)                                                                                                                                           |                                                                                                                                            CCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCC  >  NC_000913/2954933‑2955210

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