Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,525,930	C→T	*309* (TAG→TAA)	yfeR ←	transcriptional regulator of yefH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,525,930	0	C	T	100.0%	130.1 / NA	41	*309* (TAG→TAA)	yfeR	transcriptional regulator of yefH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/16);  total (25/16)

ACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAA  >  NC_000913/2525794‑2526023                                                                                                                                         |                                                                                              caacatcacacattactacggggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTta                                                                                             >  1:997616‑M2/133‑139 (MQ=255)                 tacggggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGa                                                                             >  1:32619‑M2/117‑139 (MQ=255)                     gggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATccc                                                                         >  1:890680‑M2/113‑139 (MQ=255)                        attcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCtt                                                                      >  1:115272‑M2/110‑139 (MQ=255)                        attcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCtt                                                                      >  1:708217‑M2/110‑139 (MQ=255)                              tatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGcc                                                                >  1:112219‑M2/104‑139 (MQ=255)                                        cggatgacaaactgggcgaagatttttacccggattcgataatctactccatcaaatttaacttctgatagattactgccgaggctcccgacaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCATCCATAAg                                                      >  2:429405‑M2/94‑139 (MQ=255)                                        cggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAg                                                      <  2:32619‑M2/46‑1 (MQ=255)                                          gatgacatactgggcgaagatttttaccaggatgcgataatctactccatctaatttaacttctgattgattactgacggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCTCCCAGCCATTAgcg                                                    >  1:680741‑M2/92‑139 (MQ=255)                                          gatgacaaactgggcgaggatttttaccaggatgcgataatctactccatcaaatttaacttctgattgatgactgccgggactcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAgcg                                                    >  2:636004‑M2/92‑139 (MQ=255)                                          gatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAgcg                                                    >  2:232680‑M2/92‑139 (MQ=255)                                                  actgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                            <  1:423859‑M2/56‑1 (MQ=255)                                                  actgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                            >  2:1068676‑M2/84‑139 (MQ=255)                                                  actgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                            <  2:708217‑M2/56‑1 (MQ=255)                                                  actgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                            <  2:997616‑M2/56‑1 (MQ=255)                                                       gcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATc                                                                                         >  1:99929‑M2/79‑89 (MQ=255)                                                       gcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATc                                                                                         <  2:99929‑M2/11‑1 (MQ=255)                                                          aagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGAt                                                                                     >  2:445088‑M2/76‑90 (MQ=255)                                                          aagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGAt                                                                                     <  1:445088‑M2/15‑1 (MQ=255)                                                               ttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCa                               <  1:232680‑M2/69‑1 (MQ=255)                                                                 ttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACg                             <  2:115272‑M2/71‑1 (MQ=255)                                                                       ggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGcc                                                                <  2:68863‑M2/36‑1 (MQ=255)                                                                       ggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGcc                                                                >  1:68863‑M2/63‑98 (MQ=255)                                                                         atgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTCATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCAc                              >  1:1090007‑M2/61‑130 (MQ=255)                                                                         atgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTCATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCAc                              <  2:1090007‑M2/70‑1 (MQ=255)                                                                              ataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGAACCAGGTCCCATAAAGCtt                >  2:476336‑M2/56‑139 (MQ=255)                                                                                  tctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCTTTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgc            >  1:849394‑M2/52‑139 (MQ=255)                                                                                         catcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCa                                                  >  2:227617‑M2/45‑94 (MQ=255)                                                                                         catcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCa                                                  <  1:227617‑M2/50‑1 (MQ=255)                                                                                         catcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTg     >  2:44022‑M2/45‑139 (MQ=255)                                                                                           tcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcg           <  1:522123‑M2/89‑1 (MQ=255)                                                                                           tcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcg           >  2:522123‑M2/43‑131 (MQ=255)                                                                                            caaatttgacttctgattgattactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCGCGCCCAGCCATTAGCGCATTGCCCTGCTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaa  >  1:679144‑M2/42‑139 (MQ=255)                                                                                                                ttactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcg           >  2:321386‑M2/22‑110 (MQ=255)                                                                                                                ttactgccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcg           <  1:321386‑M2/89‑1 (MQ=255)                                                                                                                     gccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCAc                              <  2:759706‑M2/70‑1 (MQ=255)                                                                                                                     gccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCAc                              >  1:759706‑M2/17‑86 (MQ=255)                                                                                                                     gccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCAc                              <  2:739258‑M2/70‑1 (MQ=255)                                                                                                                     gccggggctcccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCAc                              >  1:739258‑M2/17‑86 (MQ=255)                                                                                                                              cccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCa                                                  >  2:1098158‑M2/8‑57 (MQ=255)                                                                                                                              cccggcaCACTTTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCa                                                  <  1:1098158‑M2/50‑1 (MQ=255)                                                                                                                                         |                                                                                              ACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAA  >  NC_000913/2525794‑2526023

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