Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,147	G→A	intergenic (‑22/+72)	yjcO ← / ← fdhF	Sel1 repeat‑containing protein YjcO/formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,147	0	G	A	100.0%	158.3 / NA	45	intergenic (‑22/+72)	yjcO/fdhF	Sel1 repeat‑containing protein YjcO/formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/27);  total (18/27)

TTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTAC  >  NC_000913/4297012‑4297279                                                                                                                                        |                                                                                                                                     ttGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGt                                                                                                                                  <  1:239176/140‑1 (MQ=255) ttGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGt                                                                                                                                  <  2:260269/140‑1 (MQ=255)   gtgCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTaa                                                                                                                                <  2:899853/140‑1 (MQ=255)   gtgCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTaa                                                                                                                                <  1:9083/140‑1 (MQ=255)    tgCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTaa                                                                                                                                <  1:1138590/139‑1 (MQ=255)    tgCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTaa                                                                                                                                >  2:1138590/1‑139 (MQ=255)    tgCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAAc                                                                                                                               <  1:637877/140‑1 (MQ=255)     gcgcGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAAcc                                                                                                                              <  1:525495/140‑1 (MQ=255)     gcgcGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAAcc                                                                                                                              <  2:247145/140‑1 (MQ=255)     gcgcGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAAcc                                                                                                                              >  1:202927/1‑140 (MQ=255)              tCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAgaga                                                                                                                                    <  1:739852/125‑1 (MQ=255)              tCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAgaga                                                                                                                                    >  2:739852/1‑125 (MQ=255)                 ccGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGt                                                                                                                  >  2:597662/1‑140 (MQ=255)                        ctgctgCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATg                                                                                                           <  1:346141/140‑1 (MQ=255)                              ccTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAag                                                                                                                     >  1:339812/1‑124 (MQ=255)                              ccTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAag                                                                                                                     <  2:339812/124‑1 (MQ=255)                                  tAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATaaaaaaa                                                                                                 >  1:764411/1‑140 (MQ=255)                                        cTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCt                                                                                           >  1:952427/1‑140 (MQ=255)                                           gCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCg                                                                                        >  2:1059309/1‑140 (MQ=255)                                           gCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCg                                                                                        >  2:552727/1‑140 (MQ=255)                                           gCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCg                                                                                        >  2:633899/1‑140 (MQ=255)                                           gCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCg                                                                                        >  2:801483/1‑140 (MQ=255)                                           gCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCg                                                                                        >  2:953053/1‑140 (MQ=255)                                               ccAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAag                                                                                    <  1:923734/140‑1 (MQ=255)                                                cAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAagg                                                                                   <  1:663753/140‑1 (MQ=255)                                                   gCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAaggagg                                                                                <  2:379259/140‑1 (MQ=255)                                                   gCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAaggagg                                                                                <  1:597662/140‑1 (MQ=255)                                                    cTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGc                                                                               <  2:202927/140‑1 (MQ=255)                                                         aGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAg                                                                          >  2:531083/1‑140 (MQ=255)                                                          gTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGa                                                                         <  2:1020193/140‑1 (MQ=255)                                                          gTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGa                                                                         <  1:1059309/140‑1 (MQ=255)                                                             gTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAaggagg                                                                                >  2:1119880/1‑130 (MQ=255)                                                             gTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAaggagg                                                                                <  1:1119880/130‑1 (MQ=255)                                                             gTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAg                                                                      >  2:342388/1‑140 (MQ=255)                                                                   gTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGt                                                                >  1:392294/1‑140 (MQ=255)                                                                   gTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGt                                                                >  1:343663/1‑140 (MQ=255)                                                                     gggCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAt                                                              >  2:311607/1‑140 (MQ=255)                                                                     gggCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAt                                                              <  1:311607/140‑1 (MQ=255)                                                                          aaGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGc                                                         <  1:633899/140‑1 (MQ=255)                                                                                                 tGCGACAATTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAg                                  <  1:289695/140‑1 (MQ=255)                                                                                                       aaTTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCa                            <  1:552727/140‑1 (MQ=255)                                                                                                       aaTTTTTTTCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCa                            <  2:764411/140‑1 (MQ=255)                                                                                                            ttttCATTTTGATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttg                       <  1:378820/140‑1 (MQ=255)                                                                                                                      gATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGt             <  2:1169568/140‑1 (MQ=255)                                                                                                                      gATAATACGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGt             <  2:1056808/140‑1 (MQ=255)                                                                                                                          ataCGATTGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAAt         <  1:531083/140‑1 (MQ=255)                                                                                                                                 tGGGTAGAGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTAc  >  2:345412/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                     TTGTGCGCGTCGGTCCCCGGCCTCTGCTGCCTTTAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTAC  >  NC_000913/4297012‑4297279

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: ‑