Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,231,599	C→T	M11I (ATG→ATA)	yjbD ←	DUF3811 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,231,599	0	C	T	100.0%	108.3 / NA	32	M11I (ATG→ATA)	yjbD	DUF3811 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/16);  total (16/16)

CAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCT  >  NC_000913/4231463‑4231710                                                                                                                                         |                                                                                                                cAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTAt                                                                                                               <  1:395330/139‑1 (MQ=255)                      aTTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACgggg                                                                                          <  2:362116/139‑1 (MQ=255)                                         cttAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGccc                                                                       >  1:370572/1‑139 (MQ=255)                                           tAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTc                                                                     >  1:44012/1‑139 (MQ=255)                                              tGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTggg                                                                                                    <  1:114362/105‑1 (MQ=255)                                              tGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTggg                                                                                                    >  2:114362/1‑105 (MQ=255)                                                                tGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTCAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTggg                                                                                                    <  2:408510/87‑1 (MQ=255)                                                                tGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTggg                                                                                                    >  1:408510/1‑87 (MQ=255)                                                                  gTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGcc                                                                                   >  1:213179/1‑102 (MQ=255)                                                                  gTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGcc                                                                                   <  2:213179/102‑1 (MQ=255)                                                                            gACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGtt                                                     >  1:254833/1‑122 (MQ=255)                                                                            gACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGtt                                                     <  2:254833/122‑1 (MQ=255)                                                                                 tGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACggggg                                                                                         >  1:163326/1‑81 (MQ=255)                                                                                 tGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACggggg                                                                                         <  2:163326/81‑1 (MQ=255)                                                                                 tGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGc                               >  2:6076/1‑139 (MQ=255)                                                                                 tGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGc                               >  1:53781/1‑139 (MQ=255)                                                                                    gcAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAg                            <  2:446872/139‑1 (MQ=255)                                                                                      aaTCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCa                          <  1:208034/139‑1 (MQ=255)                                                                                        tCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGa                        >  2:14412/1‑139 (MQ=255)                                                                                        tCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGa                        >  1:344068/1‑139 (MQ=255)                                                                                            cgcCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTg                    >  1:270757/1‑139 (MQ=255)                                                                                                            gtTTTAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGc    <  2:44012/139‑1 (MQ=255)                                                                                                             ttttAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCatgatg                                                                             <  1:258769/65‑1 (MQ=255)                                                                                                             ttttAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCatgatg                                                                             >  2:258769/1‑65 (MQ=255)                                                                                                             ttttAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGtt                                                     >  1:455908/1‑89 (MQ=255)                                                                                                             ttttAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGtt                                                     <  2:455908/89‑1 (MQ=255)                                                                                                             ttttAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGcc   >  2:240936/1‑139 (MQ=255)                                                                                                             ttttAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGcc   <  2:270757/139‑1 (MQ=255)                                                                                                             ttttAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGcc   <  2:159066/139‑1 (MQ=255)                                                                                                             ttttAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGc    <  2:207684/138‑1 (MQ=255)                                                                                                             ttttAAGTTCACGCTGCTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGc    >  1:207684/1‑138 (MQ=255)                                                                                                              tttAAGTTCACGCTGGTCGCGTTCAGTTATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCt  <  1:6076/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                CAGCCTCACGCTCAACCATTAATTTATCGATGTAATCTTTCTTAATGCTGTTGGTTTCGGAGTTGGTCAATACGCGACCATGCGCAATCCGCGCCCGATCCAGCAGTGTTTTAAGTTCACGCTGCTCGCGTTCAGTCATCTCTTTTTGGGTAATACGGGGGAGTGCCATGATGGTGCCCTCAGTGAGCCAGAGAGTTCAGTGTACGGGAATTCGTTTGCGAGCAGAACTGATAATATTCATTTTGCCT  >  NC_000913/4231463‑4231710

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