Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 619,892 G→A A103V (GCG→GTG)  fepC ← ferrienterobactin ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913619,8920GA97.5% 256.8 / ‑4.2 79A103V (GCG→GTG) fepCferrienterobactin ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (34/43);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGTTATTCCCGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCGCCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTGGCGTAATGTTG  >  NC_000913/619796‑619986
                                                                                                |                                                                                              
tGAGTTATTCCCGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCAcca                                                                                             >  1:1454659/1‑100 (MQ=255)
 gAGTTATTCCCGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGcaccac                                                                                            >  1:193379/1‑100 (MQ=255)
   gTTATTCCCGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCAccacca                                                                                          <  2:1993458/100‑1 (MQ=255)
      aTTCCCGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCt                                                                                       >  1:2217201/1‑100 (MQ=255)
        tCCCGTGGCCTGCATTGCTTTCGTTCCCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTcc                                                                                     <  2:1772525/100‑1 (MQ=255)
         cccGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCt                                                                                    >  2:1146634/1‑100 (MQ=255)
         cccGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCt                                                                                    >  2:2463846/1‑100 (MQ=255)
         cccGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCt                                                                                    >  1:2234316/1‑100 (MQ=255)
         cccGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGCTTGATGCGGATAACGTCCACGCACCACCAGCTCCt                                                                                    >  2:1903000/1‑100 (MQ=255)
          ccGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTg                                                                                   <  2:1278353/100‑1 (MQ=255)
          ccGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTg                                                                                   >  1:15178/1‑100 (MQ=255)
          ccGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTg                                                                                   <  2:1773874/100‑1 (MQ=255)
          ccGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTg                                                                                   <  1:964862/100‑1 (MQ=255)
           cGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGc                                                                                  >  1:1895458/1‑100 (MQ=255)
            gTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCa                                                                                 <  2:1554757/100‑1 (MQ=255)
               gCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACgg                                                                              <  2:551624/100‑1 (MQ=255)
                ccTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGt                                                                             <  1:1723557/100‑1 (MQ=255)
                 cTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTg                                                                            <  1:646532/100‑1 (MQ=255)
                 cTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTg                                                                            >  1:1654263/1‑100 (MQ=255)
                        gCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGc                                                                     >  2:1090026/1‑100 (MQ=255)
                            tCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCgg                                                                 <  2:2396675/100‑1 (MQ=255)
                            tCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCgg                                                                 >  1:2231995/1‑100 (MQ=255)
                            tCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCgg                                                                 <  1:2374047/100‑1 (MQ=255)
                            tCGTTACCGCGTCTTCATCCTCTTGGCGCCGGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATAGCGCCCgg                                                                 <  2:950092/100‑1 (MQ=255)
                              gTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCg                                                               >  2:799280/1‑100 (MQ=255)
                                  ccGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGt                                                           <  2:1676321/100‑1 (MQ=255)
                                  ccGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGt                                                           <  1:1803395/100‑1 (MQ=255)
                                  ccGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGt                                                           >  2:1374261/1‑100 (MQ=255)
                                  ccGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGt                                                           >  2:1334716/1‑100 (MQ=255)
                                  ccGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGt                                                           >  1:1593496/1‑100 (MQ=255)
                                   cGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTa                                                          >  1:1621263/1‑100 (MQ=255)
                                    gCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAg                                                         <  1:2066403/100‑1 (MQ=255)
                                       tcttcATCATCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCAt                                                      >  1:614597/1‑100 (MQ=255)
                                         ttcATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCAttt                                                    >  1:2152983/1‑100 (MQ=255)
                                            aTCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTgc                                                 <  1:1303062/100‑1 (MQ=255)
                                            aTCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTgc                                                 <  2:2044316/100‑1 (MQ=255)
                                            aTCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACAGTGATATCGCCCGGCGTGGTAGCATTTTgc                                                 <  2:1805390/100‑1 (MQ=255)
                                             tCCTCTTTGCGCCAGTGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTgcg                                                <  1:141300/100‑1 (MQ=255)
                                             tCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTgcg                                                >  2:1369481/1‑100 (MQ=255)
                                              ccTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTgcgc                                               >  1:1469604/1‑100 (MQ=255)
                                                tctTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCa                                             <  2:1449245/100‑1 (MQ=255)
                                                tctTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCa                                             >  1:987780/1‑100 (MQ=255)
                                                  tttGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAAc                                           <  2:1813962/100‑1 (MQ=255)
                                                  tttGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAAc                                           <  1:1353994/100‑1 (MQ=255)
                                                   ttGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCGCCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACa                                          <  1:2140226/100‑1 (MQ=255)
                                                    tGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAg                                         <  2:2136047/100‑1 (MQ=255)
                                                     gcgcCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGa                                        >  1:989184/1‑100 (MQ=255)
                                                     gcgcCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGa                                        <  2:1434456/100‑1 (MQ=255)
                                                      cgcCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGAc                                       <  1:2635531/100‑1 (MQ=255)
                                                      cgcCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGAc                                       <  1:1084110/100‑1 (MQ=255)
                                                      cgcCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGAc                                       >  2:1531611/1‑100 (MQ=255)
                                                      cgcCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGAc                                       <  2:1612965/100‑1 (MQ=255)
                                                       gcCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGAcc                                      >  2:1778461/1‑100 (MQ=255)
                                                       gcCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGAcc                                      <  1:62895/100‑1 (MQ=255)
                                                                tAATCAGCGGTTGATGCGGATAACGTCCACGCATCACAAGCTACTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCAGCg                             <  2:858170/100‑1 (MQ=255)
                                                                  aaCAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTg                           <  2:872961/100‑1 (MQ=255)
                                                                    cAGCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCa                         >  1:2081893/1‑100 (MQ=255)
                                                                      gCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAGCGGACCAATCCGGCGTGCAAc                       <  1:2256844/100‑1 (MQ=255)
                                                                      gCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAAc                       <  2:2206674/100‑1 (MQ=255)
                                                                      gCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAAc                       >  1:64988/1‑100 (MQ=255)
                                                                      gCGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAAc                       <  1:2288648/100‑1 (MQ=255)
                                                                       cGGTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAAcc                      >  1:698322/1‑100 (MQ=255)
                                                                         gTTGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACctc                    <  2:732807/100‑1 (MQ=255)
                                                                          ttGATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACctct                   <  2:356617/100‑1 (MQ=255)
                                                                            gATGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCttt                 <  2:466117/100‑1 (MQ=255)
                                                                             ntGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTGTa                >  1:2113539/2‑100 (MQ=255)
                                                                              tGCGGATAACGTCCACGCGCCGCCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTAc               <  1:311282/100‑1 (MQ=255)
                                                                              tGCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGGAGCATCTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTAc               >  2:164691/1‑100 (MQ=255)
                                                                               gCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCCATCCGGCGTGCAACCTCTTTACt              >  2:1669422/1‑100 (MQ=255)
                                                                               gCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACt              <  2:1654904/100‑1 (MQ=255)
                                                                               gCGGATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACt              <  1:2187772/100‑1 (MQ=255)
                                                                                 ggATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTgg            <  2:1817429/100‑1 (MQ=255)
                                                                                 ggATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTgg            >  1:355913/1‑100 (MQ=255)
                                                                                  gATAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTGGc           <  1:2334771/100‑1 (MQ=255)
                                                                                   aTAACGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTGGCg          <  1:2379519/100‑1 (MQ=255)
                                                                                     aaCGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTGGCGTa        <  2:1216551/100‑1 (MQ=255)
                                                                                       cGTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTGGCGTAAt      >  1:1923359/1‑100 (MQ=255)
                                                                                        gTCCACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTGGCGTAATg     <  1:1927229/100‑1 (MQ=255)
                                                                                           cACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTGGCGTAATGTTg  >  1:1809168/1‑100 (MQ=255)
                                                                                           cACGCACCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTGGCGTAATGTTg  >  2:1727774/1‑100 (MQ=255)
                                                                                                |                                                                                              
TGAGTTATTCCCGTGGCCTGCATTGCTTTCGTTACCGCTTCTTCATCCTCTTTGCGCCAGCGGGTAAACAGCGGTTGATGCGGATAACGTCCACGCGCCACCAGCTCCTGCACGGTGATATCGCCCGGCGTGGTAGCATTTTGCGCCAACAGACCAATCCGGCGTGCAACCTCTTTACTGGCGTAATGTTG  >  NC_000913/619796‑619986

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