Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 612,648 G→T intergenic (‑154/‑89) fepA ← / → fes ferric enterobactin outer membrane transporter/enterochelin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913612,6480GT92.6% 317.5 / 19.5 95intergenic (‑154/‑89)fepA/fesferric enterobactin outer membrane transporter/enterochelin esterase
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base T (25/63);  total (28/67)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTCT  >  NC_000913/612601‑612695
                                               |                                               
gtgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTActct                                               <  1:4791676/50‑1 (MQ=255)
gtgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTActct                                               <  1:7365902/50‑1 (MQ=255)
gtgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACgat                                               <  1:3719240/50‑3 (MQ=255)
 tgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTa                                              <  1:6614518/50‑1 (MQ=255)
 tgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTa                                              <  1:2366257/50‑1 (MQ=255)
 tgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTa                                              <  1:6981175/50‑1 (MQ=255)
  gcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAt                                             >  1:6276871/1‑50 (MQ=255)
   cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAtt                                            >  1:6418461/1‑50 (MQ=255)
   cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAtt                                            <  1:612862/50‑1 (MQ=255)
   cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAtt                                            <  1:1994034/50‑1 (MQ=255)
   cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAtt                                            <  1:1721919/50‑1 (MQ=255)
      ttGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgca                                         <  1:4677055/50‑1 (MQ=255)
      ttGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgca                                         <  1:3814326/50‑1 (MQ=255)
        ggTTAATTAATGTCCGGGCTTCCCACGGCGCGCCATTACTCTATtgcaa                                        <  1:333100/49‑1 (MQ=255)
        ggTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaa                                       >  1:7564575/1‑50 (MQ=255)
        ggTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaa                                        <  1:7572134/49‑1 (MQ=255)
        ggTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaa                                        <  1:625638/49‑1 (MQ=255)
        ggTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaa                                        <  1:1344830/49‑1 (MQ=255)
         gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaat                                      <  1:2491567/50‑1 (MQ=255)
         gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaat                                      <  1:4592785/50‑1 (MQ=255)
         gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaa                                       <  1:827031/49‑1 (MQ=255)
         gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaa                                       <  1:7334293/49‑1 (MQ=255)
          ttaattaaTGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatg                                     >  1:5959615/1‑50 (MQ=255)
            aattaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgc                                    >  1:3444349/1‑49 (MQ=255)
             attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaa                                  <  1:7203905/50‑1 (MQ=255)
             attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaa                                  <  1:5461775/50‑1 (MQ=255)
             attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgca                                   <  1:4192406/49‑1 (MQ=255)
             attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgca                                   <  1:1137152/49‑1 (MQ=255)
             attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgca                                   >  1:952653/1‑49 (MQ=255)
              ttaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaaa                                 <  1:138946/50‑1 (MQ=255)
               taatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaaa                                 >  1:4979408/1‑49 (MQ=255)
               taatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaaa                                 >  1:1795565/1‑49 (MQ=255)
                aatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaaa                                 >  1:6357474/1‑48 (MQ=255)
                aatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTgcaaatgcaaat                                <  1:1469767/49‑1 (MQ=255)
                 atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATa                               <  1:1428899/49‑1 (MQ=255)
                  tGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGt                             <  1:847182/50‑1 (MQ=255)
                    tCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGtt                            <  1:6726440/49‑1 (MQ=255)
                    tCCGCGCTTCCCACGGCGCGCCATTACGATATTGCAAATGCAAATAGtt                            <  1:2984836/49‑1 (MQ=255)
                      cgcgCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTAt                          <  1:623862/49‑1 (MQ=255)
                      cgcgCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTAt                          <  1:3242020/49‑1 (MQ=255)
                      cgcgCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTAt                          <  1:996873/49‑1 (MQ=255)
                        cgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaa                       <  1:1421841/50‑1 (MQ=255)
                        cgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaa                       <  1:7186708/50‑1 (MQ=255)
                        cgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaa                       >  1:866068/1‑50 (MQ=255)
                         gcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaat                      <  1:2756102/50‑1 (MQ=255)
                         gcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGCTATCaa                       >  1:643549/1‑49 (MQ=39)
                          cTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaata                     <  1:1117650/50‑1 (MQ=255)
                          cTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaat                      >  1:4179058/1‑49 (MQ=255)
                          cTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaat                      >  1:67029/1‑49 (MQ=255)
                          cTTCCCACGGCGCGCCATTACGATATTGCAAATGCAAATAGTTATCaat                      >  1:4069347/1‑49 (MQ=255)
                           ttCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaataa                    >  1:3729188/1‑50 (MQ=255)
                              ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAataata                  <  1:3141831/49‑1 (MQ=255)
                              ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAataata                  <  1:5950237/49‑1 (MQ=255)
                              ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAataata                  <  1:5942136/49‑1 (MQ=255)
                              ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAatat                 <  1:6913194/50‑1 (MQ=255)
                               cACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAatat                 >  1:7222955/1‑49 (MQ=255)
                               cACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatt                <  1:1859218/50‑1 (MQ=255)
                               cACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatt                <  1:6451259/50‑1 (MQ=255)
                               cACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatt                <  1:5085471/50‑1 (MQ=255)
                                aCGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatta               <  1:2125801/50‑1 (MQ=255)
                                aCGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatta               <  1:1655972/50‑1 (MQ=255)
                                aCGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatta               <  1:1224606/50‑1 (MQ=255)
                                aCGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatta               <  1:278407/50‑1 (MQ=255)
                                 cGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtattat              <  1:6114760/50‑1 (MQ=255)
                                  ggCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATc             <  1:3746693/50‑1 (MQ=255)
                                  ggCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATc             <  1:7374720/50‑1 (MQ=255)
                                  ggCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATc             <  1:1050163/50‑1 (MQ=255)
                                   gcgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCa            <  1:2162763/50‑1 (MQ=255)
                                   gcgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCa            <  1:3617045/50‑1 (MQ=255)
                                   gcgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCa            <  1:1549950/50‑1 (MQ=255)
                                   gcgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCa            <  1:4365878/50‑1 (MQ=255)
                                   gcgcgcCATTACGATATTGCAAATGCAAATAGTTATCAATAATATTATCa            <  1:4135758/50‑1 (MQ=255)
                                    cgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCaa           >  1:490839/1‑50 (MQ=255)
                                     gcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAGTATTATCAat          <  1:7397775/50‑1 (MQ=39)
                                     gcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCaa           >  1:3291686/1‑49 (MQ=255)
                                     gcgcCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCaa           >  1:2785465/1‑49 (MQ=255)
                                      cgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAata         <  1:6448714/50‑1 (MQ=255)
                                      cgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAata         <  1:3678326/50‑1 (MQ=255)
                                       gcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatat        <  1:2960498/50‑1 (MQ=255)
                                       gcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatat        <  1:2533036/50‑1 (MQ=255)
                                       gcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatat        <  1:6035965/50‑1 (MQ=255)
                                       gcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAata         >  1:227229/1‑49 (MQ=255)
                                       gcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAata         >  1:3987046/1‑49 (MQ=255)
                                       gcCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatat        <  1:6775534/50‑1 (MQ=255)
                                        ccATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatata       <  1:1375336/50‑1 (MQ=255)
                                        ccATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatata       >  1:4416425/1‑50 (MQ=255)
                                         cATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatatat      <  1:157872/50‑1 (MQ=255)
                                         cATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatatat      <  1:146067/50‑1 (MQ=255)
                                          aTTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAGTATAtt     <  1:5716827/50‑1 (MQ=38)
                                          aTTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAtt     >  1:895473/1‑50 (MQ=255)
                                          aTTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAtt     <  1:2713015/50‑1 (MQ=255)
                                          aTTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAtt     <  1:3199041/50‑1 (MQ=255)
                                           ttACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAttt    >  1:4261458/1‑50 (MQ=255)
                                           ttACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAttt    <  1:6079596/50‑1 (MQ=255)
                                           ttACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAttt    >  1:6209807/1‑50 (MQ=255)
                                            tACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTc   >  1:7301934/1‑50 (MQ=255)
                                            tACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTc   >  1:2611416/1‑50 (MQ=255)
                                             aCGATATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTCt  >  1:4054046/1‑50 (MQ=255)
                                               |                                               
GTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTCT  >  NC_000913/612601‑612695

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