Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,955,841	T→C	T55A (ACC→GCC)	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,841	0	T	C	100.0%	123.5 / NA	37	T55A (ACC→GCC)	recB	exodeoxyribonuclease V subunit RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/15);  total (22/15)

GGCGTTCGTACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGTCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAGGGCAAGCGCAAAGG  >  NC_000913/2955706‑2955976                                                                                                                                        |                                                                                                                                        ggCGTTCGTACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAgc                                                                                                                                     >  2:209964/1‑140 (MQ=255)   cGTTCGTACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAgcgg                                                                                                                                   >  2:812758/1‑140 (MQ=255)       cGTACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAgcgggcgg                                                                                                                               >  2:450349/1‑140 (MQ=255)             gTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCGGTTCTTAAACGGCCAGCGGGCGGGGAAAg                                                                                                                         >  1:332714/1‑140 (MQ=255)              tGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAgg                                                                                                                        >  2:362659/1‑140 (MQ=255)                gATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAggcg                                                                                                                      >  1:607698/1‑140 (MQ=255)                   tGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGgcggcg                                                                                                                   <  2:607698/140‑1 (MQ=255)                     tCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGa                                                                                                                 >  1:640716/1‑140 (MQ=255)                      cGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGaa                                                                                                                <  1:786700/140‑1 (MQ=255)                          ggtTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGc                                                                                                            >  1:64602/1‑140 (MQ=255)                           gtTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGTCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGcc                                                                                                           <  1:641549/140‑1 (MQ=255)                           gtTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGcc                                                                                                           <  2:199109/140‑1 (MQ=255)                           gtTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGcc                                                                                                           <  1:812758/140‑1 (MQ=255)                                      aCAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTa                                                                                                >  1:629007/1‑140 (MQ=255)                                       cAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTaa                                                                                               <  1:132258/140‑1 (MQ=255)                                                cgcAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCaa                                                                                      >  1:515642/1‑140 (MQ=255)                                                 gcAACTCGTGGATATTGCTACGTATACGACCGCGCGATTGTGCCGTGGTAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCaaa                                                                                     <  2:332714/140‑1 (MQ=255)                                                             tatTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCa                                                                         >  2:91910/1‑140 (MQ=255)                                                               ttGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCGAGTAACGGGCGCAGATAGAGCGCCGCAAt                                                                       >  1:382226/1‑140 (MQ=255)                                                                   tACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTa                                                                   <  1:91910/140‑1 (MQ=255)                                                                   tACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTa                                                                   <  1:209964/140‑1 (MQ=255)                                                                          aCGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCa                                                                         >  2:167479/1‑127 (MQ=255)                                                                          aCGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCa                                                                         <  1:167479/127‑1 (MQ=255)                                                                                 cgcAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCtgt                                                     >  2:803606/1‑140 (MQ=255)                                                                                 cgcAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCtgt                                                     >  1:502603/1‑140 (MQ=255)                                                                                  gcAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCTAGTAACAGGCGCGACTAGAGCGCCGCAATCTCAAAGGTTTTGCCtgtg                                                    >  2:376629/1‑140 (MQ=255)                                                                                  gcAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTa                                                                   >  1:122770/1‑125 (MQ=255)                                                                                  gcAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTa                                                                   <  2:122770/125‑1 (MQ=255)                                                                                                gCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGACGCTATCCTAAAGGTTGTGCCAGTGCCGGCAGAGGCTTc                                      >  2:379782/1‑140 (MQ=255)                                                                                                    ccTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATc                                  >  2:506248/1‑140 (MQ=255)                                                                                                             aGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCa                         <  2:507046/140‑1 (MQ=255)                                                                                                             aGGTGACCACCAGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCa                         <  2:64602/140‑1 (MQ=255)                                                                                                                   ccaccaGCAGTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGt                   >  1:488661/1‑140 (MQ=255)                                                                                                                         gcagTTCTTCAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAGGGc             <  2:19704/140‑1 (MQ=255)                                                                                                                               cttcAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAGGGCAAgcgc       <  1:652491/140‑1 (MQ=255)                                                                                                                               cttcAACGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCACAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAGGGCAAgcgc       >  2:706030/1‑140 (MQ=255)                                                                                                                                    aCGGCCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAGCCTGTAAGGGCAAGCGCAAAgg  >  2:275130/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        GGCGTTCGTACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGGTCAGCGGGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAGGGCAAGCGCAAAGG  >  NC_000913/2955706‑2955976

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