Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	527,577	T→C	intergenic (+36/‑4)	rhsD → / → ybbC	protein RhsD/PF15631 family protein YbbC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	527,577	0	T	C	100.0%	186.4 / NA	54	intergenic (+36/‑4)	rhsD/ybbC	protein RhsD/PF15631 family protein YbbC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/24);  total (30/24)

CCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGATTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATACT  >  NC_000913/527440‑527713                                                                                                                                          |                                                                                                                                         ccTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTAGGGCCTATAAGGAGACTc                                                                                                                                        >  2:101062/1‑140 (MQ=255)  cTGATGGACGATTTGAGGCGCTTAATAATAAAGGACAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCa                                                                                                                                       <  2:571691/140‑1 (MQ=255)     aTGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATg                                                                                                                                    <  2:321973/140‑1 (MQ=255)          cGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAata                                                                                                                               <  2:542547/140‑1 (MQ=255)                      cTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATAtt                                                                                                                   >  2:586272/1‑140 (MQ=255)                       ttAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATAttt                                                                                                                  >  2:566619/1‑140 (MQ=255)                                 aGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCttt                                                                                                        >  1:14552/1‑140 (MQ=255)                                 aGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCttt                                                                                                        >  1:655675/1‑140 (MQ=255)                                  ggAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTc                                                                                                       <  1:566619/140‑1 (MQ=255)                                  ggAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTc                                                                                                       <  2:353613/140‑1 (MQ=255)                                     gaATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAtt                                                                                                    >  1:646325/1‑140 (MQ=255)                                      aatgtatTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAttt                                                                                                   >  2:75815/5‑140 (MQ=255)                                                  tCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTAttt                                                                                       >  2:38510/1‑140 (MQ=255)                                                   cAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTg                                                                                      <  2:53057/140‑1 (MQ=255)                                                     gAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGcc                                                                                    <  1:729966/140‑1 (MQ=255)                                                              tCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATgag                                                                           >  2:189558/1‑140 (MQ=255)                                                               cccGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTAGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATgaga                                                                          <  1:101062/140‑1 (MQ=255)                                                                ccGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGAc                                                                         >  2:477555/1‑140 (MQ=255)                                                                ccGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGAc                                                                         >  1:79341/1‑140 (MQ=255)                                                                 cGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACa                                                                        <  2:447057/140‑1 (MQ=255)                                                                       tGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGtt                                                                  <  2:655675/140‑1 (MQ=255)                                                                        gAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGttt                                                                 >  1:493502/1‑140 (MQ=255)                                                                        gAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGttt                                                                 >  2:137312/1‑140 (MQ=255)                                                                              aGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGCGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGtt                                                           >  1:42469/1‑140 (MQ=255)                                                                              aGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGtt                                                           >  2:434425/1‑140 (MQ=255)                                                                              aGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGtt                                                           <  1:137312/140‑1 (MQ=255)                                                                              aGTATGCATACCGTAGCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGtt                                                           <  2:517094/140‑1 (MQ=255)                                                                                        ccGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGaaaa                                                 >  2:213859/1‑140 (MQ=255)                                                                                          gTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACa                                               >  1:804049/1‑140 (MQ=255)                                                                                           tATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTTCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACat                                              >  2:107539/1‑140 (MQ=255)                                                                                            aTCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACatt                                             >  1:24020/1‑140 (MQ=255)                                                                                              cATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACattat                                           >  2:15551/1‑140 (MQ=255)                                                                                              cATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACattat                                           >  1:9066/1‑140 (MQ=255)                                                                                                tGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTAttt                                         >  2:17469/1‑140 (MQ=255)                                                                                                     aaTTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATg                                    >  2:671918/1‑140 (MQ=255)                                                                                                           aGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTAt                              >  1:515532/1‑140 (MQ=255)                                                                                                             tGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTAtgt                            <  1:244684/140‑1 (MQ=255)                                                                                                                           gggCCTATAAGGTGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              <  2:42469/140‑1 (MQ=255)                                                                                                                           gggCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              >  1:496657/1‑140 (MQ=255)                                                                                                                           gggCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              <  1:189558/140‑1 (MQ=255)                                                                                                                           gggCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              <  1:221394/140‑1 (MQ=255)                                                                                                                            ggCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAttt             >  1:459787/1‑140 (MQ=255)                                                                                                                              ccTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              >  1:517562/1‑137 (MQ=255)                                                                                                                              ccTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              <  2:517562/137‑1 (MQ=255)                                                                                                                              ccTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAg           >  1:624625/1‑140 (MQ=255)                                                                                                                              ccTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAg           >  2:213837/1‑140 (MQ=255)                                                                                                                                  taAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAAATaa       >  2:827074/1‑140 (MQ=255)                                                                                                                                     ggAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATa    <  1:17469/140‑1 (MQ=255)                                                                                                                                       agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  1:38510/140‑1 (MQ=255)                                                                                                                                       agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  1:434425/140‑1 (MQ=255)                                                                                                                                       agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  1:213859/140‑1 (MQ=255)                                                                                                                                       agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  2:493502/140‑1 (MQ=255)                                                                                                                                       agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  2:459787/140‑1 (MQ=255)                                                                                                                                       agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  1:477555/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         CCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGATTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATACT  >  NC_000913/527440‑527713

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