Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,449,996	T→C	E18E (GAA→GAG)	yfcO ←	DUF2544 domain‑containing protein YfcO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,449,996	0	T	C	100.0%	112.8 / NA	32	E18E (GAA→GAG)	yfcO	DUF2544 domain‑containing protein YfcO
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/18);  total (14/18)

CATTGACTACTTTTATATTCACCAACGTACACTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCTTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCA  >  NC_000913/2449865‑2450131                                                                                                                                    |                                                                                                                                        cATTGACTACTTTTATATTCACCAACGTACACTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTgtagt                                                                                                                                 <  2:744219/140‑1 (MQ=255)   ttGACTACTTTTATATTCACCAACGTACACTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTgtag                                                                                                                                  <  2:417016/137‑1 (MQ=255)   ttGACTACTTTTATATTCACCAACGTACACTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTgtag                                                                                                                                  >  1:417016/1‑137 (MQ=255)             ttATATTCACCAACGTACACTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCAt                                                                                                                     >  1:645188/1‑140 (MQ=255)             ttATATTCACCAACGTACACTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCAt                                                                                                                     >  1:766097/1‑140 (MQ=255)                tatTCACCAACGTACACTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGAc                                                                                                                  <  1:601357/140‑1 (MQ=255)                         aCGTACACTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGaa                                                                                                         <  1:99432/140‑1 (MQ=255)                          cGTACACTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGaac                                                                                                        <  2:390729/140‑1 (MQ=255)                             acacTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGaacaac                                                                                                     <  1:293741/140‑1 (MQ=255)                              cacTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGa                                                                                                          <  1:386327/134‑1 (MQ=255)                              cacTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGa                                                                                                          >  2:386327/1‑134 (MQ=255)                                 tGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATc                                                                                                 >  2:808987/1‑140 (MQ=255)                                 tGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATc                                                                                                 <  2:732383/140‑1 (MQ=255)                                      aTGGCATCATTATTTTTAAATTAAAAGTCGCATTATTACGCGAATCACCGTTGTAACCGTAGTAAACGTAAACGCTATGCCCGGCAGCCATAGCTTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCgtagt                                                                                            >  2:623038/1‑140 (MQ=255)                                           atcatTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGaacaac                                                                                                     <  2:588233/126‑1 (MQ=255)                                           atcatTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGaacaac                                                                                                     <  2:31773/126‑1 (MQ=255)                                           atcatTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGaacaac                                                                                                     >  1:588233/1‑126 (MQ=255)                                           atcatTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGaacaac                                                                                                     >  1:31773/1‑126 (MQ=255)                                                tattTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCAttt                                                                                  <  2:31697/140‑1 (MQ=255)                                                     ttaaattaaaAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAat                                                                             <  1:808987/140‑1 (MQ=255)                                                                   gCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCg                                                               <  2:435636/140‑1 (MQ=255)                                                                           aCGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTaaa                                                       >  2:960883/1‑140 (MQ=255)                                                                            cgcgAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAag                                                      >  1:92697/1‑140 (MQ=255)                                                                                aaTCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGc                                                  <  1:960883/140‑1 (MQ=255)                                                                                 aTCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCg                                                 >  2:767022/1‑140 (MQ=255)                                                                                                         gTCAACGCTATGCCCGGCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCa                         >  2:362698/1‑140 (MQ=255)                                                                                                                        ggCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGcc          <  2:92697/140‑1 (MQ=255)                                                                                                                        ggCAGCCATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCGGTAAAATCACCGAACGCCACAACACCGCCGGtttttcgc          >  2:744634/1‑136 (MQ=255)                                                                                                                             ccATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGcc          <  1:8323/135‑1 (MQ=255)                                                                                                                             ccATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGcc          >  2:8323/1‑135 (MQ=255)                                                                                                                              cATAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGc    <  1:767022/140‑1 (MQ=255)                                                                                                                               aTAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGcc   >  2:907651/1‑140 (MQ=255)                                                                                                                                tAGCCTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCa  <  2:515588/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        CATTGACTACTTTTATATTCACCAACGTACACTGCGGATGGCATCATTATTTTTAAATTAAAAGTCGCAATATTACGCGAATCACCGTTGTAACCGTAGTAAACGTCAACGCTATGCCCGGCAGCCATAGCTTCTGTAGTAGCAATAAGCATGACCAGAGCGAACAACCATCGTAGTATTTTCATTTTAAATATACTCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCA  >  NC_000913/2449865‑2450131

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