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%	166.6 / NA	47	intergenic (+36/‑4)	rhsD/ybbC	protein RhsD/PF15631 family protein YbbC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/19);  total (28/19)

ACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGATTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATACT  >  NC_000913/527448‑527713                                                                                                                                   |                                                                                                                                         aCGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAat                                                                                                                                >  1:426699/1‑140 (MQ=255)  cGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAata                                                                                                                               <  1:143744/140‑1 (MQ=255)    aTTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAAtata                                                                                                                             >  2:165422/1‑140 (MQ=255)             gcTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAatat                                                                                                                    >  2:302512/1‑140 (MQ=255)                         aGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCttt                                                                                                        >  1:46542/1‑140 (MQ=255)                           gagaATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCa                                                                                                      >  2:582557/1‑140 (MQ=255)                           gagaATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCa                                                                                                      >  1:327232/1‑140 (MQ=255)                              aaTATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAttt                                                                                                   >  1:27005/1‑140 (MQ=255)                              aaTATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAttt                                                                                                   >  1:677881/1‑140 (MQ=255)                              aaTATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAttt                                                                                                   >  2:269706/1‑140 (MQ=255)                                  tatTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTta                                                                                               <  1:173306/140‑1 (MQ=255)                                  tatTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTta                                                                                               >  2:56493/1‑140 (MQ=255)                                     ttttATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATAc                                                                                            <  2:46542/140‑1 (MQ=255)                                       ttATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTtata                                                                                             <  1:284697/138‑1 (MQ=255)                                       ttATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTtata                                                                                             >  2:284697/1‑138 (MQ=255)                                        tATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGCGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTAt                                                                                         >  2:191525/1‑140 (MQ=255)                                        tATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTAt                                                                                         >  1:310948/1‑140 (MQ=255)                                        tATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTAt                                                                                         >  2:402163/1‑140 (MQ=255)                                         aTCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTAtt                                                                                        <  1:582557/140‑1 (MQ=255)                                             gAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGcc                                                                                    <  2:426699/140‑1 (MQ=255)                                             gAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGcc                                                                                    <  2:515876/140‑1 (MQ=255)                                                      tCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGcc                                                                                    <  1:157104/131‑1 (MQ=255)                                                      tCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGcc                                                                                    >  2:157104/1‑131 (MQ=255)                                                        ccGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGAc                                                                         >  2:807608/1‑140 (MQ=255)                                                         cGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACa                                                                        <  1:165422/140‑1 (MQ=255)                                                          gggATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAg                                                                       <  2:275334/140‑1 (MQ=255)                                                          gggATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAg                                                                       <  1:807608/140‑1 (MQ=255)                                                                                  cGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAAc                                                >  1:142184/1‑140 (MQ=255)                                                                                   gTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACa                                               >  1:471738/1‑140 (MQ=255)                                                                                    tATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACat                                              <  1:269706/140‑1 (MQ=255)                                                                                     aTCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACatt                                             >  2:758466/1‑140 (MQ=255)                                                                                       cATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACattat                                           >  1:101736/1‑140 (MQ=255)                                                                                       cATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACattat                                           >  1:155114/1‑140 (MQ=255)                                                                                           aTTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTAttttt                                       >  1:695982/1‑140 (MQ=255)                                                                                            ttAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTa                                      >  2:455081/1‑140 (MQ=255)                                                                                             tAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTAt                                     <  2:421518/140‑1 (MQ=255)                                                                                              aaTTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATg                                    <  2:28728/140‑1 (MQ=255)                                                                                              aaTTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATg                                    <  1:455081/140‑1 (MQ=255)                                                                                                                    gggCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              <  1:758466/140‑1 (MQ=255)                                                                                                                     ggCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAttt             >  2:594007/1‑140 (MQ=255)                                                                                                                       ccTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAg           >  1:166165/1‑140 (MQ=255)                                                                                                                       ccTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAg           >  1:503886/1‑140 (MQ=255)                                                                                                                       ccTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAg           >  1:634480/1‑140 (MQ=255)                                                                                                                               gagaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAc   <  2:485028/140‑1 (MQ=255)                                                                                                                                agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  1:402163/140‑1 (MQ=255)                                                                                                                                agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  1:144311/140‑1 (MQ=255)                                                                                                                                agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  >  2:189338/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                         ACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATA‑CCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGATTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATACT  >  NC_000913/527448‑527713

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