Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	27,781	G→A	E163E (GAG→GAA)	rihC →	ribonucleoside hydrolase RihC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	27,781	0	G	A	100.0%	75.1 / NA	26	E163E (GAG→GAA)	rihC	ribonucleoside hydrolase RihC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/13);  total (13/13)

GTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCTGTTACCCTGGTGGCCATCGGCCCGTTAACCAATATTGCGCTGTTACTTTCACAATGCCCGGAATGCAAGCCGTATATTCGCCGTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCG  >  NC_000913/27566‑28024                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    gttgAGCACAACCGAAAGCCGCTCGGGATACCGGCGTATCTAGCGATTCGGGCTGCCCTGATGCGTGCACCAGGGCCTGTTACCCTGGTGGCCAACGGCCCCTTCCCCCATATTCCGCTGTTACTTTCACAACACCCGGAATGCCAGCCGTATATTCAACGTCTCAAGATCATGGGCGGTTCTGCCGGACGCGGCAACTGTACGCCCACCGCCGAAATTAATATTGATGACGATCCAGAAg                                                                                                                                                                                                                            >  1:249942/1‑241 (MQ=255)                                                          cTGATGCGTGCACCAGAGCCTGTTACCCTGGTGGCCATCGGCCCGTTAACCAATATTGCGATGTTCCTTTCACAATGCCCGGAATGCAAGCCGTATATTCGCCGTATGGTGATCATGGGTGGTTCTGCCGGCAGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCaa                                                                                                                                                               >  1:290608/1‑245 (MQ=255)                                                            gATGCGTGCACCAGAGCCTGTTACCCTGGTGGCCATCGGCCCGTTAACCAATATTGCGCTGTTACTTTCACAATGCCCGGAATGCAAGCCGTATATTCGCCGTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGGATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGc                                                                                                                                                         <  1:470656/249‑1 (MQ=255)                                                                                   cccTGGTGGCCATCGGCCCGTTAACCAATATTGCGCTGTTACTTTCACAATGTCCGGAATGCAAGCCGTATATTCGCCGTCTGGTGCTCATGGGTGGTTCTGCCGGACGCGGCAACTGTAAGCCAAACGCCGAATTTAATATTGCTGCCGCTCCAGAAGCTGCTGCCTGTGTCTTCCGc                                                                                                                                                                                                        >  1:350618/1‑179 (MQ=255)                                                                                   cccTGGTGGCCATCGGCCCGGTAACCAATATTGCGCTGGTACTTTCACAATGCCCGGAATGCAAGCCGTATATTCGCCGTCTGGTGATCANGGGNGGTTNTGNCGGACGCGGCAACTGNACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCNNCNNNTNTCTTCCGc                                                                                                                                                                                                        <  2:350618/179‑1 (MQ=255)                                                                                    ccTGGTGGCCATCGGCCCGTTAACCAATATTGCGCTGTTACTTTCACAATGCCCGGAATGCAAGCCGTATATTCGCCGTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAANTTANNNNngcnn                                                                                                                                                                                                                                        >  2:299075/1‑144 (MQ=255)                                                                                            ccATCGGCCCGTTAACCAATATTGCGCTGTTACTTTCACAATGCCCGGAATGCCAGCCGTATATTCGCCGTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAAATGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCCATCAgg                                                                                                                                                          >  1:202793/1‑216 (MQ=255)                                                                                            ccATCGGCCCGTTAACCAATAGTGCGCTGTTACTGTCACAATGCCCGGAATGCAAGCCGTATATTCGCCGGCTGGTGATCATGGGTGGTTCTGCCGGACGCGTCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGctgctg                                                                                                                                                                                                                      <  1:245210/156‑1 (MQ=255)                                                                                                                      ntGTTACTTTCACAATGCCCGGAATGCAAGCCGTATATTCGCCGTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGcc                                                                                                                                                                                                                                      >  1:299683/2‑114 (MQ=255)                                                                                                                      cTGTTACTTTCACAATGCCCGGAATGCAAGCCGTATATTCGCCGTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGcc                                                                                                                                                                                                                                      <  2:299683/114‑1 (MQ=255)                                                                                                                      cTGTTACGTTCACAATGCCCGGAATGCAAGCCGTATATTCGCCGTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTGCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTc                                                                                                                                              <  1:204000/202‑1 (MQ=255)                                                                                                                        gTTACGTTCACAATGGCCGGAATGCAAGCCGTATATTCGCCTTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATGGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCttn                                                                                           <  1:304490/251‑2 (MQ=255)                                                                                                                                                      gTATATTCGCCGTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCANNAgctg                                                                                                                                                                                                                         >  2:4528/1‑95 (MQ=255)                                                                                                                                                      gTATATTCGCCGTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAgctg                                                                                                                                                                                                                         <  1:4528/95‑1 (MQ=255)                                                                                                                                                                 gTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCTTAGCGGCAGTATGCAAAGc                                                  <  1:385734/251‑1 (MQ=255)                                                                                                                                                                               ggtggtTCTGCCGGACGCGGCAACTGTACGCCAAACGCCCAATTTAATATTGCTGCCGATCCAGACGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATCTTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCAc                                                                                         >  1:376315/1‑198 (MQ=255)                                                                                                                                                                                   gtTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTg                                                                                                                                                                                                    >  1:264067/1‑87 (MQ=255)                                                                                                                                                                                   gtTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAAGGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATgn                                                       <  1:52907/228‑2 (MQ=255)                                                                                                                                                                                         ccGGACGCGGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCtgctgc                                                                                                                                                                                                                     <  1:291164/64‑1 (MQ=255)                                                                                                                                                                                              cgcgGCAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGc                                                                                                                      >  1:464535/1‑154 (MQ=255)                                                                                                                                                                                                   cAACTGTACGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCct                >  1:433050/1‑251 (MQ=255)                                                                                                                                                                                                           ngCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGtt                                                                                 >  1:52442/2‑178 (MQ=255)                                                                                                                                                                                                           cGCCAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCGCTGCGCCATCGCCTGGCTGGTg        <  1:106738/251‑1 (MQ=255)                                                                                                                                                                                                             ncAAACGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCCATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTgcg      >  1:45338/2‑251 (MQ=255)                                                                                                                                                                                                               aaaCGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCGTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTGGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCgn                                                                                                                      <  1:48127/137‑2 (MQ=255)                                                                                                                                                                                                                 aCGCCGAATTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGcccn  <  1:299075/251‑2 (MQ=255)                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    GTTGAGCACAACCGAAAGCCGCTCGGGATACCGGCGTTTCTGGCGATTCGGGATGCCCTGATGCGTGCACCAGAGCCTGTTACCCTGGTGGCCATCGGCCCGTTAACCAATATTGCGCTGTTACTTTCACAATGCCCGGAATGCAAGCCGTATATTCGCCGTCTGGTGATCATGGGTGGTTCTGCCGGACGCGGCAACTGTACGCCAAACGCCGAGTTTAATATTGCTGCCGATCCAGAAGCTGCTGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCG  >  NC_000913/27566‑28024

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