Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007779	3,542,911	G→A	G168D (GGC→GAC)	rhaD →	rhamnulose‑1‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	3,542,911	0	G	A	100.0%	152.6 / NA	38	G168D (GGC→GAC)	rhaD	rhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/21);  total (17/21)

GATCGGGTGATCATGCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGGCAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGCACAAGAG  >  NC_007779/3542814‑3543008                                                                                                  |                                                                                                  gATCGGGTGATCATGCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGacag                                                                                                <  1:1401471/101‑1 (MQ=255)  aTCGGGTGATCATGCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGacagc                                                                                               <  1:1178192/101‑1 (MQ=255)        tGATCATGCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAg                                                                                         >  3:250630/1‑101 (MQ=255)               gCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCtgg                                                                                  >  1:195305/1‑101 (MQ=255)                 aCTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCtggtg                                                                                >  3:514338/1‑101 (MQ=255)                  cTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCtggtgg                                                                               >  3:714505/1‑101 (MQ=255)                   tGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTggtggt                                                                              <  1:1345748/101‑1 (MQ=255)                      cacgccacCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTAtt                                                                           >  1:1089312/1‑101 (MQ=255)                       acgccacCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTc                                                                          <  1:1509199/101‑1 (MQ=255)                             ccaACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGAt                                                                    <  3:1427942/101‑1 (MQ=255)                              caACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATg                                                                   <  3:1173373/101‑1 (MQ=255)                                  cTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGt                                                               >  3:865712/1‑101 (MQ=255)                                   tGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGtt                                                              >  1:62605/1‑101 (MQ=255)                                     aTCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTgg                                                            <  1:1204530/101‑1 (MQ=255)                                        gCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCAt                                                         >  3:816521/1‑101 (MQ=255)                                         cccTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCAtt                                                        <  3:1426617/101‑1 (MQ=255)                                           cTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCAtttt                                                      <  1:1265785/101‑1 (MQ=255)                                            tCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTg                                                     <  3:499427/101‑1 (MQ=255)                                             cACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGc                                                    <  3:710642/101‑1 (MQ=255)                                                cTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGt                                                 <  3:871105/101‑1 (MQ=255)                                                   tGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGa                                              <  1:473629/101‑1 (MQ=255)                                                   tGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGa                                              >  3:722042/1‑101 (MQ=255)                                                    gTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGAt                                             <  1:1390111/101‑1 (MQ=255)                                                      aCTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGAGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATgg                                           >  3:1319913/1‑101 (MQ=255)                                                       cTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGt                                          <  3:689276/101‑1 (MQ=255)                                                             aaCGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCgg                                    >  1:753483/1‑101 (MQ=255)                                                               cGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCa                                  <  1:1017232/101‑1 (MQ=255)                                                                   acCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGa                              >  3:590071/1‑101 (MQ=255)                                                                    ccGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGAc                             <  3:1009417/101‑1 (MQ=255)                                                                     cgcgGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACg                            <  3:175857/101‑1 (MQ=255)                                                                      gcgGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGa                           >  1:275319/1‑101 (MQ=255)                                                                        ggTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTACCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGaaa                         >  3:175705/1‑101 (MQ=255)                                                                         gTCTTCACTCGCCAACTGTGGGAAGACAGAACCGGGTGTCTGGTGGAATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGCGCCCGGCACGGACGAAAt                        >  1:179771/1‑101 (MQ=255)                                                                          tCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATc                       <  3:1407015/101‑1 (MQ=255)                                                                                   gCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCg              >  1:1130675/1‑101 (MQ=255)                                                                                         tgtgGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGca        <  3:63354/101‑1 (MQ=255)                                                                                          gtgGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGcac       <  1:580159/101‑1 (MQ=255)                                                                                          gtgGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGcac       >  3:245761/1‑101 (MQ=255)                                                                                           tgGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGcaca      <  1:102955/101‑1 (MQ=255)                                                                                               aaGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGCACAagag  <  3:1494787/101‑1 (MQ=255)                                                                                                  |                                                                                                  GATCGGGTGATCATGCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGGCAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGCACAAGAG  >  NC_007779/3542814‑3543008

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