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

ACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGATTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATACT  >  NC_000913/527439‑527713                                                                                                                                           |                                                                                                                                         acCTGATGGACGATTTGAGGCGTTTAATAATTAAGGAGATTATATTTTCATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGAtt                                                                                                                                         <  2:272653/140‑1 (MQ=255)   cTGATGGATGATTTGAGCCGCTTAATAATAATGGAGAATGTATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCa                                                                                                                                       <  2:345407/140‑1 (MQ=255)     gATGGACGATTTGAGGCGCTTAATAGTAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAAt                                                                                                                                     <  1:56543/140‑1 (MQ=255)           cGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAata                                                                                                                               <  2:538416/140‑1 (MQ=255)           cGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAata                                                                                                                               >  2:571859/1‑140 (MQ=255)                    gcgcTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAat                                                                                                                      <  2:435020/140‑1 (MQ=255)                        ttAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATAttt                                                                                                                  >  2:152639/1‑140 (MQ=255)                           ataataaAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCg                                                                                                               <  1:152639/140‑1 (MQ=255)                                 aaGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGTGGCCTATAAGGAGACTAAATGAAATCTAGTACAATATTTTCGATGCtt                                                                                                         >  2:299460/1‑140 (MQ=255)                                  aGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGCGACTCAAAGAAGTATAGTTCAGTATTTTCGATGCttt                                                                                                        >  1:209251/1‑140 (MQ=255)                                  aGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCttt                                                                                                        >  1:168145/1‑140 (MQ=255)                                   ggAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTc                                                                                                       <  1:45448/140‑1 (MQ=255)                                   ggAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTc                                                                                                       <  1:571859/140‑1 (MQ=255)                                      gaATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAtt                                                                                                    <  2:168145/140‑1 (MQ=255)                                       aaTATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTc                                                                                                       >  1:462183/1‑136 (MQ=255)                                       aaTATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTc                                                                                                       <  2:462183/136‑1 (MQ=255)                                       aaTATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAttt                                                                                                   >  1:118216/1‑140 (MQ=255)                                       aaTATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAttt                                                                                                   >  1:456023/1‑140 (MQ=255)                                         tatatTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCAttttt                                                                                                 >  2:498161/1‑140 (MQ=255)                                              ttttATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATAc                                                                                            >  2:625906/1‑140 (MQ=255)                                                 tATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTAt                                                                                         >  2:168438/1‑140 (MQ=255)                                                      gAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGcc                                                                                    <  1:168438/140‑1 (MQ=255)                                                               tCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATgag                                                                           >  1:623168/1‑140 (MQ=255)                                                                 ccGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGTTGCTTTCGTTTATTGTACTATGTGCCTGTAATGAGAc                                                                         >  1:213482/1‑140 (MQ=255)                                                                   gggATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAg                                                                       <  2:118216/140‑1 (MQ=255)                                                                          aGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTa                                                                <  2:52873/140‑1 (MQ=255)                                                                               aGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGtt                                                           <  1:625906/140‑1 (MQ=255)                                                                                         ccGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGaaaa                                                 >  1:444087/1‑140 (MQ=255)                                                                                           gTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACa                                               >  1:727503/1‑140 (MQ=255)                                                                                           gTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACa                                               <  1:498161/140‑1 (MQ=255)                                                                                             aTCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACatt                                             >  1:642402/1‑140 (MQ=255)                                                                                                    ttAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTa                                      >  2:74600/1‑140 (MQ=255)                                                                                                    ttAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTa                                      >  1:710605/1‑140 (MQ=255)                                                                                                    ttAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTCCCTGTAATGAGAAAGCGGTTTACGGTTCTGTTGAAAACATTATTTTTa                                      >  2:224945/1‑140 (MQ=255)                                                                                                    ttAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTCCCTGTAATAAGACAGCGGTTTACGGATCTGGTGAAAACATTATTTTTa                                      >  2:223021/1‑140 (MQ=255)                                                                                                      aaTTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATg                                    <  2:642402/140‑1 (MQ=255)                                                                                                      aaTTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATg                                    <  2:71772/140‑1 (MQ=255)                                                                                                            aGTGCTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGGTGAAAACATAATTTTTATGAGGTAt                              >  1:251145/1‑140 (MQ=255)                                                                                                                cTTTTATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTgg                          <  2:710605/140‑1 (MQ=255)                                                                                                                  tttATTAGTGGGGCCTATAAGGAGACTCACTGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGaa                        <  1:299460/140‑1 (MQ=255)                                                                                                                  tttATTAGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGaa                        >  2:606936/1‑140 (MQ=255)                                                                                                                      ttaGTGGGGCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAAt                    <  2:444087/140‑1 (MQ=255)                                                                                                                            gggCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACAtt              <  1:74600/140‑1 (MQ=255)                                                                                                                             ggCCTATAAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAGACATTATATTTATGGGGTATGGGGAAAAATTACAttt             >  1:262631/1‑140 (MQ=255)                                                                                                                                   taAGGAGACTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATaa       <  1:606936/140‑1 (MQ=255)                                                                                                                                        agaCTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  2:456023/140‑1 (MQ=255)                                                                                                                                        agaCTCAATGAAATATACTTCACTATTTTCGCTGCATTCCTTTTTTATGCTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATAct  <  2:209251/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                         ACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTTTTATTAGTGGGGCCTATAAGGAGATTCAATGAAATATAGTTCAATATTTTCGATGCTTTCATTTTTTATACTATTTGCCTGTAATGAGACAGCTGTTTACGGTTCTGATGAAAACATTATTTTTATGAGGTATGTGGAAAAATTACATTTAGATAAATACT  >  NC_000913/527439‑527713

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