Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,860,426	Δ1 bp	intergenic (+10/+41)	pphB → / ← ygbI	serine/threonine‑specific protein phosphatase 2/DeoR family putative transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,860,425	0	T	.	100.0%	238.7 / NA	56	intergenic (+9/+42)	pphB/ygbI	serine/threonine‑specific protein phosphatase 2/DeoR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (27/29);  total (27/29)

AATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATAATCTG  >  NC_000913/2860287‑2860556                                                                                                                                           |                                                                                                                                    aaTAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                     >  2:338947/1‑139 (MQ=255)  aTAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGctc                                                                                                                                    >  1:423436/1‑138 (MQ=255)  aTAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGctc                                                                                                                                    <  2:423436/139‑2 (MQ=255)    aaaCGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGctctt                                                                                                                                  <  1:1042736/139‑4 (MQ=255)          cgcTGATTGTTTTACATTTGGACATATGATGTTTGATAAAATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGCACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGa                                                                                                                           >  2:359590/1‑139 (MQ=255)                  attTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAg                                                                                                                   <  1:896719/139‑1 (MQ=255)                         atTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCg                                                                                                            >  1:1170665/1‑139 (MQ=255)                             ggACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGaa                                                                                                                          >  2:1122354/1‑121 (MQ=255)                             ggACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGaa                                                                                                                          <  1:1122354/121‑1 (MQ=255)                             ggACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                             >  2:545549/1‑134 (MQ=255)                             ggACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                             <  1:545549/134‑1 (MQ=255)                             ggACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCAGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAGCCCGgtgt                                                                                                        >  2:785204/1‑139 (MQ=255)                                    tgatgTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCa                                                                                                 <  1:338947/139‑1 (MQ=255)                                    tgatgTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCa                                                                                                 <  1:359590/139‑1 (MQ=255)                                    tgatgTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCa                                                                                                 <  2:1170665/139‑1 (MQ=255)                                                      aGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAg                                                                                                                   <  1:204664/103‑1 (MQ=255)                                                      aGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAg                                                                                                                   >  2:204664/1‑103 (MQ=255)                                                            tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCa                                                                                            <  2:705901/120‑1 (MQ=255)                                                            tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCa                                                                                            >  1:705901/1‑120 (MQ=255)                                                                     aGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAagag                                                                <  2:494054/139‑1 (MQ=255)                                                                      gATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAagaga                                                               <  1:514043/139‑1 (MQ=255)                                                                            tatTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCt                                                         >  1:218945/1‑139 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAct                                                                                                              >  2:1069572/1‑80 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                             <  2:84094/82‑1 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                             <  2:315416/82‑1 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                             >  1:315416/1‑82 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAccc                                                                                                             >  1:84094/1‑82 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                              <  1:1069572/81‑1 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                              >  2:809935/1‑81 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                              <  1:809935/81‑1 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTg                                                                                                   <  1:929934/92‑1 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTg                                                                                                   >  2:929934/1‑92 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCa                                                                                                 <  2:877004/94‑1 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCa                                                                                                 >  1:877004/1‑94 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                            >  1:356061/1‑131 (MQ=255)                                                                                 aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                            <  2:356061/131‑1 (MQ=255)                                                                                   aCCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCAc                                                  >  1:380851/1‑139 (MQ=255)                                                                                   aCCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCAc                                                  >  2:939536/1‑139 (MQ=255)                                                                                           gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcgcgc                                          <  1:79570/139‑1 (MQ=255)                                                                                           gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcgcgc                                          >  2:79570/1‑139 (MQ=255)                                                                                           gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcgcgc                                          >  2:761319/1‑139 (MQ=255)                                                                                                   gcgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCa                                  <  2:218945/139‑1 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgtg                                                                                                                            >  2:1182220/1‑48 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCtgtg                                                                                                                            <  1:1182220/48‑1 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGgtg                                                                                                         <  2:240738/67‑1 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGgtg                                                                                                         >  1:240738/1‑67 (MQ=255)                                                                                                    cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCAc                                 >  2:99264/1‑139 (MQ=255)                                                                                                       gcggcTGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCActct                              >  1:864968/1‑139 (MQ=255)                                                                                                            tGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                              <  2:501618/54‑1 (MQ=255)                                                                                                            tGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAAcc                                                                                                              >  1:501618/1‑54 (MQ=255)                                                                                                            tGTCATTTTATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGc                         <  1:761319/139‑1 (MQ=255)                                                                                                                  tttATAAAATAAAGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcg                   <  1:650008/139‑1 (MQ=255)                                                                                                                      taaaataaaGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt               <  2:864968/139‑1 (MQ=255)                                                                                                                           taaaGTAATCTCATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGta          <  2:380851/139‑1 (MQ=255)                                                                                                                                taatctcATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTaataat     <  1:99264/136‑1 (MQ=255)                                                                                                                                   tctcATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt               >  2:1238953/1‑126 (MQ=255)                                                                                                                                   tctcATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt               <  1:1238953/126‑1 (MQ=255)                                                                                                                                   tctcATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTaataat     >  2:181029/1‑136 (MQ=255)                                                                                                                                   tctcATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTaataat     <  1:181029/136‑1 (MQ=255)                                                                                                                                   tctcATGC‑TCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATAATCTg  >  1:4165/1‑139 (MQ=255)                                                                                                                                           |                                                                                                                                    AATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATAATCTG  >  NC_000913/2860287‑2860556

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