Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,236,119	A→G	T72A (ACC→GCC)	yjbE →	uncharacterized protein YjbE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,236,119	0	A	G	100.0%	144.4 / NA	40	T72A (ACC→GCC)	yjbE	uncharacterized protein YjbE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/20);  total (20/20)

CGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCG‑CTTGC‑GGCAACCGGTGGCGGTGATGGTTCTAATACC‑GGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAG  >  NC_000913/4235986‑4236254                                                                                                                                         |                                                                                                                                        cGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGGGGTCCTCCCCCACCAGCGGCAGCACCGTAAACTCCGCGGGGGGGGGTGCGCCTTGC‑GCCACCCGGGGGGGGGGTGGGTTTAAAAACC‑GGGGCCCCccca                                                                                                                                    >  2:1553965/1‑140 (MQ=255)                ggCAGCACCGTGGGTTGGGGGGGGGGGTCCCCCCCCGCCCGGCGCCAGCCCCGTAGGCGGGGGGGGGGGTGCCGCG‑CTGGC‑GGCACCGGGGGGGGGTG‑TGGTTTAAAACCCGGGGACCGCCACAACCACGACCACCAGTa                                                                                                                    <  2:977896/140‑1 (MQ=255)                  cagcaCCGGGGGTTGGGGGGGGGAGTCCCCCCGCGCCCGGCGCCAGCACCGTAGGCGGGGGGGGGGGTGCCGCG‑CTGGC‑GGCACCGGGGGGGGGGGATGTTTCTAATCCG‑GGGACCGCCACAACCACGACCACCAGTAcc                                                                                                                  <  2:936988/140‑1 (MQ=255)                     caACGTCGGTTTCGGCGGGGAGGTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGGGGGGGTCGCG‑CCTGCGGGCAACCGG‑GGCGGGGGGGGGTTTTAAACC‑GGGGCCGCCCCCACCCCCACCACCCGTAccccg                                                                                                               >  2:1064159/1‑138 (MQ=255)                        gggcGGTTGGGGGGGGGGGTCCCCCCGCGCCCGGCGTCGGCCCGGTAGGCCGGGGGGGGGGTGCCGCG‑CTGGG‑GGCACCGGGGGGGGGTGTTGTTTCAAATCCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGtaa                                                                                                            <  2:378761/137‑1 (MQ=255)                        ggTCGTTTGCGGGGGGTAGCCCCCCGGCGCCCGGCGTCAGCCCGGTAGGCGGGGGGGGGGGGGGCGCG‑CTTGG‑GGCACCGGGGGGGGGTGATGGTTAAAATACG‑GGGTCCGCCACAACCACGACCACCAGTACCCAGtaa                                                                                                            <  2:1457424/139‑1 (MQ=255)                                 cggcggGGGGGTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGGGGGGGGCGCG‑CTTGCGGGCAACCGG‑GGGGGGGGTGGGTCTTATACC‑CGGGCCCCCCCAACCCCGCCCCCCAGGACCCCACAATaaaaaatg                                                                                                   >  2:194364/1‑135 (MQ=255)                                  gggggggCGTCCCCCCGCGACCGCCGCCAGCCCCGTAAGCGGGGGGGGGGGTGGCGCG‑CTGGC‑GGCACCGGGGGGGGGTGATGTTTCAAATCCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAgtatgt                                                                                                  <  2:1363899/137‑1 (MQ=255)                                                 cgcgCCCAGCTCGAGCCCCGTAGGCTGGGGGGTGGGGGCCGCG‑CTGGC‑GGCGCCGGGGGGGGGTGATGTTTCAAACCCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAAtt                                                                                   <  2:1110216/140‑1 (MQ=255)                                                 cgcgACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCG‑CTTGC‑GGCAACCGGTGGCGGTGATGGTTCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAAtt                                                                                   >  2:942759/1‑140 (MQ=255)                                                   cgACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCG‑CTTGC‑GGCAACCGGTGGCGGTGATGGTTCTAATACA‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCg                                                                                 >  1:1472055/1‑140 (MQ=255)                                                        aGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCG‑CTTGC‑GGCAACCGGTGGCGGTGATGGTTCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCa                                                                            >  2:286534/1‑140 (MQ=255)                                                              cgccccGTAAGCGGGGGGGGGGGTGCCGCG‑CTGGC‑GGCACCGGGGGGGGGTGTTGTTTCAAACACG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCAt                                                                      <  2:961405/136‑1 (MQ=255)                                                              agcccggTAAGCTGGGGGGGGGGTGCCGCG‑CTGGC‑GGCACCGGGGGGGGGTGTTGTTTTAAATCCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCAt                                                                      <  2:723353/134‑1 (MQ=255)                                                                     tAAGCTCGGCGGTGGGTGTCGCG‑CTTGC‑GGCAACCGGTGGCGGTGATGGTTCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGc                                                               >  1:1098131/1‑140 (MQ=255)                                                                      aaGCTCGGCGGTGGGTGTCGCG‑CTTGC‑GGCAACCGGTGGCGGGGATGGTTCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCt                                                              >  2:377216/1‑140 (MQ=255)                                                                       ggCTGGGGGGGGGGTGCCGCG‑CTGGC‑GCCACCGGGGGGGGGTGTTGTTTCAAATACG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTa                                                             <  2:343729/139‑1 (MQ=255)                                                                           cggcggTGGGTGTCGCG‑CTTGC‑GGCAACCGGTGGCGGTGATGGTTCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGGAATAAAGTATGTATCCCCAAAATAATTCGAGGCATTGCCTCTGTGGCTAGCAg                                                         >  2:162234/1‑140 (MQ=255)                                                                             gcggTGGGTGTCGCG‑CTTGC‑GGCAACCGGTGGCGGTGATGGTTCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTa                                                       >  1:296670/1‑140 (MQ=255)                                                                                gtgggtgTCGCG‑TTTGG‑GGCATGCGGGGGGGGTGATGGTTATAATCCC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGa                                                    <  1:30130/140‑1 (MQ=255)                                                                                 ggggtgTCGCG‑CTGGC‑GGCACCGGGGGGGGGTGATGTTTTAAATCCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGaa                                                   <  1:1553965/139‑1 (MQ=255)                                                                                      gccgcg‑CTGGC‑GGCACCGGGGGGGGGTGATGGTTCTAATCCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGAt                                              <  2:666957/138‑1 (MQ=255)                                                                                       tCGCG‑CTTGC‑GGCAACCGGTGGCGGTGATGGTTCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGAtt                                             >  1:939721/1‑140 (MQ=255)                                                                                        cgcg‑CTTGC‑GGCAACCGGTGGCGGGGATGGTTCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGCATGAACGGATta                                            >  2:1158845/1‑140 (MQ=255)                                                                                               tGC‑GGCAACCGGTGGCGGTGATGGTTCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATtaaccata                                      >  2:123848/1‑140 (MQ=255)                                                                                                         gggggggggTGATGTTTCAAAACCG‑GGGCCCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATAAAGGGATTAACCATAACCACACTc                             <  2:1189633/133‑1 (MQ=255)                                                                                                            ggggggTGTTGTTTCAAATCCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCgg                          <  2:1093365/136‑1 (MQ=255)                                                                                                            ggggggTGATGTTTCTAATCCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCgg                          <  1:687367/136‑1 (MQ=255)                                                                                                             gggggTGATGGTTCTAACCCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGgt                         <  1:1562684/137‑1 (MQ=255)                                                                                                               cGGTGATGGTTCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAAAAAATCGAGTCATTGCATCTGTGGCTAGCAGTATGAAGGGATTAACCATAACCACACTCCGgtgt                       >  2:172244/1‑140 (MQ=255)                                                                                                                ggTGTTGGTTCTAATCCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGtgtg                      <  1:1338130/140‑1 (MQ=255)                                                                                                                  tGATGGTTCAAATACG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGt                    <  1:810044/140‑1 (MQ=255)                                                                                                                        ttCTAATACC‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCt              >  1:1151503/1‑140 (MQ=255)                                                                                                                          ttAATGCG‑GGGACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGc            <  2:296670/139‑1 (MQ=255)                                                                                                                                   gggACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTgg    >  2:341092/1‑140 (MQ=255)                                                                                                                                   gggACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTgg    >  2:642331/1‑140 (MQ=255)                                                                                                                                   gggACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTgg    >  2:1417894/1‑140 (MQ=255)                                                                                                                                   gggACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTgg    >  1:39854/1‑140 (MQ=255)                                                                                                                                   gggACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTgg    >  1:1464238/1‑140 (MQ=255)                                                                                                                                    ggACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGga   >  2:177975/1‑140 (MQ=255)                                                                                                                                    ggACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGga   >  1:76500/1‑140 (MQ=255)                                                                                                                                     gACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGgag  >  2:569138/1‑140 (MQ=255)                                                                                                                                     gACCGCCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGgag  <  2:947968/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        CGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCG‑CTTGC‑GGCAACCGGTGGCGGTGATGGTTCTAATACC‑GGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAG  >  NC_000913/4235986‑4236254

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