Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,184,542	C→A	P1100Q (CCG→CAG)	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,542	0	C	A	100.0%	144.5 / NA	53	P1100Q (CCG→CAG)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (29/24);  total (29/24)

GGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184408‑4184677                                                                                                                                       |                                                                                                                                        ggTTAAACTCCGGGTCTAGACTGGAGAAAAGAGAGCAGGTCGTCACGGTAACAAGGGTGTAATTTATAAGATAAAACCGATCGAAGATATGGAGTAAGAGGAAAATGGTCAGACGGTAGACATAGTACTGAAACAGCTgggag                                                                                                                                 >  2:6733/1‑141 (MQ=255)      aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCg                                                                                                                            >  2:38875/1‑143 (MQ=255)            gTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGt                                                                                                                      <  2:119638/143‑1 (MQ=255)              attcAGCGTGGTGACAAGAGGGCAGGGCGGCACGGTGAGAAGGGTGTAATTTGTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAATACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                    <  1:634880/140‑1 (MQ=255)              aTCCAGCCTGGTGACAAGAAGGCAGGTCGTCACGGGAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGAGATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                    <  2:19997/143‑1 (MQ=255)                ccAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGCTATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCCCGTagtg                                                                                                                  >  2:470724/1‑140 (MQ=255)                ccAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGGACCGTCTCGTATGa                                                                                                                  >  2:114392/1‑143 (MQ=255)                 cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGaa                                                                                                                 >  2:505204/1‑143 (MQ=255)                          gACGAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCGTACGATGAAAACGGTACGCCGGTAGACATCGTACCGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                                                        <  2:275805/143‑1 (MQ=255)                          gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                                                        >  2:44881/1‑143 (MQ=255)                           aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAACACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGGATGAACACCGGCCAg                                                                                                       >  2:304035/1‑143 (MQ=255)                                attgCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCt                                                                                                  <  1:128981/140‑1 (MQ=255)                                 tgtcAGGTCGTCACGGTAACAAGGTTGTAATTGCTAAGATCAACCCGATCGAGGATATGCCTTACGATGAAAGCGGTACGCCGGTAGATATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                 <  2:145529/140‑1 (MQ=255)                                     aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTCTGAACATCGGGCAGCTCCTCGaaa                                                                                             >  1:8720/1‑143 (MQ=255)                                          gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                        >  2:311720/1‑143 (MQ=255)                                           tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCCAAACCCAcc                                                                                       >  1:307131/1‑143 (MQ=255)                                           tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGCACATCGGTCAGATCCTCGAAGCCCAcc                                                                                       >  2:454305/1‑143 (MQ=255)                                           tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                                                       >  1:522138/1‑143 (MQ=255)                                           tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                                                       >  2:515038/1‑143 (MQ=255)                                             aCGGCAACAAGGGTGTAACTGCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTg                                                                                     <  1:439520/143‑1 (MQ=255)                                                  aaCAGGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                                                                <  2:485892/143‑1 (MQ=255)                                                           gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGCAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTAAGATCCTCGACACCCACCTGGGTATGGCTGCGaa                                                                       >  2:362556/1‑143 (MQ=255)                                                            tAATGTCTAAGATCACCCCGCTCGAAGATATGCCTTACGAGGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                      <  1:672039/143‑1 (MQ=255)                                                             aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                     <  1:428204/143‑1 (MQ=255)                                                                    aaGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                              >  1:729595/1‑143 (MQ=255)                                                                      gATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCATCGAAACCCACCTGGGTATGGCTGCGAAAGGTAACGCCg                                                            >  2:102173/1‑143 (MQ=255)                                                                            cccgaccgcctctctGCCTTAGGATGGCACCGTTACGCCGGTACATAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGACCACCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                      <  2:146631/129‑1 (MQ=255)                                                                            cccGATCGAAGATATGCCTTACGATGAAAGCGGTACGCCGGTCGACCTCGTACTGACCCAGCAGGGCGTACCGTCCGGTATGAACATCGGGCAGATGCTCGAAACCAACCTGGGTATGGCTGCGGAAGGTTGCGGCGagacgc                                                      >  1:478086/1‑138 (MQ=255)                                                                            cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGCCCAAGCTGGGCGGACCCTCTCGTACGAACAGCGCTCAGATACCCGCAACCCCCCTGGGCGTGGCGGCGACAgggcgcggcgccacta                                                      >  2:363566/1‑129 (MQ=255)                                                                            cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACGGAACCAGCTGGGCGTACCGGCTCGTATGAACGTCGGGCAGATCCTCGAAACCAACCGCGGTATGGCTGCGAAAGGAATCGGCGCCAAGa                                                      >  1:158081/1‑143 (MQ=255)                                                                              cGATCGAAGGCGTGCCGTACGATGAAAACGGTACGCCGGCAGGCATCGTACTTAACCAGCTGGGCGTACCGTCTGGGATGAACATCGGTCAGAGCCTCGAAACCCACCTGGGTATGTCTGCGAAAGGTATCGGCGACAAGATc                                                    <  2:696370/143‑1 (MQ=255)                                                                                 tCGAAGATATGCCGTACGATGAAAACGGTACGCCGGTAGGCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                 <  2:45122/143‑1 (MQ=255)                                                                                   gAAGCTATGCCTTTCGATGAAAGCGGTACGCCAGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGc                                               <  2:686170/143‑1 (MQ=255)                                                                                    aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACGAGATCAACGcc                                              >  1:54358/1‑143 (MQ=255)                                                                                    aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                              >  1:617703/1‑143 (MQ=255)                                                                                         atGCCCTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCt                                         <  1:508075/143‑1 (MQ=255)                                                                                                aCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCCACGCCATGCTGacaccg                                  >  1:439065/1‑141 (MQ=255)                                                                                                  gacgAAAACGGGGCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTAGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                <  1:450050/140‑1 (MQ=255)                                                                                                  gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                <  1:38875/143‑1 (MQ=255)                                                                                                   aTGAAATCGTCACGCCGGTAGACCTCGTACTGACCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagcag                               <  2:248177/143‑1 (MQ=255)                                                                                                   aTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagcag                               >  2:595608/1‑143 (MQ=255)                                                                                                   aTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagcag                               >  2:622895/1‑143 (MQ=255)                                                                                                   aTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCGGCGAAAGGTATCGGCGACAAGATCACCGCCATGCTGAAAcagccg                               >  1:318050/1‑141 (MQ=255)                                                                                                       aaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACCAGATCAACGCCATGCTGAAACAGCAGCaag                           >  1:545852/1‑143 (MQ=255)                                                                                                           ggTACGCCGGTAGTCATCGCACTGAACCAGCTGGGCGTACCGTCTCGTAGGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGt                       <  1:701058/143‑1 (MQ=255)                                                                                                             tACGCCGGTAGACATCGTACTGAACCAGCGGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcg                     <  1:404903/143‑1 (MQ=255)                                                                                                              aCGCCGGTGGACGTCGTACGGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTGGAGCCCCACCGGCGGATGGATGCGGTAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                    <  2:636520/143‑1 (MQ=255)                                                                                                                      aggaaTCGTGCTGTATCAGCTGTGCGTACCGTCTCGTATGAACAGCGGTCAGATCCTCGAACCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgc            <  2:78512/139‑1 (MQ=255)                                                                                                                       gACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGGCGAGATCAACGCCATGCTGCAACAGCAGCAAGGAGTCGCGAAACTgcg           >  1:685597/1‑143 (MQ=255)                                                                                                                           tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCAGAAAAAGAAGCAAGAAGTCACGAAAATGCGCGaa       >  1:98625/1‑143 (MQ=255)                                                                                                                              tccTGCACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTc    <  1:645899/141‑1 (MQ=255)                                                                                                                                cTGGACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAt  <  1:346300/143‑1 (MQ=255)                                                                                                                                cTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACAGCGCGAATTCAt  >  2:8274/1‑143 (MQ=255)                                                                                                                                       |                                                                                                                                        GGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184408‑4184677

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