Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,785,283	Δ1 bp	coding (1830/2379 nt)	ppsA ←	phosphoenolpyruvate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,785,281	0	C	.	100.0%	267.9 / NA	58	coding (1832/2379 nt)	ppsA	phosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (26/32);  total (26/32)

ATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAG  >  NC_000913/1785182‑1785376                                                                                                    |                                                                                                aTGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGc                                                                                                 >  1:809886/1‑100 (MQ=255)    aTCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTGACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGccgc                                                                                              >  2:848590/1‑98 (MQ=255)     tctcAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGccgcg                                                                                             >  2:1040744/1‑97 (MQ=255)       tcAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGcc                                                                                          >  1:348518/1‑100 (MQ=255)        cAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCa                                                                                         >  1:216430/1‑100 (MQ=255)         aaCGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCAc                                                                                        <  1:772039/100‑1 (MQ=255)         aaCGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCAc                                                                                        <  2:660318/100‑1 (MQ=255)            gTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGa                                                                                     >  1:1275779/1‑100 (MQ=255)              tGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAg                                                                                   <  1:113762/100‑1 (MQ=255)              tGGTCAGTCCCATGTCGTTGCGCACACGTGTCACTGCTTCTCACTCCAGCGCGAAACAGTCGCGGAAGCGGTCGGAAACATAGCGG‑CCGCGCCACGGAAg                                                                                   <  2:200830/100‑1 (MQ=255)                  cAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGa                                                                               >  2:425570/1‑100 (MQ=255)                      cccATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCAt                                                                           >  2:1253964/1‑100 (MQ=255)                        cATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCg                                                                         <  1:1135200/100‑1 (MQ=255)                        cATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCg                                                                         >  2:281567/1‑100 (MQ=255)                        cATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCg                                                                         >  1:1170800/1‑100 (MQ=255)                          tGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCggg                                                                       <  2:770904/100‑1 (MQ=255)                              gTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTtct                                                                   <  2:65594/100‑1 (MQ=255)                              gTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTtct                                                                   >  2:289524/1‑100 (MQ=255)                                 gcgcACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCtt                                                                <  1:1288958/100‑1 (MQ=255)                                   gcACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCa                                                              >  1:1299549/1‑100 (MQ=255)                                   gcACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCa                                                              <  2:413563/100‑1 (MQ=255)                                   gcACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCa                                                              <  2:1199853/100‑1 (MQ=255)                                    cacaCGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCAt                                                             >  1:19806/1‑100 (MQ=255)                                    cacaCGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCAt                                                             >  1:19804/1‑100 (MQ=255)                                        cGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTgg                                                         >  2:1103884/1‑100 (MQ=255)                                          tttCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCt                                                       <  1:681938/100‑1 (MQ=255)                                               cTGCTTCACACTCCAGCGCGAATCAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCTACGGACGCCGAGCATCGGGTTCTCTTCATCTGGCTCg                                                     <  2:343268/97‑1 (MQ=255)                                                    tCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTc                                             >  2:411916/1‑100 (MQ=255)                                                    tCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTc                                             >  1:179976/1‑100 (MQ=255)                                                    tCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTc                                             <  1:944970/100‑1 (MQ=255)                                                       cacTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTcacc                                          <  1:58609/100‑1 (MQ=255)                                                       cacTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTcacc                                          <  2:19804/100‑1 (MQ=255)                                                       cacTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTcacc                                          <  2:19806/100‑1 (MQ=255)                                                         cTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTcaccac                                        >  1:104436/1‑100 (MQ=255)                                                         cTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTcaccac                                        >  1:84749/1‑100 (MQ=255)                                                          tCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCaccacc                                       <  1:737745/100‑1 (MQ=255)                                                          tCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCaccacc                                       <  1:425574/100‑1 (MQ=255)                                                           ccAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCg                                      >  1:929682/1‑100 (MQ=255)                                                             aGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAc                                    >  1:392580/1‑100 (MQ=255)                                                              gcgcGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAcc                                   <  2:809881/100‑1 (MQ=255)                                                              gcgcGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAcc                                   <  1:429894/100‑1 (MQ=255)                                                               cgcgAAACAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCa                                  >  1:89903/1‑100 (MQ=255)                                                                     aCAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTgg                            >  1:522214/1‑100 (MQ=255)                                                                      cAGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGc                           <  2:824156/100‑1 (MQ=255)                                                                       aGTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCCCGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCa                          <  2:843430/100‑1 (MQ=255)                                                                        gTCGCGGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGCACAGGtt                              >  2:1147784/1‑95 (MQ=255)                                                                            cgGAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttc                     >  1:178724/1‑100 (MQ=255)                                                                             ggAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcg                    <  2:1299538/100‑1 (MQ=255)                                                                             ggAAGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcg                    <  2:30939/100‑1 (MQ=255)                                                                                aGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttc                 <  1:1040754/100‑1 (MQ=255)                                                                                aGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttc                 <  1:1241005/100‑1 (MQ=255)                                                                                aGCTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttc                 <  1:337995/100‑1 (MQ=255)                                                                                  cTGTCGTAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttc                 <  2:104434/98‑1 (MQ=255)                                                                                  cTGTCGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTGCATCTGGCTCGTAACGCTCACCACCGACCCGGTTGGCATATTCGTTCGa               <  2:867026/100‑1 (MQ=255)                                                                                      cGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTa           >  2:88065/1‑100 (MQ=255)                                                                                      cGGAAACATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTa           >  1:160798/1‑100 (MQ=255)                                                                                           aCATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATc      <  2:233547/100‑1 (MQ=255)                                                                                            cATAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCa     >  1:353168/1‑100 (MQ=255)                                                                                             aTAGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCag    <  1:631979/100‑1 (MQ=255)                                                                                               aGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCagag  >  2:54764/1‑100 (MQ=255)                                                                                               aGCGG‑CCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATGTAATATCagag  <  2:687106/100‑1 (MQ=255)                                                                                                    |                                                                                                ATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAG  >  NC_000913/1785182‑1785376

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