Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,787,349	A→G	L429P (CTG→CCG)	PP_3350 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,787,349	0	A	G	100.0%	117.5 / NA	35	L429P (CTG→CCG)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/18);  total (17/18)

CATGTGTGGCCCTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGC  >  NC_002947/3787206‑3787490                                                                                                                                                |                                                                                                                                              cATGTGTGGCCCTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGtt                                                                                                                                           <  2:234810/148‑1 (MQ=255) cATGTGTGGCCCTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTg                                                                                                                                          <  8:144638/149‑1 (MQ=255)  aTGTGTGGCCCTCTGTTGTTGTAATCGGTCGGGTGGTCAATGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCGGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTgg                                                                                                                                         <  5:111500/149‑1 (MQ=255)     tgtgGCCCTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGa                                                                                                                                       <  6:281304/148‑1 (MQ=255)          cccTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGAGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGc                                                                                                                                  >  1:82428/1‑148 (MQ=255)                gttgttgttATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAg                                                                                                                           >  1:179643/1‑149 (MQ=255)                      gttATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGt                                                                                                                      <  2:82428/148‑1 (MQ=255)                             gtcggGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGgcgc                                                                                                              <  6:226544/149‑1 (MQ=255)                             gtcggGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCACGCGCTGgctt                                                                                                              <  2:179643/149‑3 (MQ=255)                                  gtggtCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCgg                                                                                                            >  8:144437/1‑146 (MQ=255)                                              ggAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACa                                                                                             >  5:195723/1‑149 (MQ=255)                                                aaCGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATg                                                                                           <  7:144437/149‑1 (MQ=255)                                                                agaCAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCgg                                                                           <  6:195723/149‑1 (MQ=255)                                                                      acgacgAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTg                                                                     >  7:220791/1‑149 (MQ=255)                                                                          cgaATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCgg                                                                                                            >  8:299083/1‑106 (MQ=255)                                                                          cgaATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCgg                                                                                                            <  7:299083/106‑1 (MQ=255)                                                                                  cgccgcCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGa                                                                                               <  3:180307/111‑1 (MQ=255)                                                                                  cgccgcCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGa                                                                                               >  4:180307/1‑111 (MQ=255)                                                                                    ccgccAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCACGCGCTGgctt                                                                                                              <  7:339188/94‑3 (MQ=255)                                                                                    ccgccAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCACGCGCTGgctt                                                                                                              >  8:339188/1‑92 (MQ=255)                                                                                         aGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGGTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGgt                                                                         <  1:258510/126‑1 (MQ=255)                                                                                         aGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGgt                                                                         >  2:258510/1‑126 (MQ=255)                                                                                         aGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCgg                                                  <  2:111388/149‑1 (MQ=255)                                                                                           cctcctTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCggg                                                 >  5:92379/1‑148 (MQ=255)                                                                                                  gATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTaca                                         >  6:196251/1‑149 (MQ=255)                                                                                                                  tcgtAGTACTTCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACg                         <  5:196251/149‑1 (MQ=255)                                                                                                                           ttCAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGAGTGCTGGAGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCca                >  7:76156/1‑149 (MQ=255)                                                                                                                             cAGGTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCcaca              >  3:168909/1‑149 (MQ=255)                                                                                                                                gTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATg                                             >  4:144236/1‑115 (MQ=255)                                                                                                                                gTCAAGGTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATg                                             <  3:144236/115‑1 (MQ=255)                                                                                                                                      gTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTc                                  >  6:330960/1‑120 (MQ=255)                                                                                                                                      gTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTc                                  <  5:330960/120‑1 (MQ=255)                                                                                                                                      gTCCAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCg         >  8:104459/1‑145 (MQ=255)                                                                                                                                         cAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGc  <  6:92379/149‑1 (MQ=255)                                                                                                                                         cAGGGTCGGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGc  >  5:170483/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                              CATGTGTGGCCCTCTGTTGTTGTTATCGGTCGGGTGGTCAAAGGCGGAACGTGGCGCGTACTGAGACAAACGACGAATCTTCGCCGCCAGCCTCCTTGATGGCCGTGCCCACGTCGTAGTACTTCAGGTCAAGGTCCAGGGTCAGGTTGGGAGTGGCAAGCCAGCGCGTGCTGGCGCGGTAGATGGTACCGACATGGCGACCCGTGCCGCGGGTGTTGGCCAGCATGCTCATGCCGGGCATGTACACGCCGTCTTCCTCACGCTGTTTCCACACTTCGAAGATGC  >  NC_002947/3787206‑3787490

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