Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,182,817	A→G	T525A (ACG→GCG)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,182,817	0	A	G	100.0%	125.4 / NA	33	T525A (ACG→GCG)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/15);  total (18/15)

CTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATC  >  NC_000913/4182685‑4182951                                                                                                                                     |                                                                                                                                       ctGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGcaca                                                                                                                                  <  1:53739/139‑1 (MQ=255) ctGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGcaca                                                                                                                                  <  1:362107/139‑1 (MQ=255)                  ccTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCa                                                                                                                 >  2:268642/1‑139 (MQ=255)                    tGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACt                                                                                                               >  1:339685/1‑139 (MQ=255)                     gATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTc                                                                                                              >  1:201316/1‑139 (MQ=255)                        gCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGc                                                                                                           >  1:361057/1‑139 (MQ=255)                            caGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAg                                                                                                       >  2:611507/1‑139 (MQ=255)                                  atGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTc                                                                                                 <  1:383746/139‑1 (MQ=255)                                    gATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTg                                                                                               <  2:339685/139‑1 (MQ=255)                                               aGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                                                                                    >  1:216818/1‑139 (MQ=255)                                                                 tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGcaca                                                                                                                                  <  1:328670/75‑1 (MQ=255)                                                                 tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGcaca                                                                                                                                  >  2:328670/1‑75 (MQ=255)                                                                  gAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCga                                                                 <  1:268642/139‑1 (MQ=255)                                                                        gTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACgca                                                           >  1:486585/1‑137 (MQ=255)                                                                           cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACt                                                                                                               >  2:503768/1‑84 (MQ=255)                                                                           cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACt                                                                                                               <  1:503768/84‑1 (MQ=255)                                                                           cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAAGCgg                                                                                                   <  1:290949/96‑1 (MQ=255)                                                                           cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAAGCgg                                                                                                   >  2:290949/1‑96 (MQ=255)                                                                               ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGa                                                    <  2:361057/139‑1 (MQ=255)                                                                                   ccagccagCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa                                                >  2:581038/1‑139 (MQ=255)                                                                                      gccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCAcca                                             >  1:224440/1‑137 (MQ=255)                                                                                           cTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACCACGGTc                                        >  2:664653/1‑139 (MQ=255)                                                                                                        tGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAg                                                                                                       >  2:167158/1‑63 (MQ=255)                                                                                                        tGGACCAGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAg                                                                                                       <  1:167158/63‑1 (MQ=255)                                                                                                              aGAACAACCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGCACACCCGACTCACCACGGTCGCGTATGTCCAATCGAAAc                     <  2:486585/139‑1 (MQ=255)                                                                                                                aacaacCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACCACGGTCGCGTATGTCCAATCGAAAccc                   >  2:601647/1‑139 (MQ=255)                                                                                                                    acCCGCTGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACCACGGTCGCGTATGTCCAATCGAAACCCCTGa               <  2:224440/139‑1 (MQ=255)                                                                                                                          tGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTg                                                                                  >  2:25043/1‑66 (MQ=255)                                                                                                                          tGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTg                                                                                  <  1:25043/66‑1 (MQ=255)                                                                                                                          tGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACCACGGTCGCGTATGTCCAATCGa                        >  2:441061/1‑124 (MQ=255)                                                                                                                          tGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACCACGGTCGCGTATGTCCAATCGa                        <  1:441061/124‑1 (MQ=255)                                                                                                                          tGTCTGAGATTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACCACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTcc         >  1:133439/1‑139 (MQ=255)                                                                                                                                 gaTTGCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACCACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATc  <  2:216818/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       CTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATC  >  NC_000913/4182685‑4182951

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