Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,064,126	G→C	Y121* (TAC→TAG)	brpA ←	LytR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,064,126	0	G	C	100.0%	149.5 / NA	41	Y121* (TAC→TAG)	brpA	LytR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (21/20);  total (21/20)

GCTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCCTTTTTTTC  >  CP000730/1063993‑1064262                                                                                                                                      |                                                                                                                                         gCTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGc                                                                                                                                   <  2:319123/141‑1 (MQ=255)  cTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGcc                                                                                                                                  >  1:120244/1‑141 (MQ=255)      tttACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATGCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACt                                                                                                                              >  2:36589/1‑141 (MQ=255)         aCAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTa                                                                                                                           >  1:203936/1‑141 (MQ=255)            gCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGa                                                                                                                        >  2:73194/1‑141 (MQ=255)             ccTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGaa                                                                                                                       <  2:3962/141‑1 (MQ=255)              cTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAAt                                                                                                                      >  2:255836/1‑141 (MQ=255)               tCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATg                                                                                                                     >  1:85764/1‑141 (MQ=255)                cGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGt                                                                                                                    >  2:385078/1‑141 (MQ=255)                    aaaTGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCt                                                                                                                <  2:187322/141‑1 (MQ=255)                        gCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATa                                                                                                            >  1:340632/1‑141 (MQ=255)                         cTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAg                                                                                                           <  1:107447/141‑1 (MQ=255)                           ttCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTa                                                                                                         >  1:144385/1‑141 (MQ=255)                                 ttGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGt                                                                                                   <  1:251555/141‑1 (MQ=255)                                        tcACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATAc                                                                                            <  1:210779/141‑1 (MQ=255)                                            ataataATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAg                                                                                        <  1:356316/141‑1 (MQ=255)                                               ataatCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCAt                                                                                     >  1:95925/1‑141 (MQ=255)                                               ataatCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCAt                                                                                     >  1:378582/1‑141 (MQ=255)                                                     aaCAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAAt                                                                               >  2:367040/1‑141 (MQ=255)                                                     aaCAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAAt                                                                               >  1:172228/1‑141 (MQ=255)                                                     aaCAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAAt                                                                               >  1:406588/1‑141 (MQ=255)                                                            agaTTCAGGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGTGTAATGTTATGCTAATAGGCAACTTTTTGTATGTAGCTTATAGTATCACGTGGAATACGAAGCGTGCTAATTTGATGt                                                                        <  2:245151/139‑1 (MQ=255)                                                            aCATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGt                                                                        <  2:116082/141‑1 (MQ=255)                                                                  aTGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTttgat                                                                  <  1:45171/141‑1 (MQ=255)                                                                     gTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTttgattga                                                               <  2:30268/141‑1 (MQ=255)                                                                                cGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAagtagt                                                    <  1:367112/141‑1 (MQ=255)                                                                                 gAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGtt                                                   <  1:5242/141‑1 (MQ=255)                                                                                      cATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAAt                                              <  1:415041/141‑1 (MQ=255)                                                                                        ttGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCa                                            >  1:269766/1‑141 (MQ=255)                                                                                          gCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATa                                          >  2:92373/1‑141 (MQ=255)                                                                                               aaTTGGACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTc                                     <  1:367040/141‑1 (MQ=255)                                                                                                    gACCACCATATGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATc                                <  1:263459/141‑1 (MQ=255)                                                                                                     acctccATATGCATGTGCATGGGTAAACTTATCCTAATAGCCAACTTTAGAAATGTAGCTGATAGTATAACGTGGAATACAAAACATTATAATTTGCTGTTTTGATTGATTAAAAGTAGATAAAATCATAACGTCTGATCt                               >  2:172454/1‑141 (MQ=255)                                                                                                             atGCATGTGCATGGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGtt                       <  2:391322/141‑1 (MQ=255)                                                                                                                catgtgcatgGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAg                    <  2:95925/141‑1 (MQ=255)                                                                                                                catgtgcatgGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAg                    >  1:224825/1‑141 (MQ=255)                                                                                                                  tgtgcatgGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGtt                                                                       >  1:217972/1‑88 (MQ=255)                                                                                                                  tgtgcatgGGTAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGtt                                                                       <  2:217972/88‑1 (MQ=255)                                                                                                                            tAATCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCCtt        <  2:334154/141‑1 (MQ=255)                                                                                                                              aTCTTATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCCtttt      >  2:250752/1‑141 (MQ=255)                                                                                                                                  tATCCTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCCTTTTTttc  >  2:927/1‑141 (MQ=255)                                                                                                                                      |                                                                                                                                         GCTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTGCGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCCTTTTTTTC  >  CP000730/1063993‑1064262

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