Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,831 C→T 16.9% intergenic (+23/+17) tonB → / ← yciA membrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,8310CT16.9% 169.9 / 24.4 78intergenic (+23/+17)tonB/yciAmembrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base C (33/29);  new base T (8/5);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

AGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATA  >  NC_000913/1311735‑1311925
                                                                                                |                                                                                              
aGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCgg                                                                                             <  2:1892093/100‑1 (MQ=255)
 gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                            <  2:837258/100‑1 (MQ=255)
  ccGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAg                                                                                           >  1:938688/1‑100 (MQ=255)
  ccGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAg                                                                                           <  2:138611/100‑1 (MQ=255)
    gggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGc                                                                                         <  2:2103657/100‑1 (MQ=255)
       tAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCttt                                                                                      <  2:68914/100‑1 (MQ=255)
        aaGCCAGGCAGTGGGGTTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                     <  2:2017928/100‑1 (MQ=255)
        aaGCCAGGCAGTGGGATTGTGGCGAATATCCTGTTTAAAATTAACGGCACCACCGAAATGCAGTAAGCAGAAAGTCAAAAGCCCCCCACCGGAGGCt                                                                                        >  2:420885/1‑97 (MQ=255)
         aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGTCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                    <  1:2320444/100‑1 (MQ=39)
         aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                    <  1:1130993/100‑1 (MQ=255)
         aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                    <  2:1106761/100‑1 (MQ=255)
         aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                    <  1:2062103/100‑1 (MQ=255)
         aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                    >  1:2227502/1‑100 (MQ=255)
          gCCAGGCAGTGGGAGTGTGGTGAATATCCTGTTTAAAATTAACGGCACCAACGAAATTCAGTAAGTAGAAAGTCAAAAGCCTCCGAACGGAGTCtat                                                                                      >  2:1514759/1‑95 (MQ=255)
           ccAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGAc                                                                                  <  2:2131972/100‑1 (MQ=255)
           ccAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                  <  1:2070915/100‑1 (MQ=255)
           ccAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                  <  1:2071843/100‑1 (MQ=255)
            cAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                  <  1:799412/99‑1 (MQ=255)
                     ggATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTATTGACTATTACTCaa                                                                        >  1:2367392/1‑100 (MQ=255)
                         ggagatgAACATCCTCTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACAGGATGATCTTGACCATTCCTAAACa                                                                      >  2:1247567/6‑98 (MQ=255)
                               ggatatCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTaaaac                                                              >  2:921261/3‑97 (MQ=255)
                                 atatCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGgc                                                              >  1:711391/1‑98 (MQ=255)
                                    tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                         <  1:2392478/100‑1 (MQ=255)
                                     ccTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAgg                                                        <  1:1339457/100‑1 (MQ=255)
                                        gTTTATAATTGACGGGACCACCGAAATTCAGTAAGCAGTAAGTCGTAAGCCGCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                     <  2:341806/100‑1 (MQ=255)
                                        gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                     >  1:1531874/1‑100 (MQ=255)
                                        gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                     >  1:1005699/1‑100 (MQ=255)
                                         tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                                    <  1:851034/100‑1 (MQ=255)
                                           tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                  >  1:1760435/1‑100 (MQ=255)
                                           tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                  >  2:611462/1‑100 (MQ=255)
                                           tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAATCAAAACCCTCCGAACGGAGGCTTTTGACTATTACTCAACACGTAAGGCGCGAGATTTTcc                                                  >  2:962961/1‑100 (MQ=255)
                                            aaaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCt                                                 <  1:950673/100‑1 (MQ=255)
                                             aaaTTAACGTTTCCATCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCt                                                 <  1:803612/99‑1 (MQ=39)
                                               aTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGATAAGGCACAAGGTTTTCCTTCa                                              >  1:29249/1‑100 (MQ=255)
                                                 tAACGGCACCACCGAAATTCAGTAAGCAGAAATTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACACGTAAACCGCGAGGTTTTTCttta                                              >  2:1585784/1‑96 (MQ=255)
                                                 tAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAgg                                            <  1:724302/100‑1 (MQ=255)
                                                    cGGCACCACAGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATc                                         >  1:455200/1‑100 (MQ=255)
                                                     ggCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAACGGAGGCTTTTGACTATTACTCACCAGGTAAGGCACGAGGTTTTCCTTCAGGATCa                                        >  2:2223350/1‑100 (MQ=255)
                                                      gCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCaa                                       >  2:919963/1‑100 (MQ=255)
                                                        accaccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAAcc                                     >  2:469100/1‑100 (MQ=255)
                                                        accaccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAAGTAAGGCGCGAGGTTTTCCTTCAGGATCaa                                       >  2:583921/1‑98 (MQ=255)
                                                          caccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgc                                   >  1:269525/1‑100 (MQ=255)
                                                            ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGATGTTTTCCTTCAGGATCAACCGCGa                                 >  1:843765/1‑100 (MQ=255)
                                                            ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAAGCTTTTGACTATTACTCAAGAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGa                                 >  2:1304066/1‑100 (MQ=255)
                                                              gAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGAAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGATTTCGTTCAGAATCAACCGCGAc                                >  2:804729/1‑99 (MQ=255)
                                                                  ttCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACt                           <  1:90513/100‑1 (MQ=255)
                                                                     aGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaa                        >  2:1043/1‑100 (MQ=255)
                                                                     aGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaa                        >  2:1889165/1‑100 (MQ=255)
                                                                      gTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                       <  1:1378232/100‑1 (MQ=255)
                                                                      gTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                       <  1:1583756/100‑1 (MQ=255)
                                                                       tAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat                      <  1:420888/100‑1 (MQ=255)
                                                                          gCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataat                   <  1:611465/100‑1 (MQ=255)
                                                                           cAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATg                  >  1:2285927/1‑100 (MQ=255)
                                                                           cAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATg                  <  2:511037/100‑1 (MQ=255)
                                                                             gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCt                >  1:955748/1‑100 (MQ=255)
                                                                             gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCt                >  1:1355286/1‑100 (MQ=255)
                                                                              aaaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCtt               <  2:1372462/100‑1 (MQ=255)
                                                                               aaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc              <  2:765802/100‑1 (MQ=255)
                                                                                aGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTTAGGATCAACAGCGACATACTTAAATAATGAttat             >  2:565476/1‑98 (MQ=37)
                                                                                 gTCAAAAGCCTCCGGTCGGAGGATTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTg            >  2:416009/1‑100 (MQ=255)
                                                                                 gTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCACCCGCGACATACTTAAATAATGCTTCTg            >  2:1927658/1‑100 (MQ=255)
                                                                                  tCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGt           >  1:1344064/1‑100 (MQ=255)
                                                                                  tCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTTCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGt           <  2:1323052/100‑1 (MQ=255)
                                                                                   cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc          >  1:1064100/1‑100 (MQ=255)
                                                                                   cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc          >  2:1190578/1‑100 (MQ=255)
                                                                                   cAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc          >  2:853714/1‑100 (MQ=255)
                                                                                   cAAAAGCCTCCGACCGGAGGCTTTTGACTATTAATCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACAGCGACATACTTAAATAATGCTTCTGTc          >  2:82482/1‑100 (MQ=255)
                                                                                    naaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc          >  2:2175553/2‑99 (MQ=255)
                                                                                    aaaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc          >  1:1095359/1‑99 (MQ=255)
                                                                                     aaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGc        >  1:2402115/1‑100 (MQ=255)
                                                                                     aaaGCCTCCGACCGAAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGc        >  1:2176961/1‑100 (MQ=255)
                                                                                        gCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCttt     >  1:1324667/1‑100 (MQ=255)
                                                                                        gCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGATAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCttt     >  1:1690963/1‑100 (MQ=255)
                                                                                        gCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCttt     <  1:1514776/100‑1 (MQ=255)
                                                                                         ccTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTta    <  2:2385582/100‑1 (MQ=255)
                                                                                         ccTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTta    <  1:2130771/100‑1 (MQ=255)
                                                                                           tCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTtata  <  1:921272/100‑1 (MQ=255)
                                                                                              gACCGGTGGCTTATGACTATAACACAACAGGAAAGACACAAGGCTTTCCTTCAGGACCAACCGCGACATACTTAAATAATGTTTCCGTc          >  2:359478/1‑89 (MQ=255)
                                                                                                |                                                                                              
AGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATA  >  NC_000913/1311735‑1311925

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