Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,179,333	C→T	A152V (GCG→GTG)	rplA →	50S ribosomal subunit protein L1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,179,333	0	C	T	100.0%	139.4 / NA	40	A152V (GCG→GTG)	rplA	50S ribosomal subunit protein L1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/19);  total (21/19)

GCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGCGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTG  >  NC_000913/4179198‑4179466                                                                                                                                        |                                                                                                                                      gCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGtt                                                                                                                                   >  2:253705/1‑140 (MQ=255) gCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGtt                                                                                                                                   >  1:310933/1‑140 (MQ=255)    aaaTGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTaaa                                                                                                                                <  2:310933/140‑1 (MQ=255)         aaCTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGc                                                                                                                           <  1:386055/140‑1 (MQ=255)                   ttgttATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGc                                                                                                                 >  1:673350/1‑140 (MQ=255)                             ttCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGtt                                                                                                       <  1:503413/140‑1 (MQ=255)                                tcCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTAcc                                                                                                    >  1:578893/1‑140 (MQ=255)                                    gATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTaa                                                                                                >  1:631999/1‑140 (MQ=255)                                     aTGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAAc                                                                                               >  2:234581/1‑140 (MQ=255)                                            gcgcGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACaaaa                                                                                        >  2:431767/1‑140 (MQ=255)                                            gcgcGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACaaaa                                                                                        <  2:673350/140‑1 (MQ=255)                                             cgcgTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACaaaaa                                                                                       <  1:253705/140‑1 (MQ=255)                                             cgcgTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACaaaaa                                                                                       <  1:253710/140‑1 (MQ=255)                                                  tgttgGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGca                                                                                  <  2:546114/140‑1 (MQ=255)                                                     tgGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGcatca                                                                               >  1:152887/1‑140 (MQ=255)                                                            cTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCacac                                                                        >  1:413677/1‑140 (MQ=255)                                                             tGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCAcacc                                                                       <  2:460086/140‑1 (MQ=255)                                                             tGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCAcacc                                                                       <  2:506400/140‑1 (MQ=255)                                                              gggCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACaaaaa                                                                                       >  1:401096/1‑123 (MQ=255)                                                              gggCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACaaaaa                                                                                       <  2:401096/123‑1 (MQ=255)                                                               ggCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCAcaccac                                                                     <  1:274657/140‑1 (MQ=255)                                                                          gggTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAg                                                          <  2:152887/140‑1 (MQ=255)                                                                               cgcgcgGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGAc                                                     >  2:491026/1‑140 (MQ=255)                                                                                 cgcgGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACtt                                                   >  2:145173/1‑140 (MQ=255)                                                                                                aaCCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACt                                    <  2:30987/140‑1 (MQ=255)                                                                                                aaCCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACt                                    >  2:126601/1‑140 (MQ=255)                                                                                                aaCCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACt                                    <  2:62124/140‑1 (MQ=255)                                                                                                     gAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTgaaag                               >  1:93057/1‑140 (MQ=255)                                                                                                         gTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGaaaa                           <  2:703474/140‑1 (MQ=255)                                                                                                         gTGGGTACTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGaaaa                           <  2:631999/140‑1 (MQ=255)                                                                                                               aCTGTAACACCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGa                     <  2:678016/140‑1 (MQ=255)                                                                                                                      cacCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGaa                    <  2:48288/134‑1 (MQ=255)                                                                                                                      cacCGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGaa                    >  1:48288/1‑134 (MQ=255)                                                                                                                         cGAACGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTctgctg           >  2:191148/1‑140 (MQ=255)                                                                                                                           aaCGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGt         >  2:601685/1‑140 (MQ=255)                                                                                                                           aaCGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGt         >  2:664832/1‑140 (MQ=255)                                                                                                                           aaCGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGt         >  1:260653/1‑140 (MQ=255)                                                                                                                             cGTTGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTg       >  2:308177/1‑140 (MQ=255)                                                                                                                               ttGCTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTgcg     <  1:601685/140‑1 (MQ=255)                                                                                                                                  cTGAAGTGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTg  >  1:622482/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                      GCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAACGTTGCTGAAGCGGTTAAAAACGCTAAAGCTGGCCAGGTTCGTTACCGTAACGACAAAAACGGCATCATCCACACCACCATCGGTAAAGTGGACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTG  >  NC_000913/4179198‑4179466

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