Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,167,207	T→G	I7S (ATC→AGC)	oprB‑I →	carbohydrate‑selective porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,167,207	0	T	G	100.0%	79.1 / NA	22	I7S (ATC→AGC)	oprB‑I	carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/9);  total (13/9)

GAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCATCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCTTTTTCCAGCGAATCCAAATGGATGACCGGCGACTGGGGCGGCACCCGGACCGAGCTGCTGGA  >  NC_002947/1167079‑1167336                                                                                                                                 |                                                                                                                                  gAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGaca                                                                                                                            <  2:252297/136‑1 (MQ=255)    cgcCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTg                                                                                                                         >  2:192573/1‑136 (MQ=255)                  cAGTACTGGCGCAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCctt                                                                                                           >  1:143446/1‑136 (MQ=255)                  cAGTACTGGCGCAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCctt                                                                                                           >  1:67118/1‑136 (MQ=255)                  cAGTACTGGCGCAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCctt                                                                                                           >  1:309569/1‑136 (MQ=255)                  cAGTACTGGCGCAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCctt                                                                                                           <  1:192151/136‑1 (MQ=255)                    gTACTGGCGCAACCGCAAAGGACAACGGGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCcttgc                                                                                                         <  2:60834/136‑1 (MQ=255)                          gcgcAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCAc                                                                                                   >  2:266800/1‑136 (MQ=255)                                  gCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGgc                                                                                           >  2:290934/1‑136 (MQ=255)                                  gCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGgc                                                                                           >  2:37431/1‑136 (MQ=255)                                                     aCTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTg                                                                        >  2:319938/1‑136 (MQ=255)                                                                 gTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCttttt                                                            <  2:309569/136‑1 (MQ=255)                                                                            tACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCTTTTTCCAGCGAATcc                                                 <  1:290934/136‑1 (MQ=255)                                                                                      aTAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCTTTTTCCAGCGAATCCAAATGGATGa                                       >  1:71338/1‑136 (MQ=255)                                                                                                      gCAAAGGGATGGAACAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCTTTTTCCAGCGAATCCAAATGGATGACCGGCGACTGGggcgg                       <  1:37431/136‑1 (MQ=255)                                                                                                                 gAACAGCGCAAACGCAGCAAGACACTCGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCTTTTTCCAGCGAATCCAAATGGATGACCGGCGACTGGGGCGGCACCCGGACCg            >  1:209362/1‑136 (MQ=255)                                                                                                                   aCAGCGCAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCTTTTTCCAGCGAATCCAAATGGATGACCGGCGACTGGGGCGGCACCCGGACCGAg          <  2:67118/136‑1 (MQ=255)                                                                                                                       cgcAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGgca                                                                                          >  1:37539/1‑52 (MQ=255)                                                                                                                       cgcAAACGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGgca                                                                                          <  2:37539/52‑1 (MQ=255)                                                                                                                          aaaCGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCTTTTTCCAGCGAATCCAAATGGATGACCGGCGACTGGGGCGGCACCCGGACCGAGCTGCTgg   <  2:80738/136‑1 (MQ=255)                                                                                                                           aaCGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCTTTTTCCAGCGAATCCAAATGGATGACCGGCGACTGGGGCGGCACCCGGACCGAGCTGCTGGa  >  1:127065/1‑136 (MQ=255)                                                                                                                           aaCGCAGCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCTTTTTCCAGCGAATCCAAATGGATGACCGGCGACTGGGGCGGCACCCGGACCGAGCTGCTGGa  >  1:285038/1‑136 (MQ=255)                                                                                                                                 |                                                                                                                                  GAGCGCCTGCACCTGGCCAGTACTGGCGCAACCGCAAAGGACAACGTGGCTCACTTCAAAGGCCGTTGACTCGTCTACACCATCAATAAGAAAAAGAGGACGCAAAGGGATGGAACAGCGCAAACGCATCAAGACACTGGGATCGTTGGCCTTGCTTGCACTTGTAGGCAGCAGCGGTACACAGGCTGCCGAGGCTTTTTCCAGCGAATCCAAATGGATGACCGGCGACTGGGGCGGCACCCGGACCGAGCTGCTGGA  >  NC_002947/1167079‑1167336

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