Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,439	G→A	*222* (TAG→TAA)	yieH →	phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,897,439	0	G	A	100.0%	154.6 / NA	49	*222* (TAG→TAA)	yieH	phosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/24);  total (25/24)

ACGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTT  >  NC_000913/3897303‑3897564                                                                                                                                         |                                                                                                                              aCGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAAtt                                                                                                                             >  1:124854/1‑139 (MQ=255) aCGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAAtt                                                                                                                             >  1:790202/1‑139 (MQ=255)        tATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcaca                                                                                                                      <  2:790202/139‑1 (MQ=255)          tGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                       <  1:621554/136‑1 (MQ=255)          tGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                       >  2:621554/1‑136 (MQ=255)                 gtTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCa                                                                                                             >  2:827414/1‑139 (MQ=255)                 gtTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCa                                                                                                             >  2:356780/1‑139 (MQ=255)                 cttCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAAGTCTTCACACTCCCTTcg                                                                                                             >  1:515541/2‑138 (MQ=255)                  ttCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCAc                                                                                                            >  2:384947/1‑139 (MQ=255)                      aCTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTAc                                                                                                        <  1:827414/139‑1 (MQ=255)                       cTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTCCACACTCCCTTCACTTAcc                                                                                                       <  1:424729/139‑1 (MQ=255)                          cTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCg                                                                                                    >  2:223706/1‑139 (MQ=255)                            gcgcCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                       <  1:680795/118‑1 (MQ=255)                            gcgcCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                       >  2:680795/1‑118 (MQ=255)                                   cccGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTg                                                                                           <  1:592374/139‑1 (MQ=255)                                                    tCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgtaagta                                                                          >  2:673925/1‑139 (MQ=255)                                                       ttCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAg                                                                       <  1:673925/139‑1 (MQ=255)                                                       ttCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAg                                                                       <  1:223706/139‑1 (MQ=255)                                                         cACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAggg                                                                     <  2:294094/139‑1 (MQ=255)                                                                    caccacCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATg                                                          <  1:491702/139‑1 (MQ=255)                                                                          cTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCacac                                                                                                                     <  2:280693/74‑1 (MQ=255)                                                                          cTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCacac                                                                                                                     >  1:280693/1‑74 (MQ=255)                                                                                 cATCTTTCGGAGTTACCTGAGCTGTGGAAAGCGCGTGGTTGGGATATTACGACATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgtaagtaa                                                                         <  1:456348/111‑1 (MQ=255)                                                                                 cATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGACATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgtaagtaa                                                                         >  2:456348/1‑111 (MQ=255)                                                                                    cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTc                                             <  1:391186/136‑1 (MQ=255)                                                                                    cTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTc                                             >  2:391186/1‑136 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTccc                                                                                                                 >  1:385291/1‑47 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTccc                                                                                                                 <  2:385291/47‑1 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCa                                                                                    <  1:415819/76‑1 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCa                                                                                    >  2:415819/1‑76 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                     >  1:422418/1‑139 (MQ=255)                                                                                                         tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACATACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTGACGTCGGGTAATACATCACGGACTTTAtt                     >  1:567738/1‑139 (MQ=255)                                                                                                         tgGAAAACGCGTGGTTAGGATATTACCACATAATTCTTCACACTCCCTTCACTTACCCCGCTTAACTGTGCCCTCAAGGCGAAGTAAAGGATGTTTGATATCTCAGTTTCACGTCGGGTAATACATCCCCGTCTTTAtt                     >  2:401438/1‑139 (MQ=255)                                                                                                            aaaGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGATAATACATCACGGACTTTATTTTg                  <  2:567738/139‑1 (MQ=255)                                                                                                             aaGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGc                 <  2:422418/139‑1 (MQ=255)                                                                                                                cgcgTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGt                                                                  <  2:160905/87‑1 (MQ=255)                                                                                                                cgcgTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGt                                                                  >  1:160905/1‑87 (MQ=255)                                                                                                                   gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAgg                                                                      >  2:138889/1‑80 (MQ=255)                                                                                                                   gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAgg                                                                      <  1:138889/80‑1 (MQ=255)                                                                                                                   gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGt                                                         <  1:719604/93‑1 (MQ=255)                                                                                                                   gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGt                                                         >  2:719604/1‑93 (MQ=255)                                                                                                                   gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGtttt           >  2:257095/1‑139 (MQ=255)                                                                                                                            atatTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGtt  <  2:358258/139‑1 (MQ=255)                                                                                                                               ttACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgtaag                                                                            >  2:210608/1‑62 (MQ=255)                                                                                                                               ttACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgtaag                                                                            <  1:210608/62‑1 (MQ=255)                                                                                                                                 aCGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGc                 >  2:577543/1‑119 (MQ=255)                                                                                                                                 aCGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGc                 <  1:577543/119‑1 (MQ=255)                                                                                                                                     cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                <  2:150059/52‑1 (MQ=255)                                                                                                                                     cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                >  1:150059/1‑52 (MQ=255)                                                                                                                                         |                                                                                                                              ACGCAGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTT  >  NC_000913/3897303‑3897564

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