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

GTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGT  >  NC_000913/3440837‑3441083                                                                                                                  |                                                                                                                                      gtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              >  2:1600542/1‑139 (MQ=255)     acCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGc                                                                                                          <  2:514209/139‑1 (MQ=255)     acCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGc                                                                                                          <  2:625031/139‑1 (MQ=255)      ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  1:938515/1‑139 (MQ=255)       cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTc                                                                                                        >  2:1290779/1‑139 (MQ=255)         tCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              <  2:383701/131‑1 (MQ=255)         tCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              >  1:383701/1‑131 (MQ=255)         tCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGa                                                                                                      >  1:1027801/1‑139 (MQ=255)         tCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGa                                                                                                      >  2:1655562/1‑139 (MQ=255)          cggtgacAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcgtgcg                                                                                                                        <  2:1485819/114‑1 (MQ=255)          cggtgacAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcgtgcg                                                                                                                        >  1:1485819/7‑120 (MQ=255)          cAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                                     >  2:435406/1‑139 (MQ=255)                   ccTCGGTCACCGCGCAACCCGGCATCGATGAGAGCGGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  2:292254/1‑126 (MQ=255)                   ccTCGGTCACCGCGCAACCCGGCATCGATGAGAGCGGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  1:292254/126‑1 (MQ=255)                   ccTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAg                                                                                            >  1:152045/1‑139 (MQ=255)                         tCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              <  2:664500/115‑1 (MQ=255)                         tCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              >  1:664500/1‑115 (MQ=255)                             cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGata                                                                                                    <  1:412979/121‑1 (MQ=255)                             cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGata                                                                                                    >  2:412979/1‑121 (MQ=255)                              gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  2:366222/115‑1 (MQ=255)                              gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  1:366222/1‑115 (MQ=255)                                   aCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCgg                                                                                      >  2:1673044/1‑129 (MQ=255)                                   aCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCgg                                                                                      <  1:1673044/129‑1 (MQ=255)                                            cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  1:1290779/139‑1 (MQ=255)                                              aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  2:1668902/137‑1 (MQ=255)                                              aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   >  1:1668902/1‑137 (MQ=255)                                              aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGa                                                                 >  2:1648028/1‑139 (MQ=255)                                                   agCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCt                                                            >  1:1572553/1‑139 (MQ=255)                                                       gAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAt                                                        <  1:435406/139‑1 (MQ=255)                                                        aaTACGGCTCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAtt                                                       <  2:1027801/139‑1 (MQ=255)                                                        aaTACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAtt                                                       <  1:1655562/139‑1 (MQ=255)                                                             ggCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtgt                                                     <  2:663601/136‑1 (MQ=255)                                                             ggCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtgt                                                     >  1:663601/1‑136 (MQ=255)                                                               cgcAGTGCGTTACCCAGCGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                <  2:152045/139‑1 (MQ=255)                                                                              agagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   >  1:289688/1‑105 (MQ=255)                                                                              agagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  2:289688/105‑1 (MQ=255)                                                                                     ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACt                                                                                                             >  1:585820/1‑56 (MQ=255)                                                                                     ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACt                                                                                                             <  2:585820/56‑1 (MQ=255)                                                                                     ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACt                          >  1:1275824/1‑139 (MQ=255)                                                                                     ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACt                          <  2:938515/139‑1 (MQ=255)                                                                                          aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  1:1061420/55‑1 (MQ=255)                                                                                          aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  2:1061420/1‑55 (MQ=255)                                                                                           aGCCACGCTCTAAAGGCTCACGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCCACCTGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgaatt                    >  2:782945/1‑135 (MQ=255)                                                                                           aGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  1:849590/1‑54 (MQ=255)                                                                                           aGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  2:849590/54‑1 (MQ=255)                                                                                           aGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                                       >  2:298981/1‑72 (MQ=255)                                                                                           aGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                                       <  1:298981/72‑1 (MQ=255)                                                                                           aGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGt                    >  1:106315/1‑139 (MQ=255)                                                                                                 gCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   >  2:993844/1‑86 (MQ=255)                                                                                                 gCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  1:993844/86‑1 (MQ=255)                                                                                                             cAAGGATCATCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                                       <  2:504367/54‑1 (MQ=255)                                                                                                             cAAGGATCATCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                                       >  1:504367/1‑54 (MQ=255)                                                                                                             cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   >  1:361003/1‑74 (MQ=255)                                                                                                             cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  2:361003/74‑1 (MQ=255)                                                                                                             cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                <  1:24666/93‑1 (MQ=255)                                                                                                             cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                >  2:24666/1‑93 (MQ=255)                                                                                                             cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgt  >  2:797270/1‑139 (MQ=255)                                                                                                               aGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                                  <  1:1369362/89‑1 (MQ=255)                                                                                                                  |                                                                                                                                      GTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGT  >  NC_000913/3440837‑3441083

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: ‑