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,6480GT85.5% 371.5 / 56.4 125intergenic (‑154/‑89)fepA/fesferric enterobactin outer membrane transporter/enterochelin esterase
Reads supporting (aligned to +/- strand):  ref base G (3/15);  new base T (28/78);  total (31/94)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 129 total aligned reads displayed.

CAGTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTCTGC  >  NC_000913/612599‑612697
                                                 |                                                 
cagtgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACg                                                   <  1:5069512/50‑1 (MQ=255)
  gtgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTActct                                                 <  1:6433099/50‑1 (MQ=255)
  gtgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTActct                                                 <  1:5901876/50‑1 (MQ=255)
   tgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTa                                                <  1:1111827/50‑1 (MQ=255)
   tgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTa                                                >  1:3077660/1‑50 (MQ=255)
   tgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTa                                                <  1:2456881/50‑1 (MQ=255)
     cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAtt                                              <  1:851748/50‑1 (MQ=255)
     cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAtt                                              <  1:6568233/50‑1 (MQ=255)
     cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAtt                                              <  1:212754/50‑1 (MQ=255)
     cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAt                                               <  1:2907958/49‑1 (MQ=255)
      agtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtg                                             <  1:1089118/50‑1 (MQ=255)
      agtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtg                                             <  1:7808398/50‑1 (MQ=255)
      agtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtg                                             >  1:7784890/1‑50 (MQ=255)
      agtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtg                                             >  1:2152320/1‑50 (MQ=255)
      agtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATtg                                             <  1:4102387/50‑1 (MQ=255)
       gtTGGTTAATTAATGTCCGCGCTTCCCACAGCGCGCCATTACTCTATtgc                                            <  1:1965362/50‑1 (MQ=255)
        ttGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgca                                           <  1:6909217/50‑1 (MQ=255)
        ttGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgca                                           <  1:2840078/50‑1 (MQ=255)
          ggTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaa                                          <  1:4941729/49‑1 (MQ=255)
          ggTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATtgcaaa                                         >  1:5531983/1‑50 (MQ=255)
          ggTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATtgcaa                                          <  1:1934136/49‑1 (MQ=255)
           gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaat                                        >  1:3768638/1‑50 (MQ=255)
           gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaat                                        <  1:4826570/50‑1 (MQ=255)
           gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaat                                        >  1:4846691/1‑50 (MQ=255)
           gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaa                                         <  1:5067467/49‑1 (MQ=255)
           gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATtgcaaa                                         <  1:1624942/49‑1 (MQ=255)
            ttaattaaTGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatg                                       >  1:807695/1‑50 (MQ=255)
            ttaattaaTGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatg                                       <  1:4239382/50‑1 (MQ=255)
            ttaattaaTGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatg                                       <  1:4583837/50‑1 (MQ=255)
            ttaattaaTGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatg                                       >  1:2444590/1‑50 (MQ=255)
              aattaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgc                                      <  1:6552663/49‑1 (MQ=255)
               attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaa                                    <  1:2188974/50‑1 (MQ=255)
               attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgca                                     >  1:4536435/1‑49 (MQ=255)
                ttaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaaa                                   <  1:1428782/50‑1 (MQ=255)
                ttaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaaa                                   >  1:7018172/1‑50 (MQ=255)
                  aatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaaa                                   >  1:5617094/1‑48 (MQ=255)
                   atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATa                                 <  1:6393502/49‑1 (MQ=255)
                   atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATa                                 <  1:6874665/49‑1 (MQ=255)
                   atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATa                                 <  1:1083762/49‑1 (MQ=255)
                   atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAg                                >  1:6532783/1‑50 (MQ=255)
                   atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAg                                <  1:7771384/50‑1 (MQ=255)
                   atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAg                                <  1:1720080/50‑1 (MQ=255)
                   atGTCCGCGCTTCCCACGGCGCGCCATTACACTATTGCAAATGCAAATa                                 <  1:3634944/49‑1 (MQ=255)
                    tGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGt                               >  1:743545/1‑50 (MQ=255)
                    tGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGt                               <  1:4856257/50‑1 (MQ=255)
                      tCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGtt                              <  1:295227/49‑1 (MQ=255)
                      tCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGtt                              <  1:171377/49‑1 (MQ=255)
                      tCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGtt                              <  1:1192236/49‑1 (MQ=255)
                      tCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTa                             <  1:4812778/50‑1 (MQ=255)
                       ccGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTAt                            >  1:270699/1‑50 (MQ=255)
                        cgcgCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTAt                            <  1:1214975/49‑1 (MQ=255)
                        cgcgCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTAt                            >  1:8241446/1‑49 (MQ=255)
                         gcgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCa                          <  1:6930072/50‑1 (MQ=255)
                         gcgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCa                          <  1:1807076/50‑1 (MQ=255)
                         gcgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCa                          <  1:6401591/50‑1 (MQ=255)
                          cgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaa                         <  1:6200060/50‑1 (MQ=255)
                           gcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaat                        <  1:3606551/50‑1 (MQ=255)
                           gcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaat                        <  1:6079647/50‑1 (MQ=255)
                           gcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaa                         >  1:7057940/1‑49 (MQ=255)
                           gcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaa                         >  1:2483634/1‑49 (MQ=255)
                            cTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaata                       <  1:3113880/50‑1 (MQ=255)
                             ttCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaataa                      >  1:3087262/1‑50 (MQ=255)
                               cccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaataat                     <  1:4785421/49‑1 (MQ=255)
                               cccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaataat                     <  1:7706243/49‑1 (MQ=255)
                               cccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaataat                     <  1:3784275/49‑1 (MQ=255)
                               cccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaataat                     <  1:6211489/49‑1 (MQ=255)
                               cccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAataata                    >  1:4545572/1‑50 (MQ=255)
                               cccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAataata                    >  1:6941621/1‑50 (MQ=255)
                               cccACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTATCaataa                      <  1:674005/48‑1 (MQ=255)
                                ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAataata                    <  1:6278687/49‑1 (MQ=255)
                                ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAatat                   <  1:3645107/50‑1 (MQ=255)
                                ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAatat                   >  1:1780203/1‑50 (MQ=255)
                                ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAatat                   <  1:2417527/50‑1 (MQ=255)
                                 cACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatt                  <  1:6394969/50‑1 (MQ=255)
                                 cACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatt                  <  1:6671376/50‑1 (MQ=255)
                                  aCGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatta                 <  1:2698568/50‑1 (MQ=255)
                                   cGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtattat                <  1:7001602/50‑1 (MQ=255)
                                   cGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtattat                <  1:5445190/50‑1 (MQ=255)
                                   cGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtattat                <  1:453539/50‑1 (MQ=255)
                                   cGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtattat                <  1:7546949/50‑1 (MQ=255)
                                   cGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtattat                <  1:8109344/50‑1 (MQ=255)
                                   cGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTATCAATAAtattat                <  1:3414978/50‑1 (MQ=255)
                                    ggCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATc               <  1:3245112/50‑1 (MQ=255)
                                      cgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCaa             <  1:5034673/50‑1 (MQ=255)
                                       gcgcCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAat            <  1:1434945/50‑1 (MQ=255)
                                         gcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatat          <  1:3839485/50‑1 (MQ=255)
                                         gcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatat          <  1:7295955/50‑1 (MQ=255)
                                          ccATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatata         <  1:4855612/50‑1 (MQ=255)
                                          ccATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatata         <  1:1111632/50‑1 (MQ=255)
                                           cATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatatat        <  1:7220680/50‑1 (MQ=255)
                                           cATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatatat        >  1:5955331/1‑50 (MQ=255)
                                           cATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatatat        <  1:8174750/50‑1 (MQ=255)
                                           cATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatatat        <  1:3935156/50‑1 (MQ=255)
                                            aTTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAtt       >  1:6627053/1‑50 (MQ=255)
                                            aTTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAtt       >  1:1453364/1‑50 (MQ=255)
                                             ttACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAttt      <  1:1279335/50‑1 (MQ=255)
                                             ttACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAttt      <  1:4455845/50‑1 (MQ=255)
                                             ttACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAttt      <  1:537188/50‑1 (MQ=255)
                                                cGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTCt    <  1:1900935/49‑1 (MQ=255)
                                                 gCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTCTGc  <  1:8037273/50‑1 (MQ=255)
                                                 |                                                 
CAGTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTTCTGC  >  NC_000913/612599‑612697

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