Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,609,912 T→G 6.5% E144D (GAA→GAC fhuF ← ferric iron reductase involved in ferric hydroximate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,609,9120TG 6.5% 234.4 / 11.4 92E144D (GAA→GACfhuFferric iron reductase involved in ferric hydroximate transport
Reads supporting (aligned to +/- strand):  ref base T (43/43);  new base G (2/4);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 101 total aligned reads displayed.

GAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCC  >  AP009048/4609845‑4609978
                                                                   |                                                                  
gAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTca                                                                 <  1:1299682/71‑1 (MQ=255)
gAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTca                                                                 >  1:3718481/1‑71 (MQ=255)
gAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTca                                                                 <  1:3831616/71‑1 (MQ=255)
  aCCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTcaca                                                               >  1:5384058/1‑71 (MQ=255)
   ccAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTcacac                                                              <  1:4525907/71‑1 (MQ=255)
     aGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTcacaca                                                             >  1:6466774/1‑70 (MQ=255)
      gcgcCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTcacacac                                                            >  1:3148519/1‑70 (MQ=255)
       cgcCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTcacacaca                                                           >  1:3563841/1‑70 (MQ=255)
        gcCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATc                                                         <  1:1208238/71‑1 (MQ=255)
        gcCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATc                                                         <  1:3184143/71‑1 (MQ=255)
         ccTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTca                                                                 <  1:6150064/62‑1 (MQ=255)
         ccTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTca                                                                 <  1:5355938/62‑1 (MQ=255)
          cTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCt                                                                    <  1:3724004/58‑1 (MQ=255)
          cTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCt                                                                    <  1:1584719/58‑1 (MQ=255)
          cTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCt                                                                    <  1:5011507/58‑1 (MQ=255)
          cTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCt                                                                    <  1:750908/58‑1 (MQ=255)
           tGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCt                                                                    <  1:1755133/57‑1 (MQ=255)
           tGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGAc                                                      >  1:2596594/1‑71 (MQ=255)
            ggCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTcaca                                                               >  1:2518684/1‑61 (MQ=255)
             gCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGAccc                                                    <  1:7100851/71‑1 (MQ=255)
               tGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCtt                                                                   <  1:7460721/54‑1 (MQ=255)
               tGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCtt                                                                   <  1:2341359/54‑1 (MQ=255)
               tGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCtt                                                                   <  1:4246624/54‑1 (MQ=255)
               tGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTc                                                                  <  1:28236/55‑1 (MQ=255)
                gATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCa                                                   >  1:5900560/1‑69 (MQ=255)
                  ttAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATc                                                         <  1:4491983/61‑1 (MQ=255)
                  ttAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATc                                                         <  1:1976580/61‑1 (MQ=255)
                  ttAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATc                                                         <  1:4108040/61‑1 (MQ=255)
                    aaCGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTca                                                                 >  1:5345596/1‑51 (MQ=255)
                    aaCGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGaaa                                               >  1:3966561/1‑69 (MQ=255)
                    aaCGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAAc                                              <  1:4818912/70‑1 (MQ=255)
                     aCGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAg                                            >  1:3769507/1‑71 (MQ=255)
                       gTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCg                                                        <  1:3172289/57‑1 (MQ=255)
                       gTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCg                                                        <  1:3729401/57‑1 (MQ=255)
                       gTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAgg                                           <  1:503090/70‑1 (MQ=255)
                       gTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGc                                          <  1:6373998/71‑1 (MQ=255)
                       gTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGc                                          <  1:1216790/71‑1 (MQ=255)
                       gTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGc                                          <  1:4664274/71‑1 (MQ=255)
                        tttCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAgg                                           <  1:1963853/69‑1 (MQ=255)
                         ttCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCtt                                                                   >  1:7069385/1‑44 (MQ=255)
                         ttCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAgg                                           >  1:1628135/1‑68 (MQ=255)
                         ttCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAgg                                           >  1:2917468/1‑68 (MQ=255)
                         ttCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCAAGAAACAgg                                           >  1:3981574/1‑68 (MQ=255)
                          tCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGAc                                       >  1:7396678/1‑71 (MQ=255)
                           ccATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACAt                                                          >  1:3170592/1‑51 (MQ=255)
                           ccATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATc                                                         <  1:6896804/52‑1 (MQ=255)
                            cATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTcac                                                                >  1:5337382/1‑44 (MQ=255)
                            cATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGAcg                                      <  1:4446002/70‑1 (MQ=255)
                                cGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAgg                                           >  1:4419173/1‑61 (MQ=255)
                                cGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAgg                                           >  1:1597092/1‑61 (MQ=255)
                                cGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTcc                                 >  1:2899865/1‑71 (MQ=255)
                                agATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTc                                  >  1:7146194/2‑70 (MQ=255)
                                 gATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGc                                          <  1:5433765/61‑1 (MQ=255)
                                 gATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTcc                                 <  1:7932280/70‑1 (MQ=255)
                                  aTGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGaaa                                               <  1:2517375/55‑1 (MQ=255)
                                   tgctgcGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGa                                        >  1:6805122/1‑61 (MQ=255)
                                   tgctgcGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGa                                        >  1:6185399/1‑61 (MQ=255)
                                   tgctgcGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGa                                        >  1:5004355/1‑61 (MQ=255)
                                     ctgcGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGttt                            >  1:4518919/1‑71 (MQ=255)
                                     ctgcGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGtt                             <  1:5814660/70‑1 (MQ=255)
                                     ctgcGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGtt                             <  1:2390063/70‑1 (MQ=255)
                                          gcgAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGa                       <  1:3534692/71‑1 (MQ=255)
                                             aaTGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAAc                                              <  1:171017/45‑1 (MQ=255)
                                             aaTGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAAc                                              <  1:271125/45‑1 (MQ=255)
                                               tgtgGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACt                   >  1:1319513/1‑70 (MQ=255)
                                                 tggtgTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGc                                          <  1:7389618/45‑1 (MQ=255)
                                                 tggtgTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGc                                          <  1:7346643/45‑1 (MQ=255)
                                                 tggtgTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGa                                        >  1:1125342/1‑47 (MQ=255)
                                                 tggtgTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAg                >  1:2179439/1‑71 (MQ=255)
                                                   gtgtTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGt                                   <  1:4844943/50‑1 (MQ=255)
                                                   gtgtTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGt                                   <  1:6915520/50‑1 (MQ=255)
                                                   gtgtTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGc               <  1:5696097/70‑1 (MQ=255)
                                                     gtTGCGTTTTTATCTTCACACACATCGACCCAGaaa                                               <  1:4679838/36‑1 (MQ=255)
                                                     gtTGCGTTTTTATCTTCACACACATCGACCCAGaaa                                               <  1:6177513/36‑1 (MQ=255)
                                                      ttGCGTTTTTATCTTCACACACATCGACCCAGAAACAgg                                           >  1:7922548/1‑39 (MQ=255)
                                                      ttGCGTTTTTATCTTCACACACATCGACCCAGAAACAgg                                           >  1:7711778/1‑39 (MQ=255)
                                                      ttGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGa                                        >  1:3562038/1‑42 (MQ=255)
                                                      ttGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGa                                        >  1:429632/1‑42 (MQ=255)
                                                      ttGCGTTTTTATCGTCACACACATCGACCCAGAAAc                                              <  1:5171989/36‑1 (MQ=255)
                                                      ttGCGTTTTTATCGTCACACACATCGACCCAGAAAc                                              <  1:4666019/36‑1 (MQ=255)
                                                      ttGCGTTTTTATCGTCACACACATCGACCCAGAAAc                                              <  1:5259782/36‑1 (MQ=255)
                                                      ttGCGTTTTTATCGTCACACACATCGACCCAGAAAc                                              <  1:3962464/36‑1 (MQ=255)
                                                       tccGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGt                                   >  1:6756205/3‑46 (MQ=255)
                                                        gCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTc                                  <  1:3930536/46‑1 (MQ=255)
                                                        gCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTc                                  <  1:2823272/46‑1 (MQ=255)
                                                        gCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACt                   >  1:3860695/1‑61 (MQ=255)
                                                        gCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATgg           >  1:124750/1‑69 (MQ=255)
                                                         cGTTTTTATCTTCACACTCATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGaa         <  1:2754343/70‑1 (MQ=255)
                                                         cGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTc                                  <  1:7060184/45‑1 (MQ=255)
                                                          gTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCg                                >  1:7733266/1‑46 (MQ=255)
                                                          gTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAAt       >  1:3535062/1‑71 (MQ=255)
                                                           tttttATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATg      >  1:467168/1‑71 (MQ=255)
                                                           tttttATCGTCACACACATCGACCCAGAAACAGGCGGCACGGCCGGTTTCGTGAAACTCAGCAt             >  1:4313486/1‑64 (MQ=255)
                                                           tttttATCGTCACACACATCGACCCAGAAACAGGCGGCACGGCCGGTTTCGTGAAACTCAGCAt             >  1:5798674/1‑64 (MQ=255)
                                                            ttttATCTTCACACACATCGACCCAGAAACAGGCg                                         >  1:6661058/1‑35 (MQ=255)
                                                             tttATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGtt    >  1:1915113/1‑71 (MQ=255)
                                                             tttATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGtt    >  1:4894513/1‑71 (MQ=255)
                                                             tttATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGt     <  1:1308101/70‑1 (MQ=255)
                                                              ttATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGtt                             <  1:3448621/45‑1 (MQ=255)
                                                              ttATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTc   >  1:6510551/1‑71 (MQ=255)
                                                                   |                                                                  
GAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCC  >  AP009048/4609845‑4609978

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