Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,179,749	T→C	intergenic (+166/‑247)	rplA → / → rplJ	50S ribosomal subunit protein L1/50S ribosomal subunit protein L10

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,179,749	0	T	C	100.0%	133.3 / NA	38	intergenic (+166/‑247)	rplA/rplJ	50S ribosomal subunit protein L1/50S ribosomal subunit protein L10
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/24);  total (14/24)

GCAAGGCGGCAACTGAGTGAGTCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGTCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTGGCCTATCCAGGCCTCCGTCGAAGACCGC  >  NC_000913/4179673‑4179825                                                                             |                                                                             gCAAGGCGGCAACTGAGTGAGTCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCcaca                                                                          <  1:850492/81‑1 (MQ=255) gCAAGGCGGCAACTGAGTGAGTCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCcac                                                                           <  1:850493/80‑1 (MQ=255)      gcggcAACTGAGTGAGTCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTg                                                                     <  2:850494/81‑1 (MQ=255)        ggcAACTGAGTGAGTCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGaa                                                                   >  1:850544/1‑81 (MQ=255)         gcAACTGAGTGAGTCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAAt                                                                  <  1:850495/81‑1 (MQ=255)         gcAACTGAGTGAGTCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAAt                                                                  <  2:850496/81‑1 (MQ=255)             cTGAGTGAGTCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCg                                                              <  2:850497/81‑1 (MQ=255)               gagtgagtCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGAc                                                            >  1:850543/1‑81 (MQ=255)                  tgagtCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGa                                                         <  2:850498/81‑1 (MQ=255)                    agtCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGAt                                                        <  2:850500/80‑1 (MQ=255)                    agtCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGAt                                                        <  2:850499/80‑1 (MQ=255)                        gCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGaaa                                                   <  2:850501/81‑1 (MQ=255)                        gCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGaa                                                    >  1:850538/1‑80 (MQ=255)                         ccAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGaaaa                                                  <  1:850502/81‑1 (MQ=255)                          cAGGAGCATAGCTAACTATGTGACTGGTGCGGATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAg                                                 <  2:850503/81‑1 (MQ=39)                             gAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGAtt                                              <  2:850504/81‑1 (MQ=255)                              aGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGAttt                                             <  2:850505/81‑1 (MQ=255)                              aGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGAttt                                             >  2:850549/1‑81 (MQ=255)                               gCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTg                                            >  1:850540/1‑81 (MQ=255)                                cATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGt                                           <  2:850506/81‑1 (MQ=255)                                cATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGt                                           >  1:850555/1‑81 (MQ=255)                                 aTAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGtt                                          >  1:850547/1‑81 (MQ=255)                                  tAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTc                                         >  1:850546/1‑81 (MQ=255)                                   aGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCg                                        >  2:850539/1‑81 (MQ=255)                                    gCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGt                                       >  2:850551/1‑81 (MQ=255)                                     cTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGtt                                      <  2:850507/81‑1 (MQ=255)                                        aCTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTgg                                    <  1:850508/80‑1 (MQ=255)                                        aCTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTgg                                    <  1:850509/80‑1 (MQ=255)                                         cTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAg                                  <  1:850510/81‑1 (MQ=255)                                            tgtgACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCt                               <  1:850511/81‑1 (MQ=255)                                             gtgACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTg                              >  2:850559/1‑81 (MQ=255)                                              tgACTGGTGCGAATGAAGGAAGCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTg                              <  1:850512/80‑1 (MQ=255)                                                             aaggaagCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTGGCCTATCCAGGCCTcc              >  2:850545/1‑81 (MQ=255)                                                                gaagCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTGGCCTATCCAGGCCTCCGTc           >  1:850552/1‑81 (MQ=255)                                                                  agCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTGGCCTATCCAGGCCTCCGTCGa         <  2:850513/81‑1 (MQ=255)                                                                   gCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTGGCCTATCCAGGCCTCCGTCg          >  1:850554/1‑79 (MQ=255)                                                                   gCCAACGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTGGCCTATCCAGGCCTCCGTCGaa        <  1:850514/81‑1 (MQ=255)                                                                       aCGCCGCCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTGGCCTATCCAGGCCTCCGTCGAAGAcc    <  1:850515/81‑1 (MQ=255)                                                                         gccgccACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTGGCCTATCCAGGCCTCCGTCGAAGACCGc  <  1:850516/81‑1 (MQ=255)                                                                             |                                                                             GCAAGGCGGCAACTGAGTGAGTCGCCAGGAGCATAGCTAACTATGTGACTGGTGCGAATGAAGGAAGCCAACGCCGTCACAAGCTGAATAGCGACGGATAGAAAAGATTTGTTCGTTGGAGCCTGGCCTATCCAGGCCTCCGTCGAAGACCGC  >  NC_000913/4179673‑4179825

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