Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,137,229	A→G	S143S (AGT→AGC)	yijF ←	DUF1287 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,229	0	A	G	100.0%	169.0 / NA	47	S143S (AGT→AGC)	yijF	DUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/24);  total (23/24)

CTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAG  >  NC_000913/4137096‑4137353                                                                                                                                      |                                                                                                                             cttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTggg                                                                                                                         <  2:1314621/139‑1 (MQ=255)       gcgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCt                                                                                                                    >  1:897687/1‑138 (MQ=255)       gcgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCt                                                                                                                    <  2:897687/138‑1 (MQ=255)        cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTg                                                                                                                  >  2:702303/1‑139 (MQ=255)            ccTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtgg                                                                                                              <  1:786121/139‑1 (MQ=255)                             tATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCa                                                                                             >  2:607218/1‑139 (MQ=255)                                  ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                        >  2:1254551/1‑123 (MQ=255)                                  ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                        <  1:1254551/123‑1 (MQ=255)                                  ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg                                                                                        >  2:406033/1‑139 (MQ=255)                                       ggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTaa                                                                                   >  1:483548/1‑139 (MQ=255)                                            aCCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAg                                                                              <  2:483548/139‑1 (MQ=255)                                                tCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggt                                                                                                             >  1:1366226/1‑104 (MQ=255)                                                tCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggt                                                                                                             <  2:1366226/104‑1 (MQ=255)                                                     ggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGa                                                                     >  2:1042210/1‑139 (MQ=255)                                                     ggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGCTATCCCGTCGCCAGGTGACAATATCGCCCGCATGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGa                                                                     >  1:667535/1‑139 (MQ=255)                                                             cATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCAc                                                             <  1:879221/139‑1 (MQ=255)                                                                      ccaccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTg                                                    <  1:1225921/139‑1 (MQ=255)                                                                         ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCt                                                                                                                            <  2:223954/64‑1 (MQ=255)                                                                         ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCt                                                                                                                            >  1:223954/1‑64 (MQ=255)                                                                         ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                           >  1:1301025/1‑65 (MQ=255)                                                                         ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                           <  2:1301025/65‑1 (MQ=255)                                                                          cccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCg                                                <  1:406033/139‑1 (MQ=255)                                                                          cccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCg                                                >  2:268639/1‑139 (MQ=255)                                                                             aaTATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTgggg                                                                                                                        <  1:488791/64‑1 (MQ=255)                                                                             aaTATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTgggg                                                                                                                        >  2:488791/1‑64 (MQ=255)                                                                                   tgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc                                                                                >  1:372198/1‑98 (MQ=255)                                                                                   tgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc                                                                                <  2:372198/98‑1 (MQ=255)                                                                                    gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGt                                      >  1:521993/1‑139 (MQ=255)                                                                                      cAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                           >  2:1403676/1‑52 (MQ=255)                                                                                      cAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                           <  1:1403676/52‑1 (MQ=255)                                                                                        gCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCgg                                  <  1:1042210/139‑1 (MQ=255)                                                                                         cccGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                        <  1:986489/68‑1 (MQ=255)                                                                                         cccGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                        >  2:986489/1‑68 (MQ=255)                                                                                           cGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                         <  2:224222/65‑1 (MQ=255)                                                                                           cGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                         >  1:224222/1‑65 (MQ=255)                                                                                              tATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCt                                                                                                                            >  2:942039/1‑43 (MQ=255)                                                                                              tATCCAGTCGCCAGAAGACAATATCGCCCGCTTGATAGTCGCt                                                                                                                            <  1:942039/43‑1 (MQ=255)                                                                                                    gTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCg                                                                                                       >  2:1411488/1‑58 (MQ=255)                                                                                                    gTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCg                                                                                                       <  1:1411488/58‑1 (MQ=255)                                                                                                         cAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCAtt                 <  1:607218/139‑1 (MQ=255)                                                                                                           ggAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCAtttt               <  1:268639/139‑1 (MQ=255)                                                                                                                cAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCgg          <  2:521993/139‑1 (MQ=255)                                                                                                                      cGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCAtt                 <  2:33642/126‑1 (MQ=255)                                                                                                                      cGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCAtt                 >  1:33642/1‑126 (MQ=255)                                                                                                                        cccGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGtttt                                                                                                 >  2:1154983/1‑44 (MQ=255)                                                                                                                        cccGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGtttt                                                                                                 <  1:1154983/44‑1 (MQ=255)                                                                                                                        cccGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAg  >  2:899876/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                             CTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAG  >  NC_000913/4137096‑4137353

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