Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007779	4,479,514	A→G	intergenic (+117/‑28)	yjgJ → / → yjgK	transcriptional regulator/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	4,479,514	0	A	G	100.0%	257.4 / NA	61	intergenic (+117/‑28)	yjgJ/yjgK	transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/26);  total (35/26)

GGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCA  >  NC_007779/4479420‑4479609                                                                                               |                                                                                                gggAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgt                                                                                           >  3:609143/1‑101 (MQ=255) gggAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCgttgt                                                                                           <  1:1214667/101‑1 (MQ=255)     aaGCTCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTtta                                                                                       <  1:360572/101‑1 (MQ=255)     aaGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTATCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTtta                                                                                       <  1:519856/101‑1 (MQ=255)      agggcCTGCCGATAAGTTGTGATAAGACAATAATTCTCGCGTTATGGCTAGCGGAATTGGTTATCTTTTCGTATAACGGTAGAAATGAAGCGTTGTTttaa                                                                                      <  1:1580078/98‑1 (MQ=255)         gcCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatta                                                                                   >  2:155718/1‑101 (MQ=255)         gcCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatta                                                                                   >  2:73440/1‑101 (MQ=255)         cccTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTttaatta                                                                                   <  3:718057/100‑1 (MQ=255)           cTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAg                                                                                 >  2:112861/1‑101 (MQ=255)             gCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGa                                                                               >  1:1484440/1‑101 (MQ=255)               cGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGt                                                                             <  1:641119/101‑1 (MQ=255)                     gTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                                                       >  3:914828/1‑101 (MQ=255)                     gTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGaaaa                                                                       >  3:1376043/1‑101 (MQ=255)                        gtgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgat                                                                    <  3:946272/101‑1 (MQ=255)                         tgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatg                                                                   >  1:582575/1‑101 (MQ=255)                         tgATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatg                                                                   <  3:120050/101‑1 (MQ=255)                          gATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatga                                                                  >  3:1146693/1‑101 (MQ=255)                          gATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAgatga                                                                  <  3:693382/101‑1 (MQ=255)                            tAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGatc                                                                >  3:1095668/1‑101 (MQ=255)                                    taatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaaa                                                        >  4:20852/1‑101 (MQ=255)                                    taatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGaaa                                                        >  1:981298/1‑101 (MQ=255)                                     aatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAat                                                       <  3:1355874/101‑1 (MQ=255)                                     aatTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAAAAGAACAACGGgaaa                                                       >  4:7787/1‑100 (MQ=255)                                     aatTCACGAGTTCCGGCTANCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAat                                                       >  4:91905/1‑101 (MQ=255)                                         cACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTc                                                   <  3:1174283/101‑1 (MQ=255)                                         cACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTc                                                   <  1:492633/101‑1 (MQ=255)                                          aCGCGTTAAGGTTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCa                                                  <  1:1074353/101‑1 (MQ=255)                                          aCGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCa                                                  >  1:155296/1‑101 (MQ=255)                                          aCGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCa                                                  >  1:1228391/1‑101 (MQ=255)                                           cgcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCAt                                                 <  1:984725/101‑1 (MQ=255)                                           cgcgTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCAt                                                 >  1:950286/1‑101 (MQ=255)                                              gTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAAt                                              >  3:288469/1‑101 (MQ=255)                                               ttAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATc                                             <  1:990170/101‑1 (MQ=255)                                                tAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCt                                            >  3:1296223/1‑101 (MQ=255)                                                   ggCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCa                                         <  1:316070/101‑1 (MQ=255)                                                   ggCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCa                                         <  1:1121557/101‑1 (MQ=255)                                                    gCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAg                                        <  3:1357482/101‑1 (MQ=255)                                                      tAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGcc                                      >  3:1532050/1‑101 (MQ=255)                                                      tAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGcc                                      <  3:463427/101‑1 (MQ=255)                                                         cGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTg                                   >  3:459171/1‑101 (MQ=255)                                                           gAATTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCTTGGc                                 >  3:1012829/1‑101 (MQ=255)                                                            aaTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCt                                <  1:1004712/101‑1 (MQ=255)                                                            aaTTGATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCt                                <  1:20253/101‑1 (MQ=255)                                                                gATTATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTAccc                            >  3:292915/1‑101 (MQ=255)                                                                  ttATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCa                          >  3:558596/1‑101 (MQ=255)                                                                   tATCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAg                         >  1:438591/1‑101 (MQ=255)                                                                     tCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCTTGGCTACCCCAGGa                       <  4:155589/101‑1 (MQ=255)                                                                     tCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGa                       <  3:1178490/101‑1 (MQ=255)                                                                     tCTTTTCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCGTCAGCCGGGGCTACCCCAGGa                       <  4:73335/101‑1 (MQ=255)                                                                        tttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGaaga                    >  1:1167296/1‑98 (MQ=255)                                                                        tttCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGtt                    >  3:1507440/1‑101 (MQ=255)                                                                          tCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTAc                  >  1:381270/1‑101 (MQ=255)                                                                          tCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTAc                  <  1:427431/101‑1 (MQ=255)                                                                          tCGTATAACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTAc                  >  1:1117759/1‑101 (MQ=255)                                                                               taACGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAg             >  1:1135911/1‑101 (MQ=255)                                                                                 aCGATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGc           >  1:853210/1‑101 (MQ=255)                                                                                   gATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGa         >  3:1376540/1‑101 (MQ=255)                                                                                   gATAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGCTACGCCAGGCGa         >  3:1313881/1‑101 (MQ=255)                                                                                     tAGAAATGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGAtt       >  3:275422/1‑101 (MQ=255)                                                                                         aaTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGc   >  3:499310/1‑101 (MQ=255)                                                                                          aTGAAGCGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAAGCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCa  <  3:1154876/101‑1 (MQ=255)                                                                                               |                                                                                                GGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATTGATTATCTTTTCGTATAACGATAGAAATGAAACGTTGTTTTAATTAAGGAGTGGAAAAGATGATCATCGGAAATATTCATAATCTTCAGCCGTGGCTACCCCAGGAGTTACGCCAGGCGATTGAGCA  >  NC_007779/4479420‑4479609

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