Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,446,311	T→C	S79G (AGC→GGC)	rpsH ←	30S ribosomal subunit protein S8

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,446,311	0	T	C	100.0%	150.8 / NA	41	S79G (AGC→GGC)	rpsH	30S ribosomal subunit protein S8
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/14);  total (27/14)

CGCCACCAAGACCAGCCTGGCGCGCTGCACGATCAGTCATAACACCTTTAGAGGTAGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTG  >  NC_000913/3446176‑3446442                                                                                                                                        |                                                                                                                                    cGCCACCAAGACCAGCCTGGCGCGCTGCACGATCAGTCATAACACCTTTAGAGGTAGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGaca                                                                                                                                 >  2:531475/1‑140 (MQ=255) cGCCACCAAGACCAGCCTGGCGCGCTGCACGATCAGTCATAACACCTTTAGAGGTAGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGaca                                                                                                                                 >  1:98713/1‑140 (MQ=255)                                         aaCACCTTTAGAGGTAGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATAcc                                                                                         >  1:1061655/1‑139 (MQ=255)                                           cacCTTTAGAGGTAGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTc                                                                                       >  1:188520/1‑140 (MQ=255)                                              cTTTAGAGGTAGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCaga                                                                                    >  2:1072001/1‑140 (MQ=255)                                               tttAGAGGTAGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCagag                                                                                   >  2:437331/1‑140 (MQ=255)                                                 tAGAGGTAGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTa                                                                                 >  2:653047/1‑140 (MQ=255)                                                      gTAGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTc                                                                            >  1:688597/1‑140 (MQ=255)                                                        aGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCa                                                                          >  2:592030/1‑140 (MQ=255)                                                              aacTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAg                                                                    >  2:131807/1‑140 (MQ=255)                                                                 tGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCt                                                                 >  1:875451/1‑140 (MQ=255)                                                                           gcccccGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGcc                                                       <  2:571141/137‑1 (MQ=255)                                                                           aGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGtt                                                                       <  2:1054190/124‑1 (MQ=255)                                                                           aGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGtt                                                                       >  1:1054190/1‑124 (MQ=255)                                                                            gACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCt                                                      >  2:671221/1‑140 (MQ=255)                                                                             tgccGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCtt                                                     <  2:27001/138‑1 (MQ=255)                                                                             aCCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCtt                                                     <  1:131807/140‑1 (MQ=255)                                                                              cccGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTc                                                    <  1:653047/140‑1 (MQ=255)                                                                                  ccATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCt                                                                 <  1:119964/123‑1 (MQ=255)                                                                                  ccATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCt                                                                 >  2:119964/1‑123 (MQ=255)                                                                                  ccATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACt                                                >  1:527944/1‑140 (MQ=255)                                                                                  ccATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACt                                                >  2:552865/1‑140 (MQ=255)                                                                                  ccATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACt                                                >  1:331576/1‑140 (MQ=255)                                                                                           tCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCt                                       <  1:671221/140‑1 (MQ=255)                                                                                                  cTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATaa                                <  1:552865/140‑1 (MQ=255)                                                                                                     aTCTTTACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAAc                             >  2:887227/1‑140 (MQ=255)                                                                                                          tACGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACcttct                        <  2:471639/140‑1 (MQ=255)                                                                                                            cGTTTATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCGATAAAACCttcttc                      >  2:287469/1‑140 (MQ=255)                                                                                                               ttATAGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTcttcctt                   >  1:1070712/1‑140 (MQ=255)                                                                                                                   aGATGCGCAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGc               <  1:531475/140‑1 (MQ=255)                                                                                                                       gcgcAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGc               >  1:75850/1‑136 (MQ=255)                                                                                                                       gcgcAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGc               <  2:75850/136‑1 (MQ=255)                                                                                                                        cgcAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGc               <  2:657510/135‑1 (MQ=255)                                                                                                                        cgcAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGc               >  1:657510/1‑135 (MQ=255)                                                                                                                        cgcAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGtt          >  2:246590/1‑140 (MQ=255)                                                                                                                        cgcAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGtt          >  1:236544/1‑140 (MQ=255)                                                                                                                          cAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACcttctt                       <  2:665769/125‑1 (MQ=255)                                                                                                                          cAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACcttctt                       >  1:665769/1‑125 (MQ=255)                                                                                                                          cAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGc               <  2:395742/133‑1 (MQ=255)                                                                                                                          cAGACCTGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGc               >  1:395742/1‑133 (MQ=255)                                                                                                                                tGGGCGGCCGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTg  >  2:740050/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                    CGCCACCAAGACCAGCCTGGCGCGCTGCACGATCAGTCATAACACCTTTAGAGGTAGAAACAACTGCGATACCCAGACCCGCCATAACTTTCGGCAGCTCATCTTTACGTTTATAGATGCGCAGACCTGGGCGGCTGACACGCTGAATGCTTTCTACAACAGCTTTGCCCTGGAAATACTTCAGAGTAAGTTCCAGTTCAGGCTTGGTGTCGCCTTCAACTTTAAAATCTTCAATAAAACCTTCTTCCTTCAGCACGTTGGCGATTG  >  NC_000913/3446176‑3446442

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