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

ATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGATTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATACT  >  NC_000913/527444‑527713                                                                                                                                      |                                                                                                                                         aTGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATg                                                                                                                                    >  1:117528/1‑140 (MQ=255)                 gcTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAatat                                                                                                                    >  1:360087/1‑140 (MQ=255)                               gagaATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCa                                                                                                      >  1:210965/1‑140 (MQ=255)                                agaATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTGATTAGTGGGACCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAt                                                                                                     >  2:353901/1‑140 (MQ=255)                                    tatatTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAttttt                                                                                                 >  2:248690/1‑140 (MQ=255)                                         ttttATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATAc                                                                                            >  1:30259/1‑140 (MQ=255)                                            tATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTAt                                                                                         >  1:41822/1‑140 (MQ=255)                                             aTCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTAtt                                                                                        <  1:262917/140‑1 (MQ=255)                                                  aaaTAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGAGGCTTTCATTTTTTATACTATTTGCCt                                                                                   >  1:277423/1‑140 (MQ=255)                                                                     aGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTa                                                                <  1:39887/140‑1 (MQ=255)                                                                     aGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTa                                                                <  2:277432/140‑1 (MQ=255)                                                                       ttCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACg                                                              <  2:30259/140‑1 (MQ=255)                                                                          aGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGtt                                                           <  2:106376/140‑1 (MQ=255)                                                                                    ccGTCTCATGATTAATTTTAGTGCTTTTATTAATGGTGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGaaaa                                                 >  2:346430/1‑140 (MQ=255)                                                                                    ccGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGaaaa                                                 >  1:164186/1‑140 (MQ=255)                                                                                      gTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACa                                               >  2:39342/1‑140 (MQ=255)                                                                                         tCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACatta                                            >  2:41063/1‑140 (MQ=255)                                                                                         tCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACatta                                            <  1:41063/140‑1 (MQ=255)                                                                                         tCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACatta                                            <  1:353901/140‑1 (MQ=255)                                                                                                tAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTAt                                     <  2:117528/140‑1 (MQ=255)                                                                                                 aaTTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATg                                    <  2:3987/140‑1 (MQ=255)                                                                                                   ttttAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAg                                  <  1:248690/140‑1 (MQ=255)                                                                                                                       gggCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              <  1:113146/140‑1 (MQ=255)                                                                                                                          ccTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              <  2:128086/137‑1 (MQ=255)                                                                                                                          ccTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              >  1:128086/1‑137 (MQ=255)                                                                                                                                   agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  1:346430/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         ATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGATTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATACT  >  NC_000913/527444‑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: ‑