Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,952,237	G→T	I75I (ATC→ATA)	recD ←	exodeoxyribonuclease V subunit RecD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,952,237	0	G	T	100.0%	95.5 / NA	28	I75I (ATC→ATA)	recD	exodeoxyribonuclease V subunit RecD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/13);  total (15/13)

CCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCCGCGAGGGTGACGGCAGGATGTTCA  >  NC_000913/2952104‑2952372                                                                                                                                      |                                                                                                                                        ccACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCa                                                                                                                                   >  2:577021/1‑140 (MQ=255)  cACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCAc                                                                                                                                  <  1:386401/140‑1 (MQ=255)     tgtgCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGa                                                                                                                               >  2:30195/1‑140 (MQ=255)         cgcTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCA‑CTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAg                                                                                                                          <  1:686023/140‑1 (MQ=255)                              ttCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGa                                                                                                      <  1:252070/140‑1 (MQ=255)                                     aaaGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCg                                                                                               <  1:269063/140‑1 (MQ=255)                                                            atcGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCgg                                                                        >  2:38476/1‑140 (MQ=255)                                                             tcGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGc                                                                       >  1:308324/1‑140 (MQ=255)                                                                  gTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGcaaaca                                                                  >  1:620765/1‑140 (MQ=255)                                                                   tGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGcaaacaa                                                                 <  1:577021/140‑1 (MQ=255)                                                                           tctcCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTcc                                                         <  1:15728/140‑1 (MQ=255)                                                                           tctcCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTcc                                                         <  1:690651/140‑1 (MQ=255)                                                                               cccTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCctct                                                     >  1:149603/1‑140 (MQ=255)                                                                                ccTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGcaaa                                                                    <  2:117667/124‑1 (MQ=255)                                                                                ccTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGcaaa                                                                    >  1:117667/1‑124 (MQ=255)                                                                                         cGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGcatcat                                           >  1:528356/1‑140 (MQ=255)                                                                                         cGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGcatcat                                           >  2:87692/1‑140 (MQ=255)                                                                                                        caagcaTTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCTGGCATCATGACTTAACAGTgccg                            >  2:60892/1‑140 (MQ=255)                                                                                                        caagcaTTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTgccg                            >  2:668722/1‑140 (MQ=255)                                                                                                             aTTCTTCCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCcgcg                       <  2:528356/140‑1 (MQ=255)                                                                                                                  tcCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCCGCGAGGGt                  <  1:38476/140‑1 (MQ=255)                                                                                                                  tcCCAATTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCCGCGAGGGt                  <  2:308324/140‑1 (MQ=255)                                                                                                                      aaTTTTGTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCCGCGAGGGTGACg              <  1:398473/140‑1 (MQ=255)                                                                                                                            gTAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCCGCGAGGGTGACGGCAGGa        <  1:87692/140‑1 (MQ=255)                                                                                                                             tAGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCCGCGAGGGTGACGGCAGGAt       >  2:707965/1‑140 (MQ=255)                                                                                                                              aGCTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCCGCGAGGGTGACGGCAGGATg      >  2:793698/1‑140 (MQ=255)                                                                                                                                cTCACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCCGCGAGGGTGACGGCAGGATGtt    >  1:157057/1‑140 (MQ=255)                                                                                                                                  cACCTATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCCGCGAGGGTGACGGCAGGATGTTCa  >  1:573053/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                        CCACTGTGCGCTCGTTACACCACATGCGATTCAAATAAAGACGATCGCCACAGAGGATCATCGGCGTGGGTTCATCTCCCCTGCTGACCGCTTGAGAAGCCAGCAAGCATTCTTCCCAATTTTGTAGCTCACCGATTTCACTGACACAGGTCGCCAACAGCGGATGCGACGCCTCGTTATTTTCCAGTCGTGAAAGCGGCAAACAAACGTGTCCCTCTCCGGCATCATGACTTAACAGTGCCGCCGCGAGGGTGACGGCAGGATGTTCA  >  NC_000913/2952104‑2952372

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