Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,478,879	T→C	100%	V136A (GTT→GCT)	rraB →	protein inhibitor of RNase E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,478,879	0	T	C	100.0%	273.6 / ‑6.7	80	V136A (GTT→GCT)	rraB	protein inhibitor of RNase E
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (31/48);  total (31/49)

GACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCAT  >  NC_000913/4478785‑4478974                                                                                               |                                                                                                gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGcc                                                                                            <  1:1993315/100‑1 (MQ=255)  aCGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCa                                                                                           >  2:965019/1‑100 (MQ=255)   cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                          <  2:332572/100‑1 (MQ=255)   cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                          >  2:635228/1‑100 (MQ=255)   cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                          <  1:293042/100‑1 (MQ=255)   cGTTCAGTACGACGGATGGGGCACTTACTTTGAAGCTCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                          <  2:115881/100‑1 (MQ=255)     ttGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTa                                                                                        >  1:1097168/1‑100 (MQ=255)      tGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTaa                                                                                       >  2:613159/1‑100 (MQ=255)        aGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAAtt                                                                                     >  2:1028148/1‑100 (MQ=255)        aGTACGACGGATGGGGCACTTACTTTAAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTaa                                                                                       <  1:389741/98‑1 (MQ=255)         gTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTa                                                                                    >  2:1278348/1‑100 (MQ=255)          tACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATAACGGAGCTCGCCACTAATTaa                                                                                   >  2:2217733/1‑100 (MQ=255)           acgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAg                                                                                  >  1:554864/1‑100 (MQ=255)            cgacgGATGGGGCACTTACTTGGAGGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGAGGAAGACGATGACGGAGTTCGCCACTAATTAAg                                                                                  <  2:245287/99‑1 (MQ=255)             gacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCg                                                                                <  2:1304911/100‑1 (MQ=255)                 gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgc                                                                            <  1:2421146/100‑1 (MQ=255)                 gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgc                                                                            <  1:737076/100‑1 (MQ=255)                     gggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCAt                                                                        <  1:2033991/100‑1 (MQ=255)                      ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATc                                                                       >  1:1423550/1‑100 (MQ=255)                      ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATc                                                                       <  2:734284/100‑1 (MQ=255)                         acttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAgg                                                                    <  1:1290938/100‑1 (MQ=255)                          cttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGc                                                                   <  1:763937/100‑1 (MQ=255)                           ttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGcc                                                                  >  1:1642611/1‑100 (MQ=255)                           ttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGc                                                                   <  2:1157876/99‑1 (MQ=255)                              cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGTCCAg                                                               >  2:580444/1‑100 (MQ=255)                              cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                                                               <  1:1699024/100‑1 (MQ=255)                               tttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAgg                                                              >  1:2378457/1‑100 (MQ=255)                               tttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAgg                                                              <  2:1205051/100‑1 (MQ=255)                               tttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGAGGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                                                               <  2:1178961/99‑1 (MQ=255)                                ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                             <  1:1711329/100‑1 (MQ=255)                                ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                             <  1:635235/100‑1 (MQ=255)                                ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                             <  2:1529423/100‑1 (MQ=255)                                ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                             <  2:2132416/100‑1 (MQ=255)                                       tCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCt                                                      <  1:1954932/100‑1 (MQ=255)                                         cTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGc                                                    >  2:930832/1‑100 (MQ=255)                                         cTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGc                                                    <  2:1678107/100‑1 (MQ=255)                                         cTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGc                                                    >  1:221058/1‑100 (MQ=255)                                          tAACGGGGAAGAGGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGc                                                    <  2:2143247/99‑1 (MQ=255)                                          tAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCa                                                   >  1:330604/1‑100 (MQ=255)                                           aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCaa                                                  >  1:780997/1‑100 (MQ=255)                                             cGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAcc                                                >  2:1855330/1‑100 (MQ=255)                                             cGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAcc                                                >  1:1243217/1‑100 (MQ=255)                                               gCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGAGGGAGCTGGCCCCTAATAAAGCGGAGCGCAACAGACGGAGGTAGTACCTACAAcc                                                >  2:2263627/1‑98 (MQ=255)                                               gCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTg                                              >  2:2523561/1‑100 (MQ=255)                                                 gAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGt                                            <  2:1575138/100‑1 (MQ=255)                                                  aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTg                                           <  2:1955604/100‑1 (MQ=255)                                                    gATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGaa                                         <  1:1448041/100‑1 (MQ=255)                                                      tGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCATGCCCAGGTAGTTCCTGCAACCTGGTGAAtt                                       >  2:508554/1‑100 (MQ=255)                                                         cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGt                                    <  1:28683/100‑1 (MQ=255)                                                         cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGt                                    <  2:1455563/100‑1 (MQ=255)                                                               tGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCg                              <  1:375438/100‑1 (MQ=255)                                                                gAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCgg                             <  1:2002719/100‑1 (MQ=255)                                                                gAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCgg                             >  1:2154584/1‑100 (MQ=255)                                                                gAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCgg                             >  1:2154615/1‑100 (MQ=255)                                                                  aGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCCGGTAGTTCCTGCAACCTGGTGAATTTGTAGNCCGGAt                           >  1:2168603/1‑100 (MQ=255)                                                                  aGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGAt                           <  2:984654/100‑1 (MQ=255)                                                                   gATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCATGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATa                          <  2:1371409/100‑1 (MQ=255)                                                                   gATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATa                          <  2:642726/100‑1 (MQ=255)                                                                    aTTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATa                          <  2:965990/99‑1 (MQ=255)                                                                      tttGTCGATGAAGACGATTACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAgg                       <  1:1793235/100‑1 (MQ=255)                                                                         gTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCAt                    <  1:1369988/100‑1 (MQ=255)                                                                         gTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCAt                    >  1:2535948/1‑100 (MQ=255)                                                                         gTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATAAGGCCCAGGTAGTTCCTGCAACCTGGAGAATTTGTAGCCCGGGTAAGGCAt                    >  2:1906557/1‑100 (MQ=255)                                                                           cGATCAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTc                  <  2:2281413/100‑1 (MQ=255)                                                                            gATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                 >  1:608794/1‑100 (MQ=255)                                                                              tGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATg               <  1:528545/100‑1 (MQ=255)                                                                              tGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATg               <  2:1489384/100‑1 (MQ=255)                                                                              tGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATg               <  1:2490973/100‑1 (MQ=255)                                                                                aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc             <  2:1381888/100‑1 (MQ=255)                                                                                aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc             >  1:1927710/1‑100 (MQ=255)                                                                                aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc             >  1:2312898/1‑100 (MQ=255)                                                                                aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc             >  1:1748676/1‑100 (MQ=255)                                                                                aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGAGAAGGCATTCATGTc             <  2:338417/100‑1 (MQ=255)                                                                                 aGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTca            >  1:576236/1‑100 (MQ=255)                                                                                 aGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTca            <  1:192876/100‑1 (MQ=255)                                                                                     gATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATc        <  1:1846633/100‑1 (MQ=255)                                                                                     gATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATc        <  1:1847051/100‑1 (MQ=255)                                                                                       tGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCg      <  1:315889/100‑1 (MQ=255)                                                                                         aCGGAGCTCGCCACTAATTAAGCGGAGCGCATCAGGCCCAGGTAGGTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGc    <  1:815836/100‑1 (MQ=255)                                                                                           ggAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCAt  <  2:492983/100‑1 (MQ=255)                                                                                               |                                                                                                GACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCAT  >  NC_000913/4478785‑4478974

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