Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,185,540	C→T	P64L (CCG→CTG)	rpoC →	RNA polymerase, beta prime subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,185,540	0	C	T	100.0%	242.9 / NA	69	P64L (CCG→CTG)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/34);  total (35/34)

ATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGC  >  NC_000913/4185407‑4185672                                                                                                                                      |                                                                                                                                     aTCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTaa                                                                                                                                 <  2:1226962/139‑1 (MQ=255)           ctctGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGa                                                                                                                       >  1:1357416/1‑139 (MQ=255)            tctGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAg                                                                                                                      <  2:1174251/139‑1 (MQ=255)             ctGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGt                                                                                                                     >  1:45287/1‑139 (MQ=255)                gCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGcc                                                                                                                  >  2:307883/1‑139 (MQ=255)                   tCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCtgt                                                                                                               <  2:472553/139‑1 (MQ=255)                   tCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCtgt                                                                                                               <  1:1680276/139‑1 (MQ=255)                   tCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCtgt                                                                                                               <  2:1741121/139‑1 (MQ=255)                     gCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGcc                                                                                                                  >  1:360254/1‑134 (MQ=255)                     gCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGcc                                                                                                                  <  2:360254/134‑1 (MQ=255)                      ccAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCg                                                                                                            <  2:45287/139‑1 (MQ=255)                          aCATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAg                                                                                                                              >  1:1645684/1‑117 (MQ=255)                          aCATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAg                                                                                                                              <  2:1645684/117‑1 (MQ=255)                          aCATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTAc                                                                                                                         <  1:981609/122‑1 (MQ=255)                          aCATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTAc                                                                                                                         >  2:981609/1‑122 (MQ=255)                          aCATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCtgt                                                                                                               >  1:1402219/1‑132 (MQ=255)                          aCATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCtgt                                                                                                               <  2:1402219/132‑1 (MQ=255)                          aCATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTaa                                                                                                        >  1:1026147/1‑139 (MQ=255)                               aTCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAg                                                                                                                              <  1:698739/112‑1 (MQ=255)                               aTCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAg                                                                                                                              >  2:698739/1‑112 (MQ=255)                                  cGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAgc                                                                                                <  2:206619/139‑1 (MQ=255)                                    ttCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACaa                                                                                                  >  1:403106/1‑135 (MQ=255)                                    ttCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACaa                                                                                                  <  2:403106/135‑1 (MQ=255)                                    ttCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAgcgc                                                                                              >  1:33240/1‑139 (MQ=255)                                         ggTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGaa                                                                                         >  1:772718/1‑139 (MQ=255)                                          gTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCt                                                                                                                 <  1:1255842/114‑1 (MQ=255)                                          gTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCt                                                                                                                 >  2:1255842/1‑114 (MQ=255)                                                     aaGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGt                                                                             <  2:33240/139‑1 (MQ=255)                                                          aaaaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGa                                                                                                                             >  1:89049/1‑86 (MQ=255)                                                          aaaaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGa                                                                                                                             <  2:89049/86‑1 (MQ=255)                                                            aaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGcc                                                                                                                  <  2:1499220/95‑1 (MQ=255)                                                            aaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGcc                                                                                                                  >  1:1499220/1‑95 (MQ=255)                                                                 cGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGaa                                                                 <  2:1575821/139‑1 (MQ=255)                                                                 cGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGaa                                                                 >  1:1575821/1‑139 (MQ=255)                                                                  ggAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGt                                                                                                          <  2:1629416/97‑1 (MQ=255)                                                                  ggAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGt                                                                                                          >  1:1629416/1‑97 (MQ=255)                                                                        cATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCtgt                                                                                                               <  1:804750/86‑1 (MQ=255)                                                                        cATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCtgt                                                                                                               >  2:804750/1‑86 (MQ=255)                                                                         aTCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTgcggcg                                                         <  2:1490475/139‑1 (MQ=255)                                                                         aTCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTgcggcg                                                         <  1:307883/139‑1 (MQ=255)                                                                          tCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGt                                                        >  1:1707834/1‑139 (MQ=255)                                                                                  cGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTaa                                                                                                                                 <  2:1678359/58‑1 (MQ=255)                                                                                  cGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTaa                                                                                                                                 >  1:1678359/1‑58 (MQ=255)                                                                                  cGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGa                                                >  1:548950/1‑139 (MQ=255)                                                                                    tACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGa                                                                                                                             >  1:1415169/1‑60 (MQ=255)                                                                                    tACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGa                                                                                                                             <  2:1415169/60‑1 (MQ=255)                                                                                    tACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAgcgc                                                                                              >  1:1397545/1‑91 (MQ=255)                                                                                    tACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAgcgc                                                                                              <  2:1397545/91‑1 (MQ=255)                                                                                         tCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTaa                                                                                                                                 <  1:984309/51‑1 (MQ=255)                                                                                         tCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTaa                                                                                                                                 >  2:984309/1‑51 (MQ=255)                                                                                         tCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAATGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCg                                                                              <  2:1337702/102‑1 (MQ=255)                                                                                         tCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCg                                                                              >  1:1337702/1‑102 (MQ=255)                                                                                             aCCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTaaa                                     >  2:1708873/1‑139 (MQ=255)                                                                                                  aaCGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACg                                <  2:772718/139‑1 (MQ=255)                                                                                                       gACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCg                                                                                                            <  1:1145147/58‑1 (MQ=255)                                                                                                       gACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCg                                                                                                            >  2:1145147/1‑58 (MQ=255)                                                                                                       gACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTg                                                              <  1:970912/104‑1 (MQ=255)                                                                                                       gACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTg                                                              >  2:970912/1‑104 (MQ=255)                                                                                                       gACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACt                                        <  1:303726/126‑1 (MQ=255)                                                                                                       gACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACt                                        >  2:303726/1‑126 (MQ=255)                                                                                                       gACGGCCTTTTCTGCGCCCGTATCTTTGGGCTGATAAAAGATTACGAGTGCCCGTGCGGTAAGTACAAGCGCCTGAAACCCCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTg                           >  1:845811/1‑139 (MQ=255)                                                                                                         cGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGAcc                                              <  1:893/118‑1 (MQ=255)                                                                                                         cGGCCTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGAcc                                              >  2:893/1‑118 (MQ=255)                                                                                                            ccTTTTCTGCGCCCGTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACAGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGt                      >  1:1281712/1‑139 (MQ=255)                                                                                                                          gTATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGa        >  1:379384/1‑139 (MQ=255)                                                                                                                           tATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTgcggcg                                                         >  1:14124/1‑89 (MQ=255)                                                                                                                           tATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTgcggcg                                                         <  2:14124/89‑1 (MQ=255)                                                                                                                           tATCTTTGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGaa       >  1:933851/1‑139 (MQ=255)                                                                                                                                ttGGGCTGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGc  <  2:1357416/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     ATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGC  >  NC_000913/4185407‑4185672

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