Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,037,451	A→G	I219T (ATT→ACT)	recJ ←	ssDNA‑specific exonuclease RecJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,037,451	0	A	G	100.0%	131.5 / NA	37	I219T (ATT→ACT)	recJ	ssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/17);  total (20/17)

GGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAAT  >  NC_000913/3037322‑3037584                                                                                                                                  |                                                                                                                                      ggACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGGGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGtt                                                                                                                             >  2:1086901/1‑140 (MQ=255)    cGGCCCTTTCCGGCTCGGATGCGCCTCACCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACg                                                                                                                          <  2:793086/140‑1 (MQ=255)    cGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACg                                                                                                                          <  2:548094/140‑1 (MQ=255)         cTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCtcat                                                                                                                     <  1:129469/140‑1 (MQ=255)                         cGACTCACCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATc                                                                                                     <  2:200265/140‑1 (MQ=255)                                      gCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTg                                                                                        >  2:38716/1‑140 (MQ=255)                                          ggTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCa                                                                                    <  1:235317/140‑1 (MQ=255)                                          ggTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCa                                                                                    >  2:451492/1‑140 (MQ=255)                                          ggTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCa                                                                                    <  1:451492/140‑1 (MQ=255)                                           gTCAGAATGCGATTATTGGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAg                                                                                   >  2:538884/1‑140 (MQ=255)                                           gTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAg                                                                                   >  1:512987/1‑140 (MQ=255)                                           gTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAg                                                                                   >  1:769920/1‑140 (MQ=255)                                              aGAATGCGATTATTAGCGCCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAgcgc                                                                                <  2:410778/140‑1 (MQ=255)                                                 aTGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAg                                                                             >  1:814824/1‑140 (MQ=255)                                                 aTGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAg                                                                             >  1:984300/1‑140 (MQ=255)                                                   gCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCa                                                                           <  2:984300/140‑1 (MQ=255)                                                     gATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATc                                                                         >  1:993861/1‑140 (MQ=255)                                                           taGCGCCCCGCGGCCCGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATaa                                                                   <  2:769920/140‑1 (MQ=255)                                                               gTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACg                                                               <  2:814824/140‑1 (MQ=255)                                                               gTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACg                                                               >  1:1114089/1‑140 (MQ=255)                                                                tCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGc                                                              >  2:286439/1‑140 (MQ=255)                                                                     cGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCacac                                                         >  2:778906/1‑140 (MQ=255)                                                                            acgTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGcc                                                  <  1:1086901/140‑1 (MQ=255)                                                                                 cGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCa                                             >  2:485624/1‑140 (MQ=255)                                                                                           cccAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGa                                   <  2:1114089/140‑1 (MQ=255)                                                                                              aGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACgg                                <  1:538884/140‑1 (MQ=255)                                                                                                   gACCAGATCCAGCAGTTCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAAt                           >  2:538534/1‑140 (MQ=255)                                                                                                                  ttCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGt            <  1:38716/140‑1 (MQ=255)                                                                                                                   tCTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGtt           >  1:33778/1‑140 (MQ=255)                                                                                                                    cTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCgcgc                <  2:928234/134‑1 (MQ=255)                                                                                                                    cTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCgcgc                >  1:928234/1‑134 (MQ=255)                                                                                                                    cTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTa          <  1:538534/140‑1 (MQ=255)                                                                                                                    cTGCCAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTa          >  2:897372/1‑140 (MQ=255)                                                                                                                       ccAGGTTAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAggg       >  2:376309/1‑140 (MQ=255)                                                                                                                            ttAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGt      <  1:334670/136‑1 (MQ=255)                                                                                                                            ttAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGt      >  2:334670/1‑136 (MQ=255)                                                                                                                            ttAGGAGTTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAAt  >  2:681453/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                      GGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGGGTTAAT  >  NC_000913/3037322‑3037584

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