Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,824,719	T→C	M144V (ATG→GTG)	kguT ←	2‑ketogluconate transporter, putative

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,824,719	0	T	C	100.0%	199.4 / NA	54	M144V (ATG→GTG)	kguT	2‑ketogluconate transporter, putative
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/29);  total (25/29)

AGGCGCCACCAGAAGAACGCCCAGATCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCATCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCGATCAGCATGTAGACGTTGC  >  NC_002947/3824586‑3824850                                                                                                                                      |                                                                                                                                    aGGCGCCACCAGAAGAACGCCCAGATCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCa                                                                                                                                   <  2:357249/136‑1 (MQ=255) aGGCGCCACCAGAAGAACGCCCAGATCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCa                                                                                                                                   <  1:64731/136‑1 (MQ=255)     gcCACCAGAAGAACGCCCAGATCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGaa                                                                                                                               >  2:209972/1‑136 (MQ=255)     gcCACCAGAAGAACGCCCAGATCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGaa                                                                                                                               >  2:882526/1‑136 (MQ=255)     gcCACCAGAAGAACGCCCAGATCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGaa                                                                                                                               >  2:301722/1‑136 (MQ=255)      ccaccaGAAGAACGCCCAGATCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTTATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAg                                                                                                                              <  2:895746/136‑1 (MQ=255)      ccaccaGAAGAACGCCCAGATCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAg                                                                                                                              <  2:763016/136‑1 (MQ=255)                       aGATCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACg                                                                                                             >  1:755405/1‑136 (MQ=255)                          tCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTc                                                                                                          <  2:276074/136‑1 (MQ=255)                          tCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTc                                                                                                          <  1:625203/136‑1 (MQ=255)                                  ggcaggcCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCggg                                                                                                  <  2:812212/136‑1 (MQ=255)                                      ggcCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGaa                                                                                              >  1:1285365/1‑136 (MQ=255)                                      ggcCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGaa                                                                                              <  1:970130/136‑1 (MQ=255)                                        ccTTCGATGATGAACATCCAGCGCCTGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAAcc                                                                                            <  1:572314/136‑1 (MQ=255)                                          ttCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCaa                                                                                          <  2:949729/136‑1 (MQ=255)                                           tCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAAt                                                                                         <  2:872654/136‑1 (MQ=255)                                              atgatgAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGc                                                                                      <  2:45534/136‑1 (MQ=255)                                                   gaACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAgg                                                                                 <  1:301722/136‑1 (MQ=255)                                                   gaACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAgg                                                                                 <  1:209972/136‑1 (MQ=255)                                                     aCATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTa                                                                               >  2:702454/1‑136 (MQ=255)                                                            gcgcCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGc                                                                        >  1:624094/1‑136 (MQ=255)                                                            gcgcCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGc                                                                        >  1:726753/1‑136 (MQ=255)                                                            gcgcCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGc                                                                        >  1:993746/1‑136 (MQ=255)                                                            gcgcCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGc                                                                        >  1:1112527/1‑136 (MQ=255)                                                              gcCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTAGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCACGGGAGAACCAATGGCACAGGTACACCAGCaa                                                                      >  2:467610/1‑135 (MQ=255)                                                                    cGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCgg                                                                >  1:307400/1‑136 (MQ=255)                                                                    cGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCgg                                                                >  2:1132816/1‑136 (MQ=255)                                                                        aTGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATg                                                            <  2:83136/136‑1 (MQ=255)                                                                        aTGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATg                                                            <  2:980502/136‑1 (MQ=255)                                                                        aTGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATg                                                            <  1:1167743/136‑1 (MQ=255)                                                                             tGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGc                                                       <  1:790135/136‑1 (MQ=255)                                                                              gATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCa                                                      <  2:1329266/136‑1 (MQ=255)                                                                                tCAGGAAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCACAGCACCAGGAAGGTGATGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCAAGACCGCAGc                                                    >  2:356534/1‑136 (MQ=255)                                                                                       aCCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTccccc                                             >  1:514505/1‑136 (MQ=255)                                                                                            acacCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTc                                                 >  2:1108548/1‑127 (MQ=255)                                                                                            acacCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTc                                                 <  1:1108548/127‑1 (MQ=255)                                                                                            acacCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGAGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCAcccc                                        >  1:823830/1‑136 (MQ=255)                                                                                                           acaGAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGAt                         <  2:1265656/136‑1 (MQ=255)                                                                                                              taaTGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGcc                      <  1:467610/135‑1 (MQ=255)                                                                                                              gAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGcc                      <  2:624094/136‑1 (MQ=255)                                                                                                              gAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGcc                      <  2:1112527/136‑1 (MQ=255)                                                                                                              gAATGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGcc                      <  2:726753/136‑1 (MQ=255)                                                                                                               aaTGGTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCg                     <  1:1077438/136‑1 (MQ=255)                                                                                                                  ggTCACCGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCCCCACCCCGAGCAGGAAGCGGATGCCGATc                  >  2:662151/1‑136 (MQ=255)                                                                                                                       ccGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGcc                      <  1:929001/127‑1 (MQ=255)                                                                                                                       ccGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGcc                      >  2:929001/1‑127 (MQ=255)                                                                                                                       ccGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCGATCAGCAt             >  2:306457/1‑136 (MQ=255)                                                                                                                       ccGGGTTGCCCAGCACCAGGAAGCTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCGATCAGCAt             >  1:534446/1‑136 (MQ=255)                                                                                                                        cGGGTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCGATCAGCATg            >  2:1268149/1‑136 (MQ=255)                                                                                                                         gggTTGCCCAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCGATCAGCATGt           <  1:702454/136‑1 (MQ=255)                                                                                                                                ccAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCGATCAGCATGTAGACGtt    <  2:823830/136‑1 (MQ=255)                                                                                                                                 cAGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCGATCAGCATGTAGACGTTg   >  2:848985/1‑136 (MQ=255)                                                                                                                                  aGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCGATCAGCATGTAGACGTTGc  >  1:402392/1‑136 (MQ=255)                                                                                                                                  aGCACCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCGATCAGCATGTAGACGTTGc  >  1:1006669/1‑136 (MQ=255)                                                                                                                                      |                                                                                                                                    AGGCGCCACCAGAAGAACGCCCAGATCACGGCAGGCAGGCCTTCGATGATGAACATCCAGCGCCAGTCGTAATGCTTGATCAGGTAACCGGACACCACCGACATCCACAGAATGGTCACCGGGTTGCCCAGCATCAGGAAGGTGTTGGCGCGTGAACGCTCGGCCCGGGTGAACCAATGGCACAGGTACACCAGCATGGCCGGCATGACCGCAGCCTCCACCACCCCGAGCAGGAAGCGGATGCCGATCAGCATGTAGACGTTGC  >  NC_002947/3824586‑3824850

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: ‑