Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	27,519	Δ1 bp	coding (227/915 nt)	rihC →	ribonucleoside hydrolase RihC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	27,519	0	A	.	100.0%	225.6 / NA	54	coding (227/915 nt)	rihC	ribonucleoside hydrolase RihC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (28/26);  total (28/26)

AACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCTGTTACCCTGG  >  NC_000913/27384‑27653                                                                                                                                        |                                                                                                                                       aaCTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCgg                                                                                                                                   <  1:73945/140‑1 (MQ=255) aaCTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCgg                                                                                                                                   <  1:470079/140‑1 (MQ=255)    tgatgaCCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCAt                                                                                                                                <  1:937444/140‑1 (MQ=255)    tgatgaCCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCAt                                                                                                                                <  2:540469/140‑1 (MQ=255)      atgaCCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCt                                                                                                                              >  2:340561/1‑140 (MQ=255)       tgaCCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCtg                                                                                                                             >  1:155486/1‑140 (MQ=255)           caccGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCa                                                                                                                         >  2:184960/1‑140 (MQ=255)              cGTCGCGGATAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACgg                                                                                                                      >  2:203354/1‑140 (MQ=255)                       tAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCggg                                                                                                             <  1:771814/140‑1 (MQ=255)                                  ttGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATggctggct                                                                                                  >  2:969700/1‑140 (MQ=255)                                   tGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTa                                                                                                 <  2:420115/140‑1 (MQ=255)                                     agaAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACg                                                                                               >  1:1133277/1‑140 (MQ=255)                                              cccGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACt                                                                                            >  1:957852/1‑134 (MQ=255)                                              cccGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACt                                                                                            <  2:957852/134‑1 (MQ=255)                                               ccGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAAtggc                                                                                                       >  1:578902/1‑122 (MQ=255)                                               ccGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAAtggc                                                                                                       <  2:578902/122‑1 (MQ=255)                                                  cAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGca                                                                                  <  1:969700/140‑1 (MQ=255)                                                      gCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAAc                                                                              >  2:875614/1‑140 (MQ=255)                                                          tGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAg                                                                                    >  1:148979/1‑130 (MQ=255)                                                          tGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAg                                                                                    <  2:148979/130‑1 (MQ=255)                                                          tGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGaa                                                                          <  1:757364/140‑1 (MQ=255)                                                            cAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAg                                                                        <  1:129865/140‑1 (MQ=255)                                                               ctgctgCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCg                                                                     <  2:155486/140‑1 (MQ=255)                                                                 gctgcATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGca                                                                                  <  2:930980/125‑1 (MQ=255)                                                                 gctgcATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGca                                                                                  >  1:930980/1‑125 (MQ=255)                                                                  ctgcATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTc                                                                  >  2:174908/1‑140 (MQ=255)                                                                  ctgcATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTc                                                                  >  1:1077750/1‑140 (MQ=255)                                                                               aTGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt                                                     >  2:835633/1‑140 (MQ=255)                                                                               aTGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGt                                                     >  1:71419/1‑140 (MQ=255)                                                                                  cGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACTGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTc                                                  >  2:171636/1‑140 (MQ=255)                                                                                  cGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTc                                                  >  2:532657/1‑140 (MQ=255)                                                                                          ccGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGAtt                                          >  2:841857/1‑140 (MQ=255)                                                                                              tCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCggg                                      <  1:1053676/140‑1 (MQ=255)                                                                                              tCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCggg                                      <  2:106217/140‑1 (MQ=255)                                                                                              tCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCggg                                      <  1:184960/140‑1 (MQ=255)                                                                                                gCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGAt                                    <  1:875614/140‑1 (MQ=255)                                                                                                gCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGAt                                    <  2:1077750/140‑1 (MQ=255)                                                                                                 cccAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATg                                   >  2:293237/1‑140 (MQ=255)                                                                                                 cccAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATg                                   >  2:833975/1‑140 (MQ=255)                                                                                                     aGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCt                               <  1:532657/140‑1 (MQ=255)                                                                                                           cGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCg                         >  2:926795/1‑140 (MQ=255)                                                                                                           cGCTGTGCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCg                         <  1:390666/140‑1 (MQ=255)                                                                                                              tgtgCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGc                      >  1:966782/1‑140 (MQ=255)                                                                                                              tgtgCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGc                      >  1:843192/1‑140 (MQ=255)                                                                                                              tgtgCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGc                      <  1:833975/140‑1 (MQ=255)                                                                                                                tgCCACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCAc                    <  1:785956/140‑1 (MQ=255)                                                                                                                  ccACTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCa                  >  1:1103204/1‑140 (MQ=255)                                                                                                                    aCTGGTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCaga                <  1:293237/140‑1 (MQ=255)                                                                                                                        gTACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCt            >  2:426239/1‑140 (MQ=255)                                                                                                                         tACGCGCACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCTg           <  1:426239/140‑1 (MQ=255)                                                                                                                            gcgcACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCTGTTa        <  1:174908/140‑1 (MQ=255)                                                                                                                             cgcACCGCGTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCTGTTAc       >  1:26881/1‑140 (MQ=255)                                                                                                                                  cgcgTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCTGTTACCCtgg  >  2:724271/1‑140 (MQ=255)                                                                                                                                  cgcgTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCTGTTACCCtgg  >  2:813617/1‑140 (MQ=255)                                                                                                                                  cgcgTG‑TGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCTGTTACCCtgg  >  2:318597/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       AACTGATGACCACCGTCGCGGGTAATGTCTCGGTTGAGAAAACTACCCGCAATGCCCTGCAACTGCTGCATTTCTGGAATGCGGAGATTCCGCTCGCCCAAGGGGCCGCTGTGCCACTGGTACGCGCACCGCGTGATGCGGCATCTGTGCACGGCGAATCGGGAATGGCTGGCTACGACTTTGTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCTGTTACCCTGG  >  NC_000913/27384‑27653

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