Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,186,152	T→G	L268R (CTG→CGG)	rpoC →	RNA polymerase, beta prime subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,186,152	0	T	G	100.0%	205.5 / NA	54	L268R (CTG→CGG)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (28/26);  total (28/26)

AGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATC  >  NC_000913/4186005‑4186297                                                                                                                                                    |                                                                                                                                                  aGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTa                                                                                                                                                <  1:99005/151‑1 (MQ=255)  gCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTa                                                                                                                                                <  1:831617/150‑1 (MQ=255)      aTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTAtcgtc                                                                                                                                           <  1:425723/151‑1 (MQ=255)      aTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTAtcgt                                                                                                                                            <  2:254335/150‑1 (MQ=255)        cAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTAtcgtcg                                                                                                                                          <  2:667020/150‑1 (MQ=255)         aaaCTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTcgcg                                                                                                                                        <  2:439106/151‑1 (MQ=255)           aCTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTc                                                                                                                                      >  2:171864/1‑151 (MQ=255)             tgctgGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCAt                                                                                                                                    >  2:381850/1‑151 (MQ=255)                 ggAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAAc                                                                                                                                >  1:548086/1‑151 (MQ=255)                  gAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAAcc                                                                                                                               >  2:121676/1‑151 (MQ=255)                     gCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTa                                                                                                                            <  2:233596/151‑1 (MQ=255)                       gttcgttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTaac                                                                                                                          <  1:27955/151‑1 (MQ=255)                          cgttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTaacaac                                                                                                                       <  1:101367/151‑1 (MQ=255)                           gttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAAcc                                                                                                                      >  2:269606/1‑151 (MQ=255)                             tcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGt                                                                                                                    >  1:590128/1‑151 (MQ=255)                               aGTGTGGTAAAAAAGCAGAGTGGATGTTCGTGCCCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTATCGCTGCTTCTGACCTGAACGATCGGTATCGTCGCGTCATTGACAGTagcggccg                                                                                                                     <  1:17553/148‑8 (MQ=255)                                gTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTg                                                                                                                 >  2:685735/1‑151 (MQ=255)                                          aaaCCAGAGGGGATGATCCTGACCGTTCTGCCCGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTAGCGTCATTAACCGTAACAACCGTCTGAAACGTctgc                                                                                                       >  1:593890/1‑151 (MQ=255)                                             ccAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTgg                                                                                                    <  1:797372/151‑1 (MQ=255)                                              cAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTgg                                                                                                    >  2:796064/1‑150 (MQ=255)                                                  gTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTg                                                                                               <  2:351233/151‑1 (MQ=255)                                                     gatgatCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCt                                                                                            >  2:850710/1‑151 (MQ=255)                                                     gatgatCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCt                                                                                            >  2:111853/1‑151 (MQ=255)                                                       tgatCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTgc                                                                                          >  1:490884/1‑151 (MQ=255)                                                         atCCTGCCCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTgcgc                                                                                        >  1:176455/1‑151 (MQ=255)                                                          tCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGcc                                                                                       >  2:49065/1‑151 (MQ=255)                                                                   ttCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGAcatcat                                                                              >  2:630417/1‑151 (MQ=255)                                                                       gCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcg                                                                            <  1:664979/149‑1 (MQ=255)                                                                        ccGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcgtac                                                                         <  1:440444/151‑1 (MQ=255)                                                                        ccGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcgtac                                                                         >  1:672448/1‑151 (MQ=255)                                                                             aCTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTaa                                                                     <  1:418614/150‑1 (MQ=255)                                                                                gccgccAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACg                                                                   <  2:488417/149‑1 (MQ=255)                                                                                 ccgccAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGaa                                                                 >  2:84941/1‑150 (MQ=255)                                                                                   gccAGATCTGCGTCCGCTGGTTCCGCTGGATGGTTGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGaaaaa                                                              <  1:469208/151‑1 (MQ=255)                                                                                          cTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTAtgc                                                       >  1:669870/1‑151 (MQ=255)                                                                                            gCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTAtgctg                                                     <  2:805175/151‑1 (MQ=255)                                                                                            gCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTAtgc                                                       <  2:590193/149‑1 (MQ=255)                                                                                              gTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTAtgctgc                                                    >  2:122944/1‑150 (MQ=255)                                                                                                ccGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAgg                                                 <  2:828170/151‑1 (MQ=255)                                                                                                     ggTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCg                                            <  2:232044/151‑1 (MQ=255)                                                                                                         ccGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGt                                          >  2:327927/1‑149 (MQ=255)                                                                                                         ccGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTa                                         <  2:272086/150‑1 (MQ=255)                                                                                                          cGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGa                                       <  2:181699/151‑1 (MQ=255)                                                                                                                       tCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGAt                          >  2:517767/1‑151 (MQ=255)                                                                                                                        cGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGAt                          >  1:249250/1‑150 (MQ=255)                                                                                                                        cGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGAt                          >  2:819534/1‑150 (MQ=255)                                                                                                                         gTTTCGCGACTTCTGACCTGAACGATCGGCATCGTCGCGTCATTAACCGTTACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTAAGTAACGAAAAACGTATGCTGCAGGAAGTGGTAGACGCCCTGCTGGAt                          >  1:867047/1‑149 (MQ=255)                                                                                                                              gcgACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtc                   >  2:188008/1‑151 (MQ=255)                                                                                                                              gcgACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtc                   <  2:132593/151‑1 (MQ=255)                                                                                                                              gcgACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtc                   <  2:94866/151‑1 (MQ=255)                                                                                                                                gACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtcg                  <  1:243178/150‑1 (MQ=255)                                                                                                                                 aCTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtcgtc                >  1:229532/1‑151 (MQ=255)                                                                                                                                      tGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGt           >  2:38810/1‑151 (MQ=255)                                                                                                                                               cGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATc  >  1:663172/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  AGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATC  >  NC_000913/4186005‑4186297

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑