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%	129.3 / NA	40	*309* (TAG→TAA)	yfeR	transcriptional regulator of yefH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/19);  total (21/19)

ATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCA  >  NC_000913/2525796‑2526048                                                                                                                                       |                                                                                                                        acatcacacattactacggggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTtata                                                                                                                     <  1:994169‑M2/9‑1 (MQ=255)   atcacacattactacggggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑c                                                                                                                  <  2:297933‑M2/11‑1 (MQ=255)     cacacattactacggggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTg                                                                                                                >  2:640606‑M2/127‑139 (MQ=255)       cacattactacggggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGAt                                                                                                              >  1:754331‑M2/125‑139 (MQ=255)            tactacggggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAAc                                                                                                         >  1:424151‑M2/120‑139 (MQ=255)               tacggggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGa                                                                                                      <  2:754331‑M2/23‑1 (MQ=255)               tacggggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGa                                                                                                      >  2:382041‑M2/117‑139 (MQ=255)                 cggggattcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacctctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATc                                                                                                    >  1:605088‑M2/115‑139 (MQ=255)                  tacggggattcagtatcgatatgcggacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATATCCTGATACAACGGATc                                                                                                    >  2:982465‑M2/114‑139 (MQ=255)                      attcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCtt                                                                                               >  2:148402‑M2/110‑139 (MQ=255)                      attcagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCtt                                                                                               >  1:884421‑M2/110‑139 (MQ=255)                         cagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTccc                                                                                            <  1:640606‑M2/33‑1 (MQ=255)                         cagtatcgatatgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTccc                                                                                            <  1:382041‑M2/33‑1 (MQ=255)                                    tgcggatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCAtt                                                                                 >  1:372518‑M2/96‑139 (MQ=255)                                        gatgacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAgcg                                                                             >  2:782354‑M2/92‑139 (MQ=255)                                          ttacaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCa                                                                           <  1:503421‑M2/50‑1 (MQ=255)                                            acaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATcc                                                                                                   <  1:844899‑M2/26‑1 (MQ=255)                                            acaaactgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcacggcaCACTTTATAT‑CTGATACAACGGATcc                                                                                                   >  2:844899‑M2/88‑113 (MQ=255)                                                actgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATc                                                                                                    >  1:1060592‑M2/84‑108 (MQ=255)                                                actgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATc                                                                                                    <  2:1060592‑M2/25‑1 (MQ=255)                                                actgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                                                     >  1:539861‑M2/84‑139 (MQ=255)                                                actgggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGccc                                                                     <  1:148402‑M2/56‑1 (MQ=255)                                                    ggcgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGt                                                                 <  1:782354‑M2/60‑1 (MQ=255)                                                      cgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCa                                                                           <  1:37028‑M2/50‑1 (MQ=255)                                                      cgaagatttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCa                                                                           >  2:37028‑M2/78‑127 (MQ=255)                                                             ttttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCa                                                        <  2:884421‑M2/69‑1 (MQ=255)                                                              tttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCAc                                                       <  2:539861‑M2/70‑1 (MQ=255)                                                               ttaccaggatgcgataatctactccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACg                                                      <  2:424151‑M2/71‑1 (MQ=255)                                                                                      ccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTg                                                                   >  2:732803‑M2/46‑103 (MQ=255)                                                                                      ccatcaaatttaacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTg                                                                   <  1:732803‑M2/58‑1 (MQ=255)                                                                                                 aacttctgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAAc                    >  1:54829‑M2/35‑139 (MQ=255)                                                                                                       tgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                >  2:1030791‑M2/29‑73 (MQ=255)                                                                                                       tgattgattactgccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                <  1:1030791‑M2/45‑1 (MQ=255)                                                                                                                   gccggggctcccggcaCACTTTATAT‑CTGATACAACg                                                                                                        <  2:791740‑M2/21‑1 (MQ=255)                                                                                                                   gccggggctcccggcaCACTTTATAT‑CTGATACAACg                                                                                                        >  1:791740‑M2/17‑37 (MQ=255)                                                                                                                   gccggggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCa  >  2:455423‑M2/17‑139 (MQ=255)                                                                                                                        ggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGa                      <  2:1035990‑M2/103‑1 (MQ=255)                                                                                                                        ggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGa                      >  1:1035990‑M2/12‑114 (MQ=255)                                                                                                                        ggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTg         <  2:795397‑M2/116‑1 (MQ=255)                                                                                                                        ggctcccggcaCACTTTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTg         >  1:795397‑M2/12‑127 (MQ=255)                                                                                                                                       |                                                                                                                        ATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATAT‑CTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCA  >  NC_000913/2525796‑2526048

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