Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,183,097	A→T	Q618L (CAG→CTG)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,183,097	0	A	T	100.0%	188.7 / NA	47	Q618L (CAG→CTG)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (27/20);  total (27/20)

GTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCG  >  NC_000913/4182953‑4183239                                                                                                                                                 |                                                                                                                                               gTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGa                                                                                                                                           >  2:452813/1‑150 (MQ=255)     gATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTc                                                                                                                                       <  1:581224/150‑1 (MQ=255)      aTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTcc                                                                                                                                      >  2:698123/1‑150 (MQ=255)        cAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAAc                                                                                                                                   <  1:776149/151‑1 (MQ=255)             ctctGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTg                                                                                                                                <  1:4553/149‑1 (MQ=255)                         aCGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGcc                                                                                                                   <  1:215645/150‑1 (MQ=255)                            cacaGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACtt                                                                                                               >  2:837647/1‑151 (MQ=255)                             acaGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTc                                                                                                              <  1:1035678/151‑1 (MQ=255)                                    aaCGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTagaag                                                                                                       <  2:911570/151‑1 (MQ=255)                                        aaTACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCt                                                                                                   <  1:332480/151‑1 (MQ=255)                                             ggCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTa                                                                                               >  1:1014411/1‑150 (MQ=255)                                               cTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACt                                                                                            >  1:172373/1‑151 (MQ=255)                                                 tCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACtt                                                                                           <  2:942802/150‑1 (MQ=255)                                                 tCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACtt                                                                                           <  2:997149/150‑1 (MQ=255)                                                      gagaCTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTa                                                                                     <  2:117527/151‑1 (MQ=255)                                                             cGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAgg                                                                              <  2:422321/151‑1 (MQ=255)                                                               tATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCg                                                                            >  2:487294/1‑151 (MQ=255)                                                                 tCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGa                                                                           >  1:973343/1‑150 (MQ=255)                                                                 tCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGa                                                                           >  2:343325/1‑150 (MQ=255)                                                                  cGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAAt                                                                         >  2:671952/1‑151 (MQ=255)                                                                     aaaGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATcc                                                                       >  2:502146/1‑150 (MQ=255)                                                                     aaaGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCa                                                                      <  2:139549/151‑1 (MQ=255)                                                                             cGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTc                                                              >  2:831904/1‑151 (MQ=255)                                                                                   tgttgtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCcgc                                                        >  2:867654/1‑151 (MQ=255)                                                                                         aaCTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGAcc                                                    >  1:464505/1‑149 (MQ=255)                                                                                         aaCTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAg                                                  <  2:389154/151‑1 (MQ=255)                                                                                           cTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGt                                                >  1:748170/1‑151 (MQ=255)                                                                                             gACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGt                                                <  2:567439/149‑1 (MQ=255)                                                                                                      cACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATgg                                     >  2:775129/1‑151 (MQ=255)                                                                                                      cACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATg                                      <  2:107561/150‑1 (MQ=255)                                                                                                        cTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGAc                                   >  1:1005485/1‑151 (MQ=255)                                                                                                              gTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATcc                             >  2:793863/1‑151 (MQ=255)                                                                                                                 tGCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCAcc                          >  2:25029/1‑151 (MQ=255)                                                                                                                  gCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCa                            <  2:445927/148‑1 (MQ=255)                                                                                                                  gCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCAccc                         >  1:256432/1‑151 (MQ=255)                                                                                                                  gCTATCGAAGAAGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCAcc                          <  1:976380/150‑1 (MQ=255)                                                                                                                     aTCGAAGAAGGCAACTACGTTAGCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTTGAACGCCTGGTAACTTGCCGTAGCAAAGGGGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCcagc                      <  1:684603/151‑1 (MQ=255)                                                                                                                        gaagaaGGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAgg                   >  2:1023495/1‑151 (MQ=255)                                                                                                                         aagaagGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAggt                  >  2:407744/1‑151 (MQ=255)                                                                                                                         aagaagGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAgg                   >  2:894434/1‑150 (MQ=255)                                                                                                                          agaagGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAggtg                 >  1:755961/1‑151 (MQ=255)                                                                                                                            aagGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAggtggt               >  1:100412/1‑151 (MQ=255)                                                                                                                             agGCAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAggtggt               >  1:406692/1‑150 (MQ=255)                                                                                                                                cAACTACGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATcc           <  1:61627/151‑1 (MQ=255)                                                                                                                                     aCGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCg       <  2:717984/150‑1 (MQ=255)                                                                                                                                      cGTTATCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGt     >  2:773390/1‑151 (MQ=255)                                                                                                                                          aTCGCCCTGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCg  >  2:454997/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                               GTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCG  >  NC_000913/4182953‑4183239

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