Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,478,540	A→C	E23A (GAA→GCA)	rraB →	ribonuclease E inhibitor protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,478,540	0	A	C	100.0%	143.4 / NA	42	E23A (GAA→GCA)	rraB	ribonuclease E inhibitor protein B
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (24/18);  total (24/18)

AATGAGCCATGGCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGAAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAACCCTGGAAAAAG  >  NC_000913/4478465‑4478617                                                                            |                                                                              aaTGAGCCATGGCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATg                                                                          >  2:857056/1‑81 (MQ=255)  aTGAGCCATGGCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATgg                                                                         >  1:857055/1‑81 (MQ=255)  aTGAGCCATGGCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATgg                                                                         <  2:857011/81‑1 (MQ=255)   tGAGCCATGGCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGc                                                                        <  2:857012/81‑1 (MQ=255)     aGCCATGGCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCa                                                                       <  2:857013/80‑1 (MQ=255)      gCCATGGCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGc                                                                     >  2:857043/1‑81 (MQ=255)      gCCATGGCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGc                                                                     <  1:857015/81‑1 (MQ=255)         aTGGCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGAt                                                                  >  2:857045/1‑81 (MQ=255)           ggCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGAAcc                                                                >  2:857053/1‑81 (MQ=255)              aaaCCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGa                                                             >  2:857052/1‑81 (MQ=255)                  ccGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGAcgcg                                                         >  2:857060/1‑81 (MQ=255)                     gAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATCACTGGCAGATGGCAGCGATCCGGACGCGCTg                                                      >  2:857066/1‑81 (MQ=255)                      aacaacTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTg                                                      >  1:857052/1‑80 (MQ=255)                      aacaacTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGt                                                     <  1:857016/81‑1 (MQ=255)                       acaacTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTa                                                    >  2:857040/1‑81 (MQ=255)                        caacTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTa                                                    >  2:857048/1‑80 (MQ=255)                              gaagaaCAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATc                                             <  2:857017/81‑1 (MQ=255)                              gaagaaCAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATc                                             <  1:857017/81‑1 (MQ=255)                                     aGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAacac                                       >  1:857073/1‑80 (MQ=255)                                      gCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCa                                     <  1:857019/81‑1 (MQ=255)                                        gTGAAGAAACACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATc                                   <  2:857019/81‑1 (MQ=255)                                              aaaCACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCg                             >  2:857073/1‑81 (MQ=255)                                               aaCACGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCg                             >  1:857069/1‑80 (MQ=255)                                                acacGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCg                             >  1:857098/1‑79 (MQ=255)                                                  acGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGa                         >  1:857081/1‑81 (MQ=255)                                                   cGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGAc                        <  2:857020/81‑1 (MQ=255)                                                   cGTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGAc                        <  2:857021/81‑1 (MQ=255)                                                    gTTTGATTATTGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACg                       <  2:857022/81‑1 (MQ=255)                                                           tattGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTgg                 >  2:857050/1‑80 (MQ=255)                                                           tattGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTgg                 >  1:857065/1‑80 (MQ=255)                                                            attGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGa                <  2:857023/80‑1 (MQ=255)                                                            attGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGa                <  1:857026/80‑1 (MQ=255)                                                            attGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGa                <  2:857024/80‑1 (MQ=255)                                                              tGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAAc             >  1:857068/1‑81 (MQ=255)                                                              tGAAGAATTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAAc             >  2:857084/1‑81 (MQ=255)                                                                aagaaTTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAAccc           >  1:857062/1‑81 (MQ=255)                                                                aagaaTTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAAcc            <  1:857027/80‑1 (MQ=255)                                                                 agaaTTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAACCCt          <  1:857028/81‑1 (MQ=255)                                                                  gaaTTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAACCCTg         <  2:857028/81‑1 (MQ=255)                                                                   aaTTACTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAACCCTg         >  1:857075/1‑80 (MQ=255)                                                                       aCTGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAACCCTGGaaa     >  2:857065/1‑80 (MQ=255)                                                                         tGGCAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAACCCTGGAAAAag  >  1:857070/1‑81 (MQ=255)                                                                            |                                                                              AATGAGCCATGGCAAACCCGGAACAACTGGAAGAACAGCGTGAAGAAACACGTTTGATTATTGAAGAATTACTGGAAGATGGCAGCGATCCGGACGCGCTGTACACCATCGAACACCATCTTTCCGCAGACGATCTGGAAACCCTGGAAAAAG  >  NC_000913/4478465‑4478617

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