Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,098,551	A→G	G272G (GGA→GGG)	rhaS →	DNA‑binding transcriptional activator RhaS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,098,551	0	A	G	100.0%	162.5 / NA	46	G272G (GGA→GGG)	rhaS	DNA‑binding transcriptional activator RhaS
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/27);  total (19/27)

TGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTTCT  >  NC_000913/4098473‑4098627                                                                               |                                                                             tgGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCg                                                                            <  1:810627/81‑1 (MQ=255)  ggATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCgg                                                                           >  2:810699/1‑81 (MQ=255)     ttCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGAt                                                                        <  2:810628/81‑1 (MQ=255)       cAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATgg                                                                      >  1:810690/1‑81 (MQ=255)         gCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCt                                                                    <  1:810629/81‑1 (MQ=255)          cGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCtt                                                                   <  2:810630/81‑1 (MQ=255)           gACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCttt                                                                  >  2:810720/1‑81 (MQ=255)           gACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCttt                                                                  >  2:810687/1‑81 (MQ=255)              aGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTg                                                               <  1:810631/81‑1 (MQ=255)                 aaCCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCaa                                                            >  2:810697/1‑81 (MQ=255)                 aaCCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCaa                                                            >  1:810695/1‑81 (MQ=255)                   ccACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAAt                                                           >  2:810707/1‑80 (MQ=255)                   ccACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATa                                                          <  1:810632/81‑1 (MQ=255)                      cTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAAcg                                                       >  2:810689/1‑81 (MQ=255)                       ttttCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAAcgc                                                      >  2:810710/1‑81 (MQ=255)                          tCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGa                                                    <  2:810633/80‑1 (MQ=255)                           cGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAt                                                  >  2:810708/1‑81 (MQ=255)                            gACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtc                                                 <  1:810634/81‑1 (MQ=255)                            gACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAt                                                  <  2:810635/80‑1 (MQ=255)                             aCGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtct                                                <  1:810637/81‑1 (MQ=255)                             aCGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtct                                                <  2:810636/81‑1 (MQ=255)                                       gCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGt                                       <  2:810638/80‑1 (MQ=255)                                        ccGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAt                                     <  1:810639/81‑1 (MQ=255)                                          gagaGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAttt                                   >  1:810716/1‑81 (MQ=255)                                          gagaGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAttt                                   <  2:810640/81‑1 (MQ=255)                                            gagTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGt                                 <  2:810641/81‑1 (MQ=255)                                            gagTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGt                                 <  1:810642/81‑1 (MQ=255)                                             agTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTa                                >  2:810698/1‑81 (MQ=255)                                              gTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc                               >  2:810703/1‑81 (MQ=255)                                              gTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc                               >  2:810723/1‑81 (MQ=255)                                              gTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc                               >  2:810701/1‑81 (MQ=255)                                                   aCTGGTCACCGCGTGATATTGGGCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCGTCTCAACGTATTTGTACGCCat                          <  1:810643/81‑1 (MQ=255)                                                       gTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCatat                        <  1:810644/79‑1 (MQ=255)                                                       gTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCatat                        <  2:810645/79‑1 (MQ=255)                                                          aCCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCg                    <  1:810647/80‑1 (MQ=255)                                                          aCCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCg                    <  2:810646/80‑1 (MQ=255)                                                             gcgTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGaat                 >  2:810718/1‑80 (MQ=255)                                                             gcgTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACCCCATATTGCGaata                <  2:810648/81‑1 (MQ=255)                                                               gTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGaataat              >  1:810705/1‑81 (MQ=255)                                                                tGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATc             <  2:810649/81‑1 (MQ=255)                                                                     ttCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACtt        >  1:810711/1‑81 (MQ=255)                                                                       cGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTc       <  1:810650/80‑1 (MQ=255)                                                                        gCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCg      <  1:810651/80‑1 (MQ=255)                                                                        gCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGt     >  2:810702/1‑81 (MQ=255)                                                                         ccAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGtt    <  1:810652/81‑1 (MQ=255)                                                                           aGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTtct  <  1:810653/81‑1 (MQ=255)                                                                               |                                                                             TGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTTCT  >  NC_000913/4098473‑4098627

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