Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,168,461	G→T	intergenic (+14/‑180)	gltT → / → rrlB	tRNA‑Glu/23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,168,461	0	G	T	100.0%	110.1 / NA	33	intergenic (+14/‑180)	gltT/rrlB	tRNA‑Glu/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/17);  total (16/17)

CGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCT  >  NC_000913/4168353‑4168589                                                                                                             |                                                                                                                                 cgTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCCAGGGGACGCCACTTGCTGGGTTGTTAGTGAAAGTCGCCGACCTTAATATCTCaaaa                                                                                                   >  1:731871/1‑140 (MQ=12)                tCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGa                                                                                    <  1:619689/140‑1 (MQ=255)                 cGTGTCCCCTTCGTCTAGATGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGAt                                                                                   <  2:794671/140‑1 (MQ=255)                 cGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGAt                                                                                   >  2:742334/1‑140 (MQ=255)                       cccTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTga                                                                             <  1:883354/140‑1 (MQ=38)                          ttCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGat                                                                          >  2:227345/1‑140 (MQ=255)                               cTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTg                                                                     <  1:882885/140‑1 (MQ=38)                                   agGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGGTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGGGATGTTTTAGATATTTTctct                                                                 >  2:759202/1‑140 (MQ=11)                                     gCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGAGATTTTCTCttt                                                               >  1:695229/1‑140 (MQ=255)                                       ccAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCCCTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTaa                                                             >  1:288313/1‑140 (MQ=255)                                       ccAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTaa                                                             <  2:288313/140‑1 (MQ=37)                                        cAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGGGATGTTTGAGATATTTGCTCTTTaaa                                                            >  1:414766/1‑140 (MQ=255)                                          ggACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTaaaaa                                                          >  1:663794/1‑140 (MQ=37)                                           gACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAAt                                                         <  2:663794/140‑1 (MQ=37)                                               ccGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTgg                                                     >  1:572863/1‑140 (MQ=37)                                               ccGCCCTTTCACGGCGGTAACAGGGGGTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGGGATGTTTGAGATATTTTCTCTTTAAAAATCTgg                                                     >  2:525376/1‑140 (MQ=255)                                                        cACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTg                                            <  1:759202/140‑1 (MQ=37)                                                             cggTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCTAGCTGAAAAt                                       >  1:353343/1‑140 (MQ=255)                                                             cggTAACAGGGGGTCGAATCCCCTAGGGGGCGCCACTTGCTGGGTTGTTAGTGAAAGGCGCCGGCCTTAATATCTCAAAACTCATCTTTGGGGGATGTTTGAGATATTTGCTCTTTTAAAATCTGGATCAAGCTGAAAAt                                       >  2:507570/1‑140 (MQ=11)                                                               gTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTg                                     <  2:175797/140‑1 (MQ=37)                                                               gTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTg                                     <  1:431516/140‑1 (MQ=37)                                                                        gTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCTAGCTGAAAATTGAAACACTGa                            <  2:353343/140‑1 (MQ=255)                                                                          tCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGaac                          <  1:227345/140‑1 (MQ=37)                                                                          tCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGaac                          <  2:414766/140‑1 (MQ=37)                                                                                 cccTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACgaga                   >  1:397095/1‑140 (MQ=255)                                                                                  ccTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGagag                  <  1:742334/140‑1 (MQ=255)                                                                                  ccTAGGGGACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGagag                  <  2:572863/140‑1 (MQ=255)                                                                                      ggggACGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAgttgt              <  2:737050/140‑1 (MQ=255)                                                                                           cGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg         >  1:854427/1‑140 (MQ=255)                                                                                           cGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg         >  1:665309/1‑140 (MQ=255)                                                                                           cGCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg         >  2:879930/1‑140 (MQ=255)                                                                                            gCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGa        >  2:592771/1‑140 (MQ=255)                                                                                            gCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGa        >  2:613750/1‑140 (MQ=255)                                                                                            gCCACTTGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGa        >  1:29029/1‑140 (MQ=255)                                                                                                  tGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtctct  <  2:822380/140‑1 (MQ=255)                                                                                                  tGCTGGTTTGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtctct  <  2:854427/140‑1 (MQ=255)                                                                                                             |                                                                                                                                 CGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCT  >  NC_000913/4168353‑4168589

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