Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,184,514	C→T	N1090N (AAC→AAT)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,184,514	0	C	T	100.0%	175.1 / NA	48	N1090N (AAC→AAT)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/22);  total (26/22)

GGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAG  >  NC_000913/4184370‑4184656                                                                                                                                                 |                                                                                                                                               ggCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTa                                                                                                                                            <  2:583240/149‑1 (MQ=255)  gCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGc                                                                                                                                         <  1:248340/151‑1 (MQ=255)          cGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGAc                                                                                                                                 <  2:389845/151‑1 (MQ=255)              cTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATc                                                                                                                              >  1:32037/1‑150 (MQ=255)               tGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGt                                                                                                                            >  2:610645/1‑151 (MQ=255)               tGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGt                                                                                                                            >  1:836711/1‑151 (MQ=255)               tGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGt                                                                                                                            >  1:644482/1‑151 (MQ=255)                       gTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAAcc                                                                                                                    <  2:46538/151‑1 (MQ=255)                          aaGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGc                                                                                                                 >  2:78170/1‑151 (MQ=255)                          aaGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGc                                                                                                                 >  1:461722/1‑151 (MQ=255)                                     ggCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCg                                                                                                      >  1:568512/1‑151 (MQ=255)                                            aaaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGt                                                                                                >  1:544015/1‑150 (MQ=255)                                             aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTAt                                                                                              >  2:786161/1‑151 (MQ=255)                                                    tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATc                                                                                       >  1:579563/1‑151 (MQ=255)                                                     aTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATc                                                                                       <  1:533004/150‑1 (MQ=255)                                                         aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                                  <  1:78147/151‑1 (MQ=255)                                                            cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                   >  2:135300/1‑147 (MQ=255)                                                            cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                   >  1:160449/1‑147 (MQ=255)                                                            cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                   >  2:814735/1‑147 (MQ=255)                                                             tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                   >  2:301708/1‑146 (MQ=255)                                                              ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                               <  2:644524/149‑1 (MQ=255)                                                                 gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                           <  1:469175/150‑1 (MQ=255)                                                                    aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                       <  2:411826/151‑1 (MQ=255)                                                                       aTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAccc                                                                    >  2:493766/1‑151 (MQ=255)                                                                             ggTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                              <  2:58429/151‑1 (MQ=255)                                                                                         aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCg                                                   >  2:437005/1‑150 (MQ=255)                                                                                           cAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                 >  1:10160/1‑150 (MQ=255)                                                                                            aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                               <  1:205076/151‑1 (MQ=255)                                                                                              gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGt                                             <  2:85373/151‑1 (MQ=255)                                                                                                gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAt                                           >  1:146903/1‑151 (MQ=255)                                                                                                    aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                        <  1:508182/150‑1 (MQ=255)                                                                                                     aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCg                                      <  1:363094/151‑1 (MQ=255)                                                                                                        tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATACTACGCCGGTAGACATCGTACTGAACCCGCTAGGCGTACCGTCACGTATGAACATCGGTCAGATCCTCGAAACCCACCAGGGTATGGCTGCGGAAGGAAACGGAGAc                                    >  1:703763/1‑150 (MQ=255)                                                                                                          tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACa                                   >  2:567086/1‑149 (MQ=255)                                                                                                           aaGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACa                                   >  2:352002/1‑148 (MQ=255)                                                                                                           aaGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                <  2:365035/151‑1 (MQ=255)                                                                                                            aGATCAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                               <  2:414106/151‑1 (MQ=255)                                                                                                                cAACCCGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                           >  2:39030/1‑151 (MQ=255)                                                                                                                    ccGATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATGGTGCTGATCCCGCTGGGCGGACCGTCTCGTATGAACATCGGTCAGATCCTCGAATCCCACCTGGGTATGGCTGCGAAATGTAGCGGCGAGAAGATCAACggaa                       <  1:135167/151‑4 (MQ=255)                                                                                                                      gATCGAAGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATg                     >  2:737288/1‑151 (MQ=255)                                                                                                                           aaGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaa                <  2:182828/151‑1 (MQ=255)                                                                                                                           aaGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGa                 >  2:596369/1‑150 (MQ=255)                                                                                                                            aGATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCTAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa               <  1:833841/151‑1 (MQ=255)                                                                                                                             gATATGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAc              <  1:590936/151‑1 (MQ=255)                                                                                                                                atGCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag            <  1:23661/150‑1 (MQ=255)                                                                                                                                  gCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc           >  2:815868/1‑149 (MQ=255)                                                                                                                                  gCCTTACGATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc           >  1:203031/1‑149 (MQ=255)                                                                                                                                         gATGAAAATGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaag  <  2:693265/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                               GGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAG  >  NC_000913/4184370‑4184656

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