Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,472,706	A→G	intergenic (‑173/+108)	pyrL ← / ← yjgH	pyrBI operon leader peptide/RutC family protein YjgH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,472,706	0	A	G	100.0%	205.7 / NA	57	intergenic (‑173/+108)	pyrL/yjgH	pyrBI operon leader peptide/RutC family protein YjgH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/30);  total (27/30)

TCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATGCGC  >  NC_000913/4472572‑4472832                                                                                                                                       |                                                                                                                               tCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTc                                                                                                                           >  1:146345/1‑140 (MQ=255)  cAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTcc                                                                                                                          >  1:751104/1‑140 (MQ=255)  cAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTcc                                                                                                                          <  1:746974/140‑1 (MQ=255)  cAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTcc                                                                                                                          >  1:638841/1‑140 (MQ=255)   aGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTccc                                                                                                                         >  2:319328/1‑140 (MQ=255)           cgcAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGcttttt                                                                                                                 >  2:1035394/1‑140 (MQ=255)           cgcAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGcttttt                                                                                                                 >  1:442675/1‑140 (MQ=255)                          cAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTtttcactttca                                                                                                  <  2:456227/140‑1 (MQ=255)                          cAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTtttcactttca                                                                                                  >  1:456227/1‑140 (MQ=255)                             tatTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGc                                                                                               <  1:1035394/140‑1 (MQ=255)                              atTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCa                                                                                              <  2:751104/140‑1 (MQ=255)                                  tGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCa                                                                                          >  1:435238/1‑140 (MQ=255)                                       gcgcAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCaa                                                                                             >  2:819683/1‑132 (MQ=255)                                       gcgcAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCaa                                                                                             <  1:819683/132‑1 (MQ=255)                                       gcgcAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGatat                                                                                     <  2:202596/140‑1 (MQ=255)                                         gcAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGc                                                                                   >  2:414829/1‑140 (MQ=255)                                                     tttGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCAt                                                                       >  2:539269/1‑140 (MQ=255)                                                            tAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCg                                                                <  1:980267/140‑1 (MQ=255)                                                              aaaGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGAc                                                              <  1:259505/140‑1 (MQ=255)                                                                 gTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTta                                                           <  1:780549/140‑1 (MQ=255)                                                                    aaGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAg                                                        >  1:846715/1‑140 (MQ=255)                                                                      ggTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAggg                                                      <  2:638841/140‑1 (MQ=255)                                                                      ggTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAggg                                                      <  2:435238/140‑1 (MQ=255)                                                                      ggTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAggg                                                      <  1:604911/140‑1 (MQ=255)                                                                        tataTGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCg                                                    <  2:2675/140‑1 (MQ=255)                                                                          tatGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGa                                                  <  1:248724/140‑1 (MQ=255)                                                                            tGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAgt                                                >  1:959329/1‑140 (MQ=255)                                                                             gCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAgtg                                               <  2:105365/140‑1 (MQ=255)                                                                                tttttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGAc                                                              <  2:148017/122‑1 (MQ=255)                                                                                tttttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGAc                                                              >  1:148017/1‑122 (MQ=255)                                                                                tttttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGa                                            >  2:736502/1‑140 (MQ=255)                                                                                 ttttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGaa                                           >  1:602211/1‑140 (MQ=255)                                                                                 ttttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGaa                                           >  1:211858/1‑140 (MQ=255)                                                                                  tttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGa                                                  >  1:214715/1‑132 (MQ=255)                                                                                  tttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGa                                                  <  2:214715/132‑1 (MQ=255)                                                                                  tttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGaaa                                          >  2:495906/1‑140 (MQ=255)                                                                                  tttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGaaa                                          >  2:1077871/1‑140 (MQ=255)                                                                                  tttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGaaa                                          >  2:312699/1‑140 (MQ=255)                                                                                      tttGATGTGATTCTGCAGGGGATCTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGa                                      <  2:557362/140‑1 (MQ=255)                                                                                      tttGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGa                                      <  1:414829/140‑1 (MQ=255)                                                                                       ttGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGaa                                     >  1:539946/1‑140 (MQ=255)                                                                                        tGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAAc                                    <  1:1077871/140‑1 (MQ=255)                                                                                         gATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACg                                   <  1:312699/140‑1 (MQ=255)                                                                                           tgtgATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGc                                 >  2:770190/1‑140 (MQ=255)                                                                                              gATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCACCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTaa                              <  2:959329/140‑1 (MQ=255)                                                                                                ttCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAaca                            <  2:602211/140‑1 (MQ=255)                                                                                                ttCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAaca                            <  2:846715/140‑1 (MQ=255)                                                                                                     cAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTAtt                       >  2:592176/1‑140 (MQ=255)                                                                                                     cAGGGGAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTAtt                       >  1:726461/1‑140 (MQ=255)                                                                                                          gAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                  <  2:442675/140‑1 (MQ=255)                                                                                                          gAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                  <  2:211858/140‑1 (MQ=255)                                                                                                          gAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                  <  1:826237/140‑1 (MQ=255)                                                                                                          gAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                  <  1:770190/140‑1 (MQ=255)                                                                                                          gAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                  <  1:736502/140‑1 (MQ=255)                                                                                                          gAACTGTCCTTCGGTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                  >  1:725977/1‑140 (MQ=255)                                                                                                                       gTATCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATg     <  2:38979/140‑1 (MQ=255)                                                                                                                          tCAATAATTGCAGGCATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATgcgc  >  1:389019/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                               TCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATGCGC  >  NC_000913/4472572‑4472832

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