Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,153,408	A→G	intergenic (‑310/+288)	ppc ← / ← argE	phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,153,408	0	A	G	100.0%	242.4 / NA	65	intergenic (‑310/+288)	ppc/argE	phosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/36);  total (29/36)

ACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACA  >  NC_000913/4153285‑4153533                                                                                                                            |                                                                                                                              acgCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATgcgg                                                                                                                >  2:463173/1‑139 (MQ=255)   gCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATgcggcg                                                                                                              >  1:1081847/1‑139 (MQ=255)     tttAAATAAGCGTAACTTCTGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGcgcg                                                                                                            <  2:886253/139‑1 (MQ=255)     tttAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGcgcg                                                                                                            <  1:1295859/139‑1 (MQ=255)              gCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGa                                                                                                           >  1:686070/1‑131 (MQ=255)              gCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGa                                                                                                           <  2:686070/131‑1 (MQ=255)               cGTCACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                  >  2:893386/1‑139 (MQ=255)               cGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                  <  1:527379/139‑1 (MQ=255)               cGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                  <  1:893386/139‑1 (MQ=255)                gTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                 <  1:67994/139‑1 (MQ=255)                 tAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                >  1:1138311/1‑139 (MQ=255)                   aCTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                              >  1:165376/1‑139 (MQ=255)                   aCTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                              >  2:580925/1‑139 (MQ=255)                   aCTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTAGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                 >  2:742263/1‑136 (MQ=255)                   aCTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTAGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                 <  1:742263/136‑1 (MQ=255)                        tGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTAc                                                                                         >  1:340881/1‑139 (MQ=255)                              tGTTAAAATATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                   >  2:856906/1‑139 (MQ=255)                              tGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                   >  1:235667/1‑139 (MQ=255)                               gTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGcgcg                                                                                                            <  1:164827/113‑1 (MQ=255)                               gTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGcgcg                                                                                                            >  2:164827/1‑113 (MQ=255)                                 tAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCa                                                                                <  2:235667/139‑1 (MQ=255)                                 tAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCa                                                                                <  2:340881/139‑1 (MQ=255)                                      aaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCTGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGtt                                                                           <  1:928032/139‑1 (MQ=255)                                      aaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGtt                                                                           <  1:1089154/139‑1 (MQ=255)                                        tCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTg                                                                         >  1:527359/1‑139 (MQ=255)                                         cGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGAACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCATTTGa                                                                        >  2:817893/1‑139 (MQ=255)                                               aaGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACt                                                                  <  2:1138311/139‑1 (MQ=255)                                                   aaCACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcg                                                              <  2:852287/139‑1 (MQ=255)                                                    acacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcgg                                                             <  1:127569/139‑1 (MQ=255)                                                     cacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggt                                                            <  1:1268760/139‑1 (MQ=255)                                                      acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGttt                                                                          >  1:1148404/1‑124 (MQ=255)                                                      acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGttt                                                                          <  2:1148404/124‑1 (MQ=255)                                                         aaaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGa                                                        <  1:463173/139‑1 (MQ=255)                                                         aaaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACACTATTGCGGCACTGGAGTTTGGCAGCAGTGCCAGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGa                                                        <  1:856906/139‑1 (MQ=255)                                                          aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                          <  1:772817/104‑1 (MQ=38)                                                          aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                          >  2:772817/1‑104 (MQ=38)                                                          aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGc                                                                                 <  1:139067/113‑1 (MQ=255)                                                          aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGc                                                                                 >  2:139067/1‑113 (MQ=255)                                                          aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCa                                                                                <  2:437143/114‑1 (MQ=255)                                                          aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCa                                                                                >  1:437143/1‑114 (MQ=255)                                                          aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGAt                                                       >  2:658773/1‑139 (MQ=255)                                                          aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGAt                                                       >  1:1066195/1‑139 (MQ=255)                                                              tgtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGa                                                   <  2:1081847/139‑1 (MQ=37)                                                               gtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGAt                                                  <  1:381451/139‑1 (MQ=38)                                                                       gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                >  2:583806/1‑133 (MQ=37)                                                                       gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                <  1:583806/133‑1 (MQ=25)                                                                       gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCaa                                            <  1:494464/137‑1 (MQ=21)                                                                       gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCaa                                            >  2:494464/1‑137 (MQ=37)                                                                            gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggt                                                            >  2:1115348/1‑116 (MQ=255)                                                                            gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggt                                                            <  1:1115348/116‑1 (MQ=255)                                                                            gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGc                                     >  1:874784/1‑139 (MQ=255)                                                                             aTGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCa                                    <  1:658773/139‑1 (MQ=255)                                                                               gcgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATc                                  <  2:341949/139‑1 (MQ=255)                                                                                cgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATcc                                 <  2:165376/139‑1 (MQ=255)                                                                                                  cATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGGTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTg               <  2:527359/139‑1 (MQ=255)                                                                                                     tattGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAg            <  1:817893/139‑1 (MQ=38)                                                                                                             cACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa    >  2:12226/1‑139 (MQ=255)                                                                                                               cTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGt                                                                            >  2:470450/1‑65 (MQ=25)                                                                                                               cTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGt                                                                            <  1:470450/65‑1 (MQ=255)                                                                                                               cTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACa  >  1:457285/1‑139 (MQ=255)                                                                                                               cTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACa  <  2:972698/139‑1 (MQ=37)                                                                                                                      ttGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTgg              >  2:370172/1‑120 (MQ=255)                                                                                                                      ttGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTgg              <  1:370172/120‑1 (MQ=255)                                                                                                                       tGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGtt                                                                           <  2:423374/58‑1 (MQ=255)                                                                                                                       tGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGtt                                                                           >  1:423374/1‑58 (MQ=255)                                                                                                                            |                                                                                                                              ACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACA  >  NC_000913/4153285‑4153533

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