Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,088	+C	coding (2/624 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,088	1	.	C	100.0%	192.3 / NA	77	S133W (TCG→TGG)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (46/31);  total (46/31)

TGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCAC  >  NC_000913/3554013‑3554155                                                                             |                                                                     tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCg                                                                       >  1:2791283/1‑76 (MQ=255)   gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                     <  2:4749636/76‑3 (MQ=255)   gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                     <  1:1633737/76‑3 (MQ=255)   gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCgca                                                                     >  1:3466719/1‑74 (MQ=255)    cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                    >  2:3470276/1‑73 (MQ=255)       aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCGCCCGGCGGATAAAATCCGCAggcg                                                                 <  2:1489186/76‑1 (MQ=255)       aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                 <  2:3331533/76‑1 (MQ=255)       aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                 <  2:3579467/76‑1 (MQ=255)        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcgg                                                                <  2:4249639/76‑1 (MQ=255)         cggCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                               <  2:1375485/74‑1 (MQ=255)         cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                               <  1:883847/76‑1 (MQ=255)             gTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCg                                                           <  1:4227599/76‑1 (MQ=255)             gTGGTTGGCTGATGAATTCCTATTTCCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCg                                                           <  2:415880/76‑1 (MQ=255)              tGGTTTGCTGCTGCATTCCTATTTTCGCCAGCAGCGGCTCCCGCCCACAGCGGCTAAAACCCGCAGGCGGCACCgc                                                          >  1:1972899/1‑75 (MQ=255)               ggTTTGCTGATGAATTCCTATTTTCGCCGGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGgcggcggcggc                                                         >  1:939795/1‑74 (MQ=255)               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAc                                                         >  1:4248764/1‑76 (MQ=255)                gTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                        <  1:626923/76‑1 (MQ=255)                gTTTGCTGATGAATTCCTATTTTCGCAAGCAGCGGCTCCAGCACCCAGCGGATAAAATCCGCAGGCGGCGCCggtc                                                        >  1:1461057/1‑73 (MQ=255)                  ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                      >  1:1098895/1‑76 (MQ=255)                  ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                      >  1:1300062/1‑76 (MQ=255)                   tGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGt                                                     >  1:2296956/1‑76 (MQ=255)                   tGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGt                                                     >  1:3913858/1‑76 (MQ=255)                   tGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGt                                                     <  1:3895678/76‑1 (MQ=255)                    gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggggt                                                    >  1:3316020/1‑74 (MQ=255)                    gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGtt                                                    >  2:247615/1‑76 (MQ=255)                     cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTa                                                   >  1:4179155/1‑76 (MQ=255)                      tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTAt                                                  >  1:4262146/1‑76 (MQ=255)                         tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                               >  1:2840926/1‑76 (MQ=255)                          gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCggg                                              >  1:4091728/1‑75 (MQ=255)                          gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt                                              >  1:2714287/1‑76 (MQ=255)                          gaATTCCTATTGTCGCCAGCAGCGACTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt                                              >  1:3053425/1‑76 (MQ=255)                           aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCgggg                                             >  1:666993/1‑74 (MQ=255)                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                            >  1:215415/1‑74 (MQ=255)                               ccTATTTTCGCCAGCAGCGGCTCCAGTACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                       <  1:4394423/76‑1 (MQ=255)                                cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                      <  1:1807443/76‑2 (MQ=255)                                cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                      <  1:4604619/76‑2 (MQ=255)                                 tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccg                                     <  1:2747906/76‑3 (MQ=255)                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGTCGGCGCCGACGGGTTATCGGt‑gccg‑ccgc                                    >  1:3859406/1‑73 (MQ=255)                                      tCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                 <  1:4099206/76‑1 (MQ=255)                                      tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                 >  1:3617561/1‑76 (MQ=255)                                      tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCCACGGGGTATTGGTGGCGC‑CGCTccc                                 >  1:475139/1‑74 (MQ=255)                                       cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                <  1:2630716/76‑1 (MQ=255)                                       cGCCAGCAGCGGCTCCAGCACCCCGCGGAGAAAATCAGCAGGCGGCGCCCCCGGGTTATCGGTCCCCC‑CGCCCACt                                >  2:1718441/1‑76 (MQ=255)                                        cccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACtt                               <  1:700046/75‑1 (MQ=255)                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCa                             >  1:2950094/1‑76 (MQ=255)                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATAGGT‑GCCGCCGCTCACTTCaa                            >  1:4850711/1‑76 (MQ=255)                                             cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGGCGACGGGTTATCGGT‑GCCGCCGCTCCCTTCCacc                          >  2:4316987/1‑75 (MQ=255)                                             cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                          >  1:4030622/1‑76 (MQ=255)                                              agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCCCTTCAacac                         >  1:2387443/1‑76 (MQ=255)                                              agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGGCCCCC‑CGCCCCCCTCAaaac                         >  2:916362/1‑76 (MQ=255)                                               gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcg                        <  2:3019909/76‑1 (MQ=255)                                                 ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                      >  1:3458170/1‑76 (MQ=255)                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGa                     >  1:1516797/1‑76 (MQ=255)                                                    tCCAGCACCCAGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCCC‑AGCTCACTGCAACACGCCAAg                   >  2:3098578/1‑76 (MQ=255)                                                    cccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                   <  1:3094290/75‑1 (MQ=255)                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                 <  1:598005/76‑1 (MQ=255)                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                 <  1:247615/76‑1 (MQ=255)                                                       aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTc                >  1:3708937/1‑76 (MQ=255)                                                        ggaCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc               <  2:2430666/74‑1 (MQ=255)                                                        gCACCCGGCGGATAAAATCCGCGGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc               <  1:4245024/76‑1 (MQ=255)                                                        gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCCACACGCGAAGGTcc               >  1:1152671/1‑76 (MQ=255)                                                        gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc               <  2:2585498/76‑1 (MQ=255)                                                        gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc               <  1:2674116/76‑1 (MQ=255)                                                         cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa              >  2:2679099/1‑76 (MQ=255)                                                         cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGGCCa              >  2:3695984/1‑76 (MQ=255)                                                           cccGGCGGATAGAATCCGCGGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc            <  2:1149211/76‑1 (MQ=255)                                                           cccGGCGGATAAAATCCGCGGGCGGCACCGAGGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc            <  1:601504/76‑1 (MQ=255)                                                           cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc            <  2:1255810/76‑1 (MQ=255)                                                           cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc            <  1:4690778/76‑1 (MQ=255)                                                              ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc         >  1:706242/1‑76 (MQ=255)                                                               gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCa        >  1:3364904/1‑76 (MQ=255)                                                               gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCa        >  1:1905830/1‑76 (MQ=255)                                                               gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCa        >  1:3994208/1‑76 (MQ=255)                                                                cggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATAGGCaa       >  1:1313397/1‑76 (MQ=255)                                                                cggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCACCGGCaa       >  2:3412423/1‑76 (MQ=255)                                                                 ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa      <  1:2835090/76‑1 (MQ=255)                                                                 ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa      >  1:4496734/1‑76 (MQ=255)                                                                  gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc     >  2:1637919/1‑76 (MQ=255)                                                                   aTAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCCCTTCAACACGCGAAGGTCCATCGGCACAcc    >  1:4762489/1‑76 (MQ=255)                                                                   aTAAAACCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc    <  2:4287181/76‑1 (MQ=255)                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac  >  2:4113008/1‑76 (MQ=255)                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGCCCATCGGCAAACcac  >  2:1161618/1‑76 (MQ=255)                                                                             |                                                                     TGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCAC  >  NC_000913/3554013‑3554155

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