Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,640,785	C→T	intergenic (+14/‑157)	ytjE → / → yjtD	protein YtjE/putative tRNA/rRNA methyltransferase YjtD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,640,785	0	C	T	100.0%	106.6 / NA	36	intergenic (+14/‑157)	ytjE/yjtD	protein YtjE/putative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/15);  total (21/15)

GAGGGGCCGGACACCCAATAAAACTACGCTTCGTTGACATATATCAAGTTCAATTGTAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAG  >  NC_000913/4640551‑4641027                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   gagGGGCCGGACACCCAATAAAACTACGCTTCGTTGACATATATCAAGTTCAATTGTAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTAc                                                                                                                                                                                                                                      >  2:249780/1‑249 (MQ=255)     ggacGGACACCCAATAAAACTACGCTTCGTTGACATATATCAAGTTCAATTGTAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAg                                                                                                                                                                                                                                <  2:330671/248‑1 (MQ=255)                    aaaaCTACGCTTCGTTGACATATATCAAGTTCAATTGTAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGcc                                                                                                                                                                                                                  <  1:358951/250‑1 (MQ=255)                     aaaCTACGCTTCGTTGACATATATCAAGTTCAATTGTAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGCTATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTAAGGAAATAGTTTAGCCtt                                                                                                                                                                                                                >  2:676506/1‑251 (MQ=255)                          aCGCTTCGTTGACATATATCAAGATCAATTGTAGCACGTTAACCGTTTGCTGAAATCATCGAATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCtttttt                                                                                                                                                                                                            >  2:434726/1‑250 (MQ=255)                                   tGACATATATCAAGTTCAATTGTAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGt                                                                                                                                                                                                   <  1:631446/250‑1 (MQ=255)                                      cATATATCAAGTTCAATTGAAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTCGCATTTTTTAACAATTCAACCGTTAGTACAGTAAGGAAATAGTTTAGCCTTTTTTAAGCTAAgaaa                                                                                                                                                                                                 >  2:451272/1‑246 (MQ=255)                                            tCAAGTTCAATTGTAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCttt                                                                                                                                                                                         >  1:428681/1‑251 (MQ=255)                                                  tCAATTGTAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTg                                                                                                                                                                                    <  1:266467/250‑1 (MQ=255)                                                     aTTGTAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGATATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCAGTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGAc                                                                                                                                                                                <  2:22974/251‑1 (MQ=255)                                                           gCACGTTAACAGATTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGTTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTCCAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAGGAAAAGGGCTTTTTCTGCGACTTaagt                                                                                                                                                                          >  2:318058/1‑248 (MQ=255)                                                                          gatgaAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAAt                                                                                                                                                           >  1:612083/1‑251 (MQ=255)                                                                                         tCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGtaat                                                                                                                                             >  1:649273/1‑250 (MQ=255)                                                                                                               attaatTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCaacctgtt                                                                                                                      <  2:428681/251‑8 (MQ=255)                                                                                                                       aataatCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTa                                                                                                              <  1:641435/251‑1 (MQ=255)                                                                                                                             ccAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTgg                                                                                                         <  1:418381/250‑2 (MQ=255)                                                                                                                             ccAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAA‑CAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTAAAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGAAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTggtc                                                                                                       >  2:401881/1‑248 (MQ=255)                                                                                                                                            atACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTCCGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCa                                                                                          >  1:423243/1‑250 (MQ=255)                                                                                                                                               cAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAA‑CAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATg                                                                                     <  1:401881/251‑1 (MQ=255)                                                                                                                                                               tgAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACCTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGACGAGTGGCAATATGCGTATAACGattatt                                                                      >  1:359485/1‑251 (MQ=255)                                                                                                                                                                  aTCCAATTGTAGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTATGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGCAGAGTGGCAATATGCGTATAACGATTATTCTg                                                                   >  2:194548/1‑251 (MQ=255)                                                                                                                                                                     cAATTGTCGAGATTTATTTTTTATAAAATTATCCTACGTAAACAGAAGGATATGTAGCATTTTTTAACAATTAAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTATGTTATGGTTTTGTAAATTCGCACCGCGTACTAAGTTGACAGTGATCACCCTGCTCGCGGTTCTTTGGACAAGCAGAGTGGCAATATGCGTATAACGATTATTCTGGTc                                                                >  2:419067/1‑251 (MQ=255)                                                                                                                                                                     cAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATCAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGtt                                                                >  2:507843/1‑250 (MQ=255)                                                                                                                                                                            cGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGAACCCg                                                         >  1:428986/1‑251 (MQ=255)                                                                                                                                                                                                 ttATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGCTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGTGCATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAACATATTgggg                                    >  1:588738/1‑251 (MQ=255)                                                                                                                                                                                                  tATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGc                                   >  1:235677/1‑251 (MQ=255)                                                                                                                                                                                                   aTCCTAAGGAAACAGAAGGATAGTTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATTGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACATTGATCCCCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCa                                  <  2:574766/251‑1 (MQ=255)                                                                                                                                                                                                        aaggaaaCAGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTCTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAgcg                               <  1:15524/245‑1 (MQ=255)                                                                                                                                                                                                            aaaCAGACGGATATTTAGAATATTTTAACACTTCAACCGTTAGTAACGTCGGGAAATAGTTAATCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGGAAATACGGACCGCGTGATAAGTTGACAGTGATCACACTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCATTATGCGGATCACGATTATTCTGGCCGCACCCGCCAGAGCAGAAGATATTgtagcagaggggagg                         >  2:193851/1‑237 (MQ=255)                                                                                                                                                                                                              aCAGAAGGATATGTAGCATTTTTTAACAGTTCAACCGTTAGTACAGTCCGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAACGGGCTTTTTCTGCGACTTa                                                                                                                                                                             >  2:41084/1‑101 (MQ=255)                                                                                                                                                                                                              aCAGAAGGATATGTAGCATGTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTa                                                                                                                                                                             <  1:41084/101‑1 (MQ=255)                                                                                                                                                                                                                aGAAGGATATGTAGCATTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCa                      >  1:42098/1‑250 (MQ=255)                                                                                                                                                                                                                    ggATATGTATCATTTTTTAACAATTCAACCGTTAGGACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGTTGAAAGTGAACACCCTGTTCGCGGTTATTTGGTCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGTGCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGaa                 <  2:649273/251‑1 (MQ=255)                                                                                                                                                                                                                               aTTTTTTAACAATTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGtt                                                                                                                                        >  1:418582/1‑121 (MQ=255)                                                                                                                                                                                                                               aTTTTTTAACAATGCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGGATTTGTAAATTCGCACCGCGTAATAAGtt                                                                                                                                        <  2:418582/121‑1 (MQ=255)                                                                                                                                                                                                                                    ttAACAATTCACCCGTTAGTAAAGTCAGGAAATAGTTTAGCTTTGTGTAAGCTAAGTAAAGGGCGTTTGCTTCGATTTACGTTAAGGTTTGGTAAATTCGCACCGCGTAATAAGTGGCCAGTGATCACCCTGTACGCGGTTATTTGGTCAAGACGAGTGGAAATATGTGTATAACGATGATTCTGGCCGCACCCGCTAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAg  <  2:102547/250‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   GAGGGGCCGGACACCCAATAAAACTACGCTTCGTTGACATATATCAAGTTCAATTGTAGCACGTTAACAGTTTGATGAAATCATCGTATCTAAATGCTAGCTTTCGTCACATTATTTTAATAATCCAACTAGTTGCATCATACAACTAATAAACGTGGTGAATCCAATTGTCGAGATTTATTTTTTATAAAATTATCCTAAGTAAACAGAAGGATATGTAGCATTTTTTAACAACTCAACCGTTAGTACAGTCAGGAAATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAG  >  NC_000913/4640551‑4641027

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