Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,955,349	C→T	A219T (GCT→ACT)	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,349	0	C	T	100.0%	50.8 / NA	18	A219T (GCT→ACT)	recB	exodeoxyribonuclease V subunit RecB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/8);  total (10/8)

CGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTC  >  NC_000913/2955103‑2955587                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               cGGCCTTCGTGCGATCTCCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCAATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGCAGTCAGc                                                                                                                                                                                                                                            >  1:69136/1‑251 (MQ=255)                                           ccAGCGACTCCGGCAACTGATCACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCg                                                                                                                                                                                                  >  1:408871/1‑251 (MQ=255)                                                                         nccnttcntnnGNNCAGNCGCTGATCTTGTCGATCCATTNNGCCTGATTGNTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTNCACCCACTGNNTCNCGCCANTNNNGNNNTANCGTATCAATACGCGCCACANTTNNCGCGTGNCGGGAAGTCAGCGTTTCATCATCggg                                                                                                                                                                                                                              <  2:117014/189‑1 (MQ=255)                                                                                                        gATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCATCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTTTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGTCGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGa                                                                                                                                                                                                                <  1:335055/176‑1 (MQ=255)                                                                                                                                gTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCCCAATTTGCGCGTGACGGGAAGTCAGCGTTTCATAATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGc                                                                                                              >  1:325368/1‑250 (MQ=255)                                                                                                                                     aCTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCggg                                                                                                        <  1:33945/251‑1 (MQ=255)                                                                                                                                                                        gcATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGTGTTGATAACCGGCGCTGCGTCTTGCAGTTAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTccn                                                                                                                                                <  1:57768/176‑2 (MQ=255)                                                                                                                                                                               gTTCCCCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGAGAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGgtagcgt                                                              <  1:449703/251‑1 (MQ=255)                                                                                                                                                                                    cccACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGgtagcgtagc                                                           >  1:194315/1‑249 (MQ=255)                                                                                                                                                                                               gcgcCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCtt                                              >  1:280031/1‑251 (MQ=255)                                                                                                                                                                                                     ctgctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGTCAGCGTTTCATCCTCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATc                                         >  1:151412/1‑250 (MQ=255)                                                                                                                                                                                                          gTTTTAACGTATCAATACACGCCACAATTTGCGCGTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCATCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGGGGCGACGCCAGAAATCGGCGCAGGCCTGGGAGCGTAGCAGAGACTCATCTTCAATCAgctgc                                   >  1:115488/1‑251 (MQ=255)                                                                                                                                                                                                                    aTCAATACGCTCCACAATTTGCGCGTGACTGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCATTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCGTCAATCAGCTGCTTCTCAAAc                          <  1:390557/250‑1 (MQ=255)                                                                                                                                                                                                                     tCAATACGCGCCACAATTTGCGCGTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTCATTTGATAACCGGCGCTTCGCCTTGCAGATAACGAATAATATCGCGCAGCAAc                                                                                                                                                                 >  1:271889/1‑114 (MQ=255)                                                                                                                                                                                                                                cacaATTTGCGCGTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGTCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGGTGCAAAGACGACCTGGGCTATTTCACGCGGCAGCGTGTAGCAGTGGCGACGCGAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCAGCTTCAATCAgc                                      <  1:370543/226‑1 (MQ=255)                                                                                                                                                                                                                                      gtGCGCGTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTACGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCATCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCTTAGCAGAGACTCATCTTCAATCAGCTGCCGCTCAAACAGCATGCCGGATTCAAAg        <  1:486466/249‑1 (MQ=255)                                                                                                                                                                                                                                         ngngnNACGCGAAGNCAGCGTTTCATCATCGGGCGGCGGt                                                                                                                                                                                                                       <  2:138372/39‑1 (MQ=25)                                                                                                                                                                                                                                         cgcgTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCAt                     >  1:353683/1‑234 (MQ=255)                                                                                                                                                                                                                                           cgTGACGGGAAGTCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGCAATCGGCGCAGGCCTGGAAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTc  >  1:268020/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               CGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGACTCATCTTCAATCAGCTGCTGCTCAAACAGCATGCCGGATTCAAAGGCATTC  >  NC_000913/2955103‑2955587

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