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%	147.0 / NA	46	T27I (ACC→ATC)	rpoA	RNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/22);  total (24/22)

ACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGT  >  NC_000913/3440839‑3441083                                                                                                                |                                                                                                                                      acacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGaa                                                                                                                   >  2:249995/1‑132 (MQ=255) acacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGaa                                                                                                                   <  1:249995/132‑1 (MQ=255)   acCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGc                                                                                                          <  2:961365/139‑1 (MQ=255)    ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  2:558947/1‑139 (MQ=255)    ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  1:740742/1‑139 (MQ=255)    ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGCTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  1:180604/139‑1 (MQ=255)    ccATCAATCTCAACCTCGGTCACCACGCAACCCGGCATAGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGCGTATGGCCAAAGCCACGCTTTAAAGGCTCAAGGATCCCCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  2:983875/139‑1 (MQ=255)     cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTc                                                                                                        >  1:281670/1‑139 (MQ=255)     cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTc                                                                                                        >  2:870170/1‑139 (MQ=255)      aTCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCg                                                                                                       >  2:936851/1‑139 (MQ=255)          aTCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                                                   >  1:508167/1‑139 (MQ=255)               aaCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc                                                                                              >  1:567621/1‑139 (MQ=255)                   tCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgc                                                                                          <  2:664844/139‑1 (MQ=255)                    cGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                                     >  1:613407/1‑127 (MQ=255)                    cGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTAACTTGcgccat                                                                                                     <  2:613407/127‑6 (MQ=255)                     ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcgc                                                                                        >  1:568917/1‑139 (MQ=255)                     ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAAGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcgc                                                                                        >  2:37870/1‑139 (MQ=255)                      gTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                                       <  2:568917/139‑1 (MQ=255)                                 aCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGGTCTAAAGGCTCAAGGATCACCTTGGCGTGCGACGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAAc                                                                            >  2:686003/1‑139 (MQ=255)                                  cccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCACTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAAct                                                                           >  2:688138/1‑139 (MQ=255)                                   ccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAGActc                                                                          <  1:530767/139‑1 (MQ=255)                                          cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  1:870170/139‑1 (MQ=255)                                            aTGAGGGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAACGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGa                                                                 >  1:2609/1‑139 (MQ=255)                                            aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActc                                                                          <  2:404737/130‑1 (MQ=255)                                            aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActc                                                                          >  1:404737/1‑130 (MQ=255)                                                    aGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  2:838405/91‑1 (MQ=255)                                                    aGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  1:838405/1‑91 (MQ=255)                                                    aGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCa                                                         <  1:545026/139‑1 (MQ=255)                                                      aaTACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAtt                                                       <  1:936851/139‑1 (MQ=255)                                                                aGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctctc                                             >  1:205588/1‑139 (MQ=255)                                                                           cAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTaa                                  <  2:281670/139‑1 (MQ=255)                                                                            agagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAg                                 >  1:518048/1‑139 (MQ=255)                                                                              agTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActc                                                                          <  2:1035145/96‑1 (MQ=255)                                                                              agTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActc                                                                          >  1:1035145/1‑96 (MQ=255)                                                                              agTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCt                               >  2:4757/1‑139 (MQ=255)                                                                                   ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCt                               >  2:140151/1‑134 (MQ=255)                                                                                   ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCt                               <  1:140151/134‑1 (MQ=255)                                                                                        aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                <  2:235575/112‑1 (MQ=255)                                                                                        aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                >  1:235575/1‑112 (MQ=255)                                                                                          gCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTa                   <  2:508167/139‑1 (MQ=255)                                                                                               gCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCt              <  2:567621/139‑1 (MQ=255)                                                                                                         cTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                >  1:906156/1‑95 (MQ=255)                                                                                                         cTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                <  2:906156/95‑1 (MQ=255)                                                                                                           cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgt  <  2:205588/139‑1 (MQ=255)                                                                                                           cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgt  <  1:81726/139‑1 (MQ=255)                                                                                                           cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgt  >  1:131888/1‑139 (MQ=255)                                                                                                                |                                                                                                                                      ACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGT  >  NC_000913/3440839‑3441083

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