Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,440,950	G→A	T27I (ACC→ATC)	rpoA ←	RNA polymerase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,440,950	0	G	A	100.0%	186.9 / NA	59	T27I (ACC→ATC)	rpoA	RNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/32);  total (27/32)

GTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGT  >  NC_000913/3440826‑3441083                                                                                                                             |                                                                                                                                      gTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcgtgc                                                                                                                         <  2:972712/139‑1 (MQ=255)  tACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcgtgcg                                                                                                                        <  1:1234608/139‑1 (MQ=255)         gtagtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTTCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACt                                                                                                                 <  1:799618/139‑1 (MQ=255)               cacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTg                                                                                                           >  2:244966/1‑139 (MQ=255)                   aTCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCg                                                                                                       >  1:786915/1‑139 (MQ=255)                   aTCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCg                                                                                                       >  1:774746/1‑139 (MQ=255)                    tCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGCCGAACTCACTTGCTCGa                                                                                                      >  1:1264797/1‑139 (MQ=255)                     cAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                                     >  1:1327966/1‑139 (MQ=255)                       aTCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                                                   >  2:776614/1‑139 (MQ=255)                                 cGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         <  1:223999/139‑1 (MQ=255)                                   gTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              >  1:1153889/1‑116 (MQ=255)                                   gTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              <  2:1153889/116‑1 (MQ=255)                                             aaCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGaaa                                                                             <  2:1184453/139‑1 (MQ=255)                                                ccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActc                                                                          <  1:776614/139‑1 (MQ=255)                                                       cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  1:244966/139‑1 (MQ=255)                                                       cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  2:786915/139‑1 (MQ=255)                                                       cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  2:611921/139‑1 (MQ=255)                                                       cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   >  2:479138/1‑139 (MQ=255)                                                         aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGa                                                                 >  1:475901/1‑139 (MQ=255)                                                                     tACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtgt                                                     <  1:479138/139‑1 (MQ=255)                                                                          cgcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  1:22402/82‑1 (MQ=255)                                                                          cgcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  2:22402/1‑82 (MQ=255)                                                                          cgcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                <  2:885401/139‑1 (MQ=255)                                                                             aGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   >  2:1219477/1‑117 (MQ=255)                                                                             aGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  1:1219477/117‑1 (MQ=255)                                                                               tGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCtt                                           <  2:774746/139‑1 (MQ=255)                                                                               tGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCATCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCtt                                           <  2:475901/139‑1 (MQ=255)                                                                                         agagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         >  2:997412/1‑83 (MQ=255)                                                                                         agagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         <  1:997412/83‑1 (MQ=255)                                                                                           agTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                <  1:1305455/122‑1 (MQ=255)                                                                                           agTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                >  2:1305455/1‑122 (MQ=255)                                                                                               tGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTAc                           >  2:1259084/1‑139 (MQ=255)                                                                                                ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTg                                                                        <  1:331870/93‑1 (MQ=255)                                                                                                ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTg                                                                        >  2:331870/1‑93 (MQ=255)                                                                                                     aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         <  2:828152/71‑1 (MQ=255)                                                                                                     aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         >  1:828152/1‑71 (MQ=255)                                                                                                     aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTa                                                                                 >  1:588426/1‑79 (MQ=255)                                                                                                     aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTa                                                                                 <  2:588426/79‑1 (MQ=255)                                                                                                     aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAg                     <  2:1133776/139‑1 (MQ=255)                                                                                                      aGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcgtgcg                                                                                                                        <  1:720836/39‑1 (MQ=255)                                                                                                      aGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcgtgcg                                                                                                                        >  2:720836/1‑39 (MQ=255)                                                                                                       gCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGaa                                                                                                                   <  1:655675/43‑1 (MQ=255)                                                                                                       gCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGaa                                                                                                                   >  2:655675/1‑43 (MQ=255)                                                                                                       gCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         <  2:98756/69‑1 (MQ=255)                                                                                                       gCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         >  1:98756/1‑69 (MQ=255)                                                                                                              tctAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcg            <  2:1264797/139‑1 (MQ=255)                                                                                                                       tCAAGGATCACCTTTGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                                       <  2:833320/55‑1 (MQ=255)                                                                                                                       tCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                                       >  1:833320/1‑55 (MQ=255)                                                                                                                        cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  2:996702/36‑1 (MQ=255)                                                                                                                        cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  1:996702/1‑36 (MQ=255)                                                                                                                        cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                                  <  2:29572/91‑1 (MQ=255)                                                                                                                        cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                                  >  1:29572/1‑91 (MQ=255)                                                                                                                        cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTAc                                     <  2:213810/104‑1 (MQ=255)                                                                                                                        cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTAc                                     >  1:213810/1‑104 (MQ=255)                                                                                                                        cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACtt                         >  1:1250212/1‑116 (MQ=255)                                                                                                                        cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACtt                         <  2:1250212/116‑1 (MQ=255)                                                                                                                        cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgt  >  1:755874/1‑139 (MQ=255)                                                                                                                          aGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAtt                                                       >  1:324734/1‑84 (MQ=255)                                                                                                                          aGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAtt                                                       <  2:324734/84‑1 (MQ=255)                                                                                                                             |                                                                                                                                      GTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGT  >  NC_000913/3440826‑3441083

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