Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_tpiA	4,182,833	C→T	P520L (CCG→CTG)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_vch_tpiA	4,182,833	0	C	T	100.0%	126.7 / NA	39	P520L (CCG→CTG)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/12);  total (27/12)

GTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGT  >  NC_000913_3_vch_tpiA/4182692‑4182972                                                                                                                                              |                                                                                                                                            gTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCt                                                                                                                                      >  2:50024/1‑149 (MQ=255)       tGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGAtt                                                                                                                                >  2:293650/1‑149 (MQ=255)       tGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTAAGAtt                                                                                                                                >  1:152163/1‑149 (MQ=255)          aaGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTTCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACg                                                                                                                             <  2:633701/149‑1 (MQ=255)                 tctgtctCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAAc                                                                                                                      <  2:627132/149‑1 (MQ=255)                    gtctCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAAcgtc                                                                                                                   >  2:525400/1‑149 (MQ=255)                       tctGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTa                                                                                                                >  1:202027/1‑149 (MQ=255)                          gggCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGCTTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGGACAGCCTGCTGTCTGAGATTACCCACAAACGTCGTAtct                                                                                                             >  2:312558/1‑149 (MQ=255)                               aTCTGGATACCCTGAGGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCa                                                                                                        >  1:507979/1‑149 (MQ=255)                                  tGGATACCCTGATGCCACAGGATATGATCGACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGATGTTCGGTTCCAGCCAGCTGTCTCAGTTGATGGACCAGAACAACCTGCTGCCTGAGATTACGCACAAACGTCGTATCTCCGCACTc                                                                                                     >  1:43690/1‑149 (MQ=255)                                  tGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTc                                                                                                     >  1:170826/1‑149 (MQ=255)                                      tACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCGGTTTATGGACCAGACCAACCTGCTGTCTGAGTTTGCGCACAACCATCGTGTCTCCGCACTCGGcc                                                                                                 >  2:375359/1‑149 (MQ=255)                                                 cacaTGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGCTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTg                                                                                      >  2:535725/1‑149 (MQ=255)                                                      gATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                                                                                 <  1:18738/149‑1 (MQ=255)                                                       atatGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCAGCCCAGGCGGTCTGACCCGt                                                                                >  1:426055/1‑149 (MQ=255)                                                         atGATCACCGCCAAGGCGATTTCCGCAGCAGTGAAAGAGTTCGTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGa                                                                              <  1:262366/149‑1 (MQ=255)                                                            aTCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTATATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                                                                           >  1:586918/1‑149 (MQ=255)                                                            aTCAACGCCAAGCCGATTTCCGCAGAAGTGAAAGAGTTCTTCGGTTGCAGCCAGCTGTCTCAGTTTGCGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCTGCCCAGGCGGTCTGACCCCTGAACg                                                                           >  1:609758/1‑149 (MQ=255)                                                                  gCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAgg                                                                     >  1:653078/1‑149 (MQ=255)                                                                       gCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGACCGTGGAGCCTTCg                                                                >  1:572479/1‑149 (MQ=255)                                                                         cGATTTCCGCAGCAGTGAAAGAGTTCGTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGaa                                                              <  1:293650/149‑1 (MQ=255)                                                                               ccGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCga                                                        <  1:388726/149‑1 (MQ=255)                                                                                          aaaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAccc                                             >  2:442073/1‑149 (MQ=255)                                                                                           aaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTAAACCCg                                            >  1:626190/1‑149 (MQ=255)                                                                                            agagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGa                                           <  1:205143/149‑1 (MQ=255)                                                                                            agagTTCTGCGGTTCCAGCCAGCTGTCTCAGTTTATGGCCCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGa                                           <  1:204651/149‑1 (MQ=255)                                                                                                   ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGCTTACGCCAAAACGTAGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCATGAACGAGCAGGCTTCGAAGTTCGAGACGTGCACCCGACTCACTa                                    >  2:659151/1‑149 (MQ=255)                                                                                                   ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTa                                    >  2:111743/1‑149 (MQ=255)                                                                                                       gTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGt                                >  1:72289/1‑149 (MQ=255)                                                                                                       gTTCCAGCCACCTGTCCCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGt                                <  1:334503/149‑1 (MQ=255)                                                                                                          ccagccagCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgc                             <  1:76970/149‑1 (MQ=255)                                                                                                               cagcTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATg                        <  2:653078/149‑1 (MQ=255)                                                                                                                      ctcAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATc                 >  2:354222/1‑149 (MQ=255)                                                                                                                         aGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa              >  2:410109/1‑149 (MQ=255)                                                                                                                             tATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTATGACCCGTGAACGTGCAGGCTTCGAAGTTCGAAACGTACACACGGCTCACTACGGTCGCGTATGTCCGATcagatcgg          >  1:111306/1‑142 (MQ=255)                                                                                                                                 gACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGa      <  1:398884/149‑1 (MQ=255)                                                                                                                                  aCCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGaa     >  1:190954/1‑149 (MQ=255)                                                                                                                                    cAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAgt   >  2:354893/1‑148 (MQ=255)                                                                                                                                     aGAACAACCTGCTGTCTGCGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGTTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGt  >  2:61780/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                            GTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGT  >  NC_000913_3_vch_tpiA/4182692‑4182972

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑