Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,187,071	C→T	R331K (AGA→AAA)	rsbU ←	sigma factor B regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,187,071	0	C	T	100.0%	107.2 / NA	31	R331K (AGA→AAA)	rsbU	sigma factor B regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/25);  total (6/25)

AAGATTCATTTCATTACACTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTCTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCATATGTTTATGTTTTTTAATATATTCTAAAAGT  >  CP000730/2186934‑2187199                                                                                                                                          |                                                                                                                                 aaGATTCATTTCATTACACTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTtttt                                                                                                                               >  2:221427/1‑141 (MQ=255)         tttcattACACTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATc                                                                                                                       <  2:325181/141‑1 (MQ=255)         tttcattACACTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATc                                                                                                                       <  1:140959/141‑1 (MQ=255)         tttcattACACTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATc                                                                                                                       <  1:350911/141‑1 (MQ=255)                acacTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAg                                                                                                                <  1:214385/141‑1 (MQ=255)                acacTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAg                                                                                                                <  1:214390/141‑1 (MQ=255)                acacTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATAACCATCGAAATAATTTATAAACTTCTTCATTAAATTAAATTACTTTTTTTATAATCAAAATAg                                                                                                                >  1:57764/1‑141 (MQ=255)                 cacTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGt                                                                                                               >  2:80473/1‑141 (MQ=255)                   cTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCa                                                                                                             <  2:184144/141‑1 (MQ=255)                       tACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATatc                                                                                                         <  1:221427/141‑1 (MQ=255)                                                 tttAAGCTCTTTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCt                                                                               <  2:57764/141‑1 (MQ=255)                                                  ttAAGCGNTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCtt                                                                              <  2:296714/141‑1 (MQ=255)                                                     aGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTaa                                                                           >  2:152950/1‑141 (MQ=255)                                                       cgcTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTaaaa                                                                         >  2:157889/1‑141 (MQ=255)                                                                caccaTATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCAt                                                                <  1:80473/141‑1 (MQ=255)                                                                              ttCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACa                                                  <  1:82153/141‑1 (MQ=255)                                                                                   attattATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATc                                             <  2:244382/141‑1 (MQ=255)                                                                                   attattATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATc                                             <  1:69233/141‑1 (MQ=255)                                                                                     tattatATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCtt                                           <  1:77604/141‑1 (MQ=255)                                                                                      attatATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTtg                                          <  1:368849/141‑1 (MQ=255)                                                                                                          atTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCAtatgtttatgtt                      <  1:351850/141‑1 (MQ=255)                                                                                                          atTTCTAATCTTCTTTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCAtatgtttatgtt                      <  1:155703/141‑1 (MQ=255)                                                                                                                   cttcttTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCATATGTTTATGTTTTTTAatat             <  2:1031/141‑1 (MQ=255)                                                                                                                   cttcttTNTTCAATTTCTTTACTTTTTTTATACTCAAAATAGTCATATCATTTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATTTATCATTTGAACATTACCTTGTGGGTGCATATGTTCATCTTTTTTAatat             <  2:295985/141‑1 (MQ=255)                                                                                                                      cttTTTTAAATTATTTTACTTTTTATATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAATCTTTAAAATTGCTTAATATATAATTTAAACAATATCTTGTGTGTGCATATGTTTATGTTTTTTAATATAtt          <  1:227700/141‑1 (MQ=255)                                                                                                                      cttTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCATATGTTTATGTTTTTTAATATAtt          <  1:167034/141‑1 (MQ=255)                                                                                                                      cttTTTTAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCATATGTTTATGTTTTTTAATATAtt          <  1:287869/141‑1 (MQ=255)                                                                                                                        tttttAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTTTTTGTTTTTTGAAGCTTTAAAATTGCTTCATATATAATTTTAACAATATCTTGTGGTTGCATATGTTTATGTTTTTTAATATAtttt        >  2:298410/1‑139 (MQ=255)                                                                                                                        tttttAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCATATGTTTATGTTTTTTAATATATTCt        <  1:150257/141‑1 (MQ=255)                                                                                                                         ttttAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCATATGTTTATGTTTTTTAATATATTCTa       <  2:222628/141‑1 (MQ=255)                                                                                                                            tAAATTAATTTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCATATGTTTATGTTTTTTAATATATTCTaaaa    <  1:253146/141‑1 (MQ=255)                                                                                                                              aattaattTACTTTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCATATGTTTATGTTTTTTAATATATTCTAAAAGt  >  1:31402/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                                 AAGATTCATTTCATTACACTCCTACTTATTAAAAATATTTATCTTCATTTTAAGCGCTGTATCCACCATATTTATATTTCAAATTATTATATACCCATCGAAATAATTTCTAATCTTCTTTTTTAAATTAATTTACTCTTTTTATAATCAAAATAGTCATATCATCTTTTTTATTTGGGTTTTGAAGCTTTAAAATTGCTTCATAGATAATTTGAACAATATCTTGTGGGTGCATATGTTTATGTTTTTTAATATATTCTAAAAGT  >  CP000730/2186934‑2187199

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