Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,955,714	T→C	Y97C (TAC→TGC)	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,714	0	T	C	100.0%	143.0 / NA	41	Y97C (TAC→TGC)	recB	exodeoxyribonuclease V subunit RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/19);  total (22/19)

ATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGTACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGG  >  NC_000913/2955584‑2955840                                                                                                                                   |                                                                                                                               aTTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCAc                                                                                                                              >  2:569377/1‑133 (MQ=255)                 cgcTGGCAAAAGCCGTGAATAGTAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCggtgg                                                                                                              >  2:421997/1‑133 (MQ=255)                  gcTGGCAAAACCCGTCAATAGTAACGACTGCCGCTTCATCCCTCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCggtggt                                                                                                             <  2:162754/133‑1 (MQ=255)                      gCAAAAGCCGTGAATAGTAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCa                                                                                                         <  1:421997/133‑1 (MQ=255)                                tGAATAGTAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAgg                                                                                               >  2:66/1‑133 (MQ=255)                                tGAATAGTAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAgg                                                                                               >  2:399603/1‑133 (MQ=255)                                tGAATAGTAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAgg                                                                                               <  1:23839/133‑1 (MQ=255)                                tGAATAGTAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGGGGTTTCACGCAGACAgg                                                                                               >  1:154940/1‑133 (MQ=255)                                  aaTAGTAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCg                                                                                             <  1:399603/133‑1 (MQ=255)                                       tAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATgcg                                                                                        <  2:602/133‑1 (MQ=255)                                        aaaGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATgcgc                                                                                       <  1:494495/133‑1 (MQ=255)                                            aCTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTACACGCAGATAGGCCCTGCGCaagt                                                                                   >  2:294464/1‑131 (MQ=255)                                              tGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATgcgc                                                                                       >  2:156492/1‑127 (MQ=255)                                              tGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATgcgc                                                                                       <  1:156492/127‑1 (MQ=255)                                              tGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCg                                                                                 >  1:586116/1‑133 (MQ=255)                                                ccGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCGGACTGGCGATGCGTAACTCGTg                                                                               >  1:291831/1‑133 (MQ=255)                                                ccGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTg                                                                               >  2:439934/1‑133 (MQ=255)                                                ccGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTg                                                                               >  1:357065/1‑133 (MQ=255)                                                ccGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTg                                                                               >  1:35281/1‑133 (MQ=255)                                                  gCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGa                                                                             >  2:593924/1‑133 (MQ=255)                                                         ccatcTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCt                                                                      <  1:440306/133‑1 (MQ=255)                                                         ccatcTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCt                                                                      <  1:593924/133‑1 (MQ=255)                                                                ccGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATa                                                               >  1:611133/1‑133 (MQ=255)                                                                     cGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGAcc                                                          >  1:45242/1‑133 (MQ=255)                                                                                     cTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACGCGTGGATATTGCGACGGGTACGACCTCGCACATCTGCCGTg                                          >  2:300610/1‑133 (MQ=255)                                                                                     cTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCACAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTg                                          >  2:143576/1‑133 (MQ=255)                                                                                      tGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGTTACGGATGCGACCGTGCAATTATGCCGTgg                                         >  1:234065/1‑133 (MQ=255)                                                                                        cgcgGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCa                                       >  2:221318/1‑133 (MQ=255)                                                                                          cgGCTTGCTCTTTATAGTAGATCTCTTCCAACAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGc                                     <  1:294464/133‑1 (MQ=255)                                                                                          cgGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGTGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGc                                     <  2:465352/133‑1 (MQ=255)                                                                                          cgGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGc                                     <  1:221318/133‑1 (MQ=255)                                                                                          cgGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGc                                     <  1:501268/133‑1 (MQ=255)                                                                                          cgGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGc                                     <  2:586116/133‑1 (MQ=255)                                                                                              ttGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCg                                 <  2:45242/133‑1 (MQ=255)                                                                                                 cgcTTTATCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTa                              <  2:35281/133‑1 (MQ=255)                                                                                                       aTCGTCGATCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTg                        <  1:569377/133‑1 (MQ=255)                                                                                                              aTCTCTTCCAGCAGGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACTGCGCAATTTTGCGGAGCTAGCCTCGATAAATTTGAccacca                 >  2:296599/1‑133 (MQ=255)                                                                                                               tctcTTCCAGCAGGCGTTCGCACACTGGTTTGGCGGTGGTTGCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACcag                <  2:234065/133‑1 (MQ=255)                                                                                                                    tCCAGCAGCCGTTCGCACCGTGAATTGTCGCTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGtt           <  1:296599/133‑1 (MQ=255)                                                                                                                         cagGCGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTGGCTAAAGGTGACCACCAGTAGttcttaa      >  2:288528/1‑131 (MQ=255)                                                                                                                             cGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACgg  >  1:418454/1‑133 (MQ=255)                                                                                                                             cGTTCGCACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACgg  >  1:1208/1‑133 (MQ=255)                                                                                                                                   |                                                                                                                               ATTCAGGTTGAGCATGCGCTGGCAAAAGCCGTGAATAGTAAAGACTGCCGCTTCATCCATCTGCCGTTCGGCTAACAACAACCACTGCGCGGCTTGCGCTTTATCGTCGATCTCTTCCAGCAGGCGTTCGTACAGTGGATTGTCGGTGGTTTCACGCAGACAGGCGATGCGCAACTCGTGGATATTGCTACGGATACGACCGCGCAATTCTGCCGTGGCAGCCTCGGTAAAGGTGACCACCAGCAGTTCTTCAACGG  >  NC_000913/2955584‑2955840

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