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%	146.8 / NA	43	G272G (GGA→GGG)	rhaS	DNA‑binding transcriptional activator RhaS
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/26);  total (17/26)

GGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTTCTC  >  NC_000913/4098474‑4098628                                                                              |                                                                              ggATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCgg                                                                            >  1:828729/1‑81 (MQ=255) ggATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCgg                                                                            >  1:828694/1‑81 (MQ=255)   aTTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGa                                                                          >  1:828705/1‑81 (MQ=255)      cAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATg                                                                        <  2:828638/80‑1 (MQ=255)       aGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGc                                                                      >  2:828710/1‑81 (MQ=255)         cGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCtt                                                                    <  1:828640/81‑1 (MQ=255)         cGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCtt                                                                    <  1:828639/81‑1 (MQ=255)          gACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCttt                                                                   <  2:828641/81‑1 (MQ=255)           aCAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTc                                                                  <  2:828642/81‑1 (MQ=255)              gTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTg                                                                >  1:828686/1‑80 (MQ=255)              gTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGc                                                               <  2:828643/81‑1 (MQ=255)               tAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGc                                                               >  1:828708/1‑80 (MQ=255)                 aCCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCaa                                                             <  1:828644/80‑1 (MQ=255)                  ccACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATa                                                           <  1:828645/81‑1 (MQ=255)                    aCTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAAc                                                         <  2:828646/81‑1 (MQ=255)                     cTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAAcg                                                        <  1:828647/81‑1 (MQ=255)                       tttCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAAcgc                                                       >  1:828717/1‑80 (MQ=255)                         tCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGa                                                     <  1:828648/80‑1 (MQ=255)                           gACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtc                                                  >  2:828700/1‑81 (MQ=255)                           gACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtc                                                  <  1:828649/81‑1 (MQ=255)                           gACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtc                                                  >  1:828720/1‑81 (MQ=255)                            aCGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtct                                                 >  2:828726/1‑81 (MQ=255)                            aCGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtct                                                 <  2:828650/81‑1 (MQ=255)                                  tttCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTtctc                                             <  1:828653/79‑1 (MQ=255)                                  tttCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCaa                                           >  1:828711/1‑81 (MQ=255)                                  tttCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCa                                            <  2:828652/80‑1 (MQ=255)                                  tttCGCCGAGAGTTTAACTGGTCACCGCGAGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTtctc                                             <  2:828651/79‑1 (MQ=255)                                    tCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAAc                                          <  2:828654/80‑1 (MQ=255)                                     cGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACg                                         >  1:828718/1‑80 (MQ=255)                                           gagTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTg                                   >  1:828712/1‑80 (MQ=255)                                           gagTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGt                                  <  2:828655/81‑1 (MQ=255)                                            agTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTa                                 <  1:828656/81‑1 (MQ=255)                                             gTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc                                <  2:828657/81‑1 (MQ=255)                                              tttAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACg                               <  2:828658/81‑1 (MQ=255)                                                 aaCTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGTCa                            >  2:828713/1‑81 (MQ=255)                                                  aCTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCat                           <  1:828659/81‑1 (MQ=255)                                                    tGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCatat                         >  1:828676/1‑81 (MQ=255)                                                         aCCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCg                     <  2:828660/80‑1 (MQ=255)                                                            gcgTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGaat                  >  2:828707/1‑80 (MQ=255)                                                            gcgTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTCCTCAACGTATTTGTACGCCATATTGCGaat                  >  2:828709/1‑80 (MQ=255)                                                                 atatTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCaa            <  1:828661/81‑1 (MQ=255)                                                                    ttCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACtt         <  1:828662/81‑1 (MQ=255)                                                                        ccAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGtt     <  1:828663/81‑1 (MQ=255)                                                                           ggggCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTtctc  <  1:828664/81‑1 (MQ=255)                                                                              |                                                                              GGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTTCTC  >  NC_000913/4098474‑4098628

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