Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,784,778	T→G	T779P (ACC→CCC)	ppsA ←	phosphoenolpyruvate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,784,778	0	T	G	100.0%	183.4 / NA	52	T779P (ACC→CCC)	ppsA	phosphoenolpyruvate synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (23/29);  total (23/29)

GCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGC  >  NC_000913/1784641‑1784914                                                                                                                                          |                                                                                                                                         gcgTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGGTTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGgggt                                                                                                                                        >  1:415575/1‑140 (MQ=255)  cgTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTc                                                                                                                                       <  1:955467/140‑1 (MQ=255)   gTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTcc                                                                                                                                      <  2:1132278/140‑1 (MQ=255)        aCGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGtt                                                                                                                                 <  1:990719/140‑1 (MQ=255)                  tCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAgg                                                                                                                       >  1:112071/1‑140 (MQ=255)                      tgAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTAt                                                                                                                   <  2:222440/140‑1 (MQ=255)                          gagaTTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGAt                                                                                                               <  1:888496/140‑1 (MQ=255)                           agaTTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATc                                                                                                              >  1:29912/1‑140 (MQ=255)                            gaTTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATcc                                                                                                             >  2:95997/1‑140 (MQ=255)                            gaTTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATcc                                                                                                             <  1:359227/140‑1 (MQ=255)                             aTTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAgaga                                                                                                                           >  1:320068/1‑125 (MQ=255)                             aTTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAgaga                                                                                                                           <  2:320068/125‑1 (MQ=255)                             aTTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATccc                                                                                                            >  2:551739/1‑140 (MQ=255)                               tGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCct                                                                                                          <  2:991821/140‑1 (MQ=255)                                        ttCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCa                                                                                                 >  1:182969/1‑140 (MQ=255)                                         tCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCaa                                                                                                <  2:112071/140‑1 (MQ=255)                                           aTCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCct                                                                                                          >  1:1046643/1‑128 (MQ=255)                                           aTCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCct                                                                                                          <  2:1046643/128‑1 (MQ=255)                                           aTCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAAcc                                                                                              >  2:682252/1‑140 (MQ=255)                                           aTCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAAcc                                                                                              >  2:96676/1‑140 (MQ=255)                                                     aTCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCaaa                                                                                    <  2:29912/140‑1 (MQ=255)                                                     aTCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCaaa                                                                                    <  1:930965/140‑1 (MQ=255)                                                       cacaTAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGt                                                                                  >  2:773558/1‑140 (MQ=255)                                                                 cGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTc                                                                        >  2:498207/1‑140 (MQ=255)                                                                 cGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTc                                                                        >  1:964372/1‑140 (MQ=255)                                                                    cGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTcgga                                                                     <  2:182969/140‑1 (MQ=255)                                                                    cGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTcgga                                                                     <  2:925980/140‑1 (MQ=255)                                                                    cGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTcgga                                                                     <  2:998553/140‑1 (MQ=255)                                                                    cGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTcgga                                                                     <  1:682252/140‑1 (MQ=255)                                                                      aCTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTcggacg                                                                   <  1:627387/140‑1 (MQ=255)                                                                            cgccgcCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCt                                                             <  1:699697/140‑1 (MQ=255)                                                                                gccgGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGAcc                                                         >  2:526606/1‑140 (MQ=255)                                                                                 ccgGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAATCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCg                                                        >  2:139237/1‑140 (MQ=255)                                                                                 ccgGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCg                                                        <  1:773558/140‑1 (MQ=255)                                                                                 ccgGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCg                                                        >  1:536796/1‑140 (MQ=255)                                                                                   ggggATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAATCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCa                                                      <  1:139237/140‑1 (MQ=255)                                                                                    gggATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCaa                                                     <  1:397172/140‑1 (MQ=255)                                                                                      gATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAAt                                                   <  1:188567/140‑1 (MQ=255)                                                                                          tattttatttCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCg                                               <  2:964372/140‑1 (MQ=255)                                                                                               tatttCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACata                                          >  2:112958/1‑140 (MQ=255)                                                                                                      tcagttcagCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGccc                                   >  2:275308/1‑140 (MQ=255)                                                                                                            cagCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTttc                             >  2:70966/1‑140 (MQ=255)                                                                                                              gCCAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTttctt                           <  1:1077849/140‑1 (MQ=255)                                                                                                                cAGGCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTcg                         <  1:70966/140‑1 (MQ=255)                                                                                                                   gCTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCgg                      <  1:96676/140‑1 (MQ=255)                                                                                                                    cTTAACCAGGTTTGCACCACGGGGTCCGGGTTCAGAGACATGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGc                     >  1:932568/1‑140 (MQ=255)                                                                                                                             gTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGc            >  2:564442/1‑140 (MQ=255)                                                                                                                             gTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGc            >  1:661074/1‑140 (MQ=255)                                                                                                                             gTTTGCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGc            >  1:26898/1‑140 (MQ=255)                                                                                                                                 gCACCACGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATc        >  2:375878/1‑140 (MQ=255)                                                                                                                                   accacGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGa      <  2:661074/140‑1 (MQ=255)                                                                                                                                     cacGGGGTCCGGGTTCAGAGACATGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGAca    <  2:932568/140‑1 (MQ=255)                                                                                                                                     cacGGGGTCCGGGTTCAGAGACAGGGTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGAca    <  2:246725/140‑1 (MQ=255)                                                                                                                                       cGGGGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGAcagc  <  2:415575/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         GCGTTTAACGCAGGATGTCTGTGAAGAGATTGAATAAGTTTCATCTTCGGGGATCACATAACCCCGGCGACTAAACGCCGCCGGGGATTTATTTTATTTCTTCAGTTCAGCCAGGCTTAACCAGGTTTGCACCACGGTGTCCGGGTTCAGAGACAGGCTATCGATCCCCTCTTCCATCAACCATGCGGCAAAGTCTTCGTGGTCGGACGGACCCTGACCGCAAATCCCGACATATTTGCCCTGTTTCTTCGCGGCACGGATAGCCATCGACAGC  >  NC_000913/1784641‑1784914

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