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

TGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTTCTC  >  NC_000913/4098473‑4098628                                                                               |                                                                              tgGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCg                                                                             <  2:760871/81‑1 (MQ=255)  ggATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCgg                                                                            >  2:760931/1‑81 (MQ=255)          cGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCtt                                                                    <  1:760875/81‑1 (MQ=255)          cGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCtt                                                                    <  2:760875/81‑1 (MQ=255)           gACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCttt                                                                   >  1:760947/1‑81 (MQ=255)             cAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCt                                                                 >  2:760924/1‑81 (MQ=255)                tAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGc                                                               >  1:760949/1‑80 (MQ=255)                  aCCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAAt                                                            <  1:760877/81‑1 (MQ=255)                      cTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAAcg                                                        >  2:760945/1‑81 (MQ=255)                         ttCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAAcgcg                                                      <  1:760878/80‑1 (MQ=255)                         ttCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAAcgc                                                       <  2:760878/79‑1 (MQ=255)                          tCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGa                                                     <  1:760881/80‑1 (MQ=255)                          tCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGa                                                     <  2:760879/80‑1 (MQ=255)                            gACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtc                                                  >  2:760942/1‑81 (MQ=255)                              cGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtctt                                                >  2:760940/1‑81 (MQ=255)                               gCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtcttc                                               <  2:760881/81‑1 (MQ=255)                               gCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtcttc                                               >  2:760938/1‑81 (MQ=255)                                cTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAAtcttct                                              >  2:760936/1‑81 (MQ=255)                                    ttCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCaa                                           <  1:760884/80‑1 (MQ=255)                                    ttCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCaa                                           <  1:760883/80‑1 (MQ=255)                                    ttCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAAc                                          >  1:760935/1‑81 (MQ=255)                                       gCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGt                                        <  2:760884/80‑1 (MQ=255)                                         cGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAtt                                     <  1:760886/81‑1 (MQ=255)                                          gagaGTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAttt                                    >  2:760939/1‑81 (MQ=255)                                            gagTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTg                                   >  1:760951/1‑80 (MQ=255)                                             agTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTa                                 >  1:760931/1‑81 (MQ=255)                                              gTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc                                >  2:760929/1‑81 (MQ=255)                                              gTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc                                >  1:760934/1‑81 (MQ=255)                                              gTTTAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc                                <  2:760886/81‑1 (MQ=255)                                               tttAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACg                               >  1:760953/1‑81 (MQ=255)                                               tttAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACg                               <  1:760888/81‑1 (MQ=255)                                                 tAACTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGc                              >  2:760937/1‑80 (MQ=255)                                                  aaCTGGTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGcc                             <  1:760889/80‑1 (MQ=255)                                                       gTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCatat                         <  1:760890/79‑1 (MQ=255)                                                       gTCACCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCatat                         <  1:760891/79‑1 (MQ=255)                                                          aCCGCGTGATATTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGa                    <  2:760891/81‑1 (MQ=255)                                                                  atatTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCaa            >  1:760952/1‑81 (MQ=255)                                                                  atatTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCaa            <  1:760893/81‑1 (MQ=255)                                                                  atatTCGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCa             <  1:760894/80‑1 (MQ=255)                                                                       cGCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTc        <  2:760894/80‑1 (MQ=255)                                                                        gCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGt      <  2:760895/81‑1 (MQ=255)                                                                        gCCAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGt      >  2:760947/1‑81 (MQ=255)                                                                         ccAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGtt     <  2:760896/81‑1 (MQ=255)                                                                          cAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTtc    <  2:760897/81‑1 (MQ=255)                                                                          cAGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTtc    >  2:760932/1‑81 (MQ=255)                                                                           aGGGGCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGCACGCCATATTGCGAATAAGCAACTTCGTtct   >  2:760935/1‑81 (MQ=255)                                                                            ggggCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTtctc  <  1:760899/81‑1 (MQ=255)                                                                               |                                                                              TGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTTCTC  >  NC_000913/4098473‑4098628

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: ‑