Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,478,943	C→T	intergenic (+54/+91)	rraB → / ← yjgM	ribonuclease E inhibitor protein B/putative acetyltransferase YjgM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,478,943	0	C	T	100.0%	222.1 / NA	62	intergenic (+54/+91)	rraB/yjgM	ribonuclease E inhibitor protein B/putative acetyltransferase YjgM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/33);  total (29/33)

ACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGG  >  NC_000913/4478809‑4479073                                                                                                                                       |                                                                                                                                   acttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGa                                                                                                                               <  1:261312/140‑1 (MQ=255)   ttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATa                                                                                                                             <  2:488479/140‑1 (MQ=255)        ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCa                                                                                                                        >  1:466662/1‑140 (MQ=255)         tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCAt                                                                                                                       >  2:946522/1‑140 (MQ=255)          gAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCAtt                                                                                                                      <  2:195826/140‑1 (MQ=255)           aaGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTc                                                                                                                     <  1:674830/140‑1 (MQ=255)           aaGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTc                                                                                                                     >  1:557677/1‑140 (MQ=255)               tCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCa                                                                                                                        <  1:1111921/133‑1 (MQ=255)               tCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCa                                                                                                                        >  2:1111921/1‑133 (MQ=255)               tCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGt                                                                                                                 <  2:514546/140‑1 (MQ=255)                ccTAACGGCGAAGATGGCGAC‑‑‑GAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTcaca                                                                                                             >  1:1066908/1‑140 (MQ=255)                   aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTcaca                                                                                                             >  2:925869/1‑140 (MQ=255)                   aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTcaca                                                                                                             >  2:415204/1‑140 (MQ=255)                     cGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATc                                                                                                           <  1:521828/140‑1 (MQ=255)                     cGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATc                                                                                                           >  2:1045715/1‑140 (MQ=255)                          aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACAt                                                                                                            <  1:373945/134‑1 (MQ=255)                          aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACAt                                                                                                            >  2:373945/1‑134 (MQ=255)                          aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCa                                                                                                      <  1:382426/140‑1 (MQ=255)                          aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCa                                                                                                      <  2:445976/140‑1 (MQ=255)                          aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCa                                                                                                      >  1:935440/1‑140 (MQ=255)                               ggCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAAc                                                                                                 >  2:561493/1‑140 (MQ=255)                                gCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACa                                                                                                <  2:1122545/140‑1 (MQ=255)                                      aTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCa                                                                                                      <  2:926827/128‑1 (MQ=255)                                      aTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCa                                                                                                      >  1:926827/1‑128 (MQ=255)                                       tGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAacgta                                                                                         <  1:56404/140‑1 (MQ=255)                                       tGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAacgta                                                                                         <  2:935440/140‑1 (MQ=255)                                       tGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAacgta                                                                                         <  1:380785/140‑1 (MQ=255)                                       tGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAacgta                                                                                         <  1:747381/140‑1 (MQ=255)                                       tGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAacgta                                                                                         >  2:228079/1‑140 (MQ=255)                                       tGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAacgta                                                                                         <  1:925869/140‑1 (MQ=255)                                             ttttGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACgttgt                                                                                   <  1:415204/140‑1 (MQ=255)                                               ttGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTa                                                                                 <  1:1045715/140‑1 (MQ=255)                                                          gACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATg                                                                      <  1:228079/140‑1 (MQ=255)                                                             gATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGccc                                                                   <  1:561493/140‑1 (MQ=255)                                                             gATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGccc                                                                   >  2:356971/1‑140 (MQ=255)                                                              aTGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcg                                                                  <  1:356971/140‑1 (MQ=255)                                                                gACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCcgcg                                                                <  1:385616/140‑1 (MQ=255)                                                                  cGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAg                                                              <  2:557677/140‑1 (MQ=255)                                                                         cGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                       <  2:287163/140‑1 (MQ=255)                                                                                           gcgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTaga                                     >  1:588962/1‑140 (MQ=255)                                                                                           gcgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTaga                                     >  2:1147004/1‑140 (MQ=255)                                                                                            cgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTagag                                    <  1:235906/140‑1 (MQ=255)                                                                                            cgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTagag                                    >  2:771689/1‑140 (MQ=255)                                                                                            cgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTagag                                    >  2:318334/1‑140 (MQ=255)                                                                                            cgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTagag                                    >  2:207618/1‑140 (MQ=255)                                                                                             gcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGt                                   <  1:911256/140‑1 (MQ=255)                                                                                                tCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcg                                >  1:1099441/1‑140 (MQ=255)                                                                                                 cAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgc                               >  1:736471/1‑140 (MQ=255)                                                                                                         ggTAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATc                       >  1:609160/1‑140 (MQ=255)                                                                                                           tAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATccc                     <  1:140850/140‑2 (MQ=255)                                                                                                           tAGTTCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCg                     <  1:207618/140‑1 (MQ=255)                                                                                                              ttCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                  >  2:609788/1‑140 (MQ=255)                                                                                                              ttCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                  >  2:180849/1‑140 (MQ=255)                                                                                                              ttCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                  >  2:111296/1‑140 (MQ=255)                                                                                                               tCCTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACt                 <  1:318334/140‑1 (MQ=255)                                                                                                                ccTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCgcgc                             <  1:91393/127‑1 (MQ=255)                                                                                                                ccTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCgcgc                             >  2:91393/1‑127 (MQ=255)                                                                                                                ccTGCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACtt                <  1:1147004/140‑1 (MQ=255)                                                                                                                   gCAACCTGGTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcac             >  2:199277/1‑140 (MQ=255)                                                                                                                          ggTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCg                     <  2:440562/125‑1 (MQ=255)                                                                                                                          ggTGAATTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCg                     >  1:440562/1‑125 (MQ=255)                                                                                                                              aaTTTGTAGTCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATgg  >  2:668440/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                   ACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGG  >  NC_000913/4478809‑4479073

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