Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 609,476 A→G 5.5% intergenic (‑165/+1) entD ← / ← fepA phosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/iron‑enterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048609,4760AG 5.5% 422.1 / 11.8 146intergenic (‑165/+1)entD/fepAphosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/iron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (67/71);  new base G (5/3);  total (72/74)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 152 total aligned reads displayed.

GGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCG  >  AP009048/609407‑609546
                                                                     |                                                                      
ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCAt                                                                       <  1:4467367/71‑1 (MQ=255)
ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCAt                                                                       >  1:3662887/1‑71 (MQ=255)
 gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCAt                                                                       >  1:2214797/1‑70 (MQ=255)
 gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATc                                                                      <  1:6263926/71‑1 (MQ=255)
 gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATc                                                                      <  1:217419/71‑1 (MQ=255)
  ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCAt                                                                       >  1:4119157/1‑69 (MQ=255)
     gATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATcgg                                                                    >  1:4416089/1‑67 (MQ=255)
     gATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGaa                                                                  <  1:586060/71‑1 (MQ=255)
       tAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGt                                                                <  1:5034899/71‑1 (MQ=255)
        aaGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                                                               >  1:6451364/1‑71 (MQ=255)
        aaGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                                                               >  1:2358680/1‑71 (MQ=255)
           gCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgt                                                            >  1:1223834/1‑71 (MQ=255)
            cGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGaa                                                                  >  1:6922265/1‑64 (MQ=255)
             gTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAg                                                                    <  1:2178768/61‑1 (MQ=255)
             gTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgtgt                                                          >  1:7474475/1‑71 (MQ=255)
               tCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                         >  1:1908663/1‑70 (MQ=255)
                cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                        >  1:2336040/1‑70 (MQ=255)
                cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                         >  1:2026877/1‑69 (MQ=255)
                  cgccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACg                                                     >  1:1718169/1‑71 (MQ=255)
                   gccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGc                                                    <  1:7654689/71‑1 (MQ=255)
                    ccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAg                                                                    <  1:1963804/54‑1 (MQ=255)
                     cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCt                                                   <  1:4021075/70‑1 (MQ=255)
                       cATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGTTGTTTACGCTCat                                                <  1:4239660/71‑1 (MQ=255)
                       cATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCat                                                <  1:3817595/71‑1 (MQ=255)
                       cATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCat                                                <  1:369017/71‑1 (MQ=255)
                        aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCat                                                <  1:4559861/70‑1 (MQ=255)
                         tCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                 <  1:5316977/68‑1 (MQ=255)
                         tCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCa                                                 <  1:5631879/68‑1 (MQ=255)
                          ccGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTa                                                       >  1:4833292/1‑61 (MQ=255)
                          ccGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTa                                                       >  1:7610779/1‑61 (MQ=255)
                           cGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCatat                                              >  1:4351272/1‑69 (MQ=255)
                              cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                        <  1:5056891/56‑1 (MQ=255)
                              cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                        <  1:3068111/56‑1 (MQ=255)
                              cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                         <  1:528128/71‑1 (MQ=255)
                              cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                         <  1:3423054/71‑1 (MQ=255)
                                tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGa                                                                   <  1:3238108/43‑1 (MQ=255)
                                tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGa                                                                   <  1:329285/43‑1 (MQ=255)
                                tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGa                                                                   <  1:1456338/43‑1 (MQ=255)
                                 gAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                         <  1:5664479/52‑1 (MQ=255)
                                  aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCa                                                                     >  1:1603819/1‑39 (MQ=255)
                                  aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCa                                                                     >  1:4123520/1‑39 (MQ=255)
                                  aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCa                                                                     >  1:2854344/1‑39 (MQ=255)
                                  aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                         <  1:3159263/51‑1 (MQ=255)
                                   acgaggCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTgg                                                              >  1:1769419/1‑45 (MQ=255)
                                   acgacgCGCACTTTGTCAACAATCTGACGTTAGCATCAGaa                                                                  >  1:5731396/1‑41 (MQ=255)
                                   acgacgCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTgg                                                              >  1:6992482/1‑45 (MQ=255)
                                   acgacgCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTgg                                                              >  1:4927997/1‑45 (MQ=255)
                                    cgacgCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgtacgt                                   >  1:110934/1‑71 (MQ=255)
                                     gacgCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTc                                  <  1:3271129/71‑1 (MQ=255)
                                        gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCg                               <  1:2465694/71‑1 (MQ=255)
                                        gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCg                               <  1:5412467/71‑1 (MQ=255)
                                         cgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                              <  1:5534423/71‑1 (MQ=255)
                                          gcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGc                             <  1:3011917/71‑1 (MQ=255)
                                           cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTAc                                                      <  1:6953943/45‑1 (MQ=255)
                                           cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTAc                                                      <  1:925556/45‑1 (MQ=255)
                                            aCTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgtgt                                                          >  1:6388943/1‑40 (MQ=255)
                                             cTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTAc                                                      <  1:707033/43‑1 (MQ=255)
                                             cTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCt                                                   <  1:680149/46‑1 (MQ=255)
                                             cTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCg                          <  1:1309926/71‑1 (MQ=255)
                                             cTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCg                          >  1:7366537/1‑71 (MQ=255)
                                              tttGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacg                                        <  1:4084180/56‑1 (MQ=255)
                                               ttGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTAc                                                      >  1:5555082/1‑41 (MQ=255)
                                                tGTCAACAATCTGACGTTAGCATCAGAAGTGGgtg                                                           >  1:3002776/1‑35 (MQ=255)
                                                tGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGc                                                    <  1:4191500/42‑1 (MQ=255)
                                                tGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacg                                        >  1:4598344/1‑54 (MQ=255)
                                                tGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacg                                        >  1:7589959/1‑54 (MQ=255)
                                                tGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGtt                        <  1:7372040/70‑1 (MQ=255)
                                                 gTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTta                       <  1:6970738/70‑1 (MQ=255)
                                                  tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGc                                                    <  1:4844260/40‑1 (MQ=255)
                                                  tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtat                      <  1:5388897/70‑1 (MQ=255)
                                                   caacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCt                                                   >  1:1302467/1‑40 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCg                               <  1:1379949/59‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg                   <  1:2904530/71‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg                   <  1:2022847/71‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg                    <  1:4172284/70‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg                    <  1:1472404/70‑1 (MQ=255)
                                                       aaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                         >  1:2278659/1‑46 (MQ=255)
                                                       aaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                         >  1:5460437/1‑46 (MQ=255)
                                                       aaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                         >  1:3237159/1‑46 (MQ=255)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCacgta                                      >  1:7492396/2‑48 (MQ=21)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCacgta                                      >  1:7257973/2‑48 (MQ=21)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCacgta                                      >  1:6642135/2‑48 (MQ=21)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTc                                  >  1:7937092/2‑52 (MQ=25)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCg                          <  1:7248046/59‑1 (MQ=37)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCg                          <  1:3411010/59‑1 (MQ=37)
                                                         tCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg                    >  1:775177/1‑65 (MQ=255)
                                                         tCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg                    >  1:2593560/1‑65 (MQ=255)
                                                           tGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                              <  1:7089220/53‑1 (MQ=255)
                                                            gACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGcc           >  1:1596888/1‑71 (MQ=255)
                                                             aCGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGc            <  1:7773303/69‑1 (MQ=255)
                                                              cGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGcccc         >  1:813759/1‑71 (MQ=255)
                                                              cGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGcccc         >  1:1416213/1‑71 (MQ=255)
                                                              cGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGcccc         >  1:4299405/1‑71 (MQ=255)
                                                               gTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg                   >  1:7690807/1‑60 (MQ=255)
                                                               gTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg                   >  1:6324707/1‑60 (MQ=255)
                                                                     aTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                                 >  1:5651371/1‑40 (MQ=255)
                                                                     aTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGc                             <  1:2386621/44‑1 (MQ=255)
                                                                     aTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGc                             <  1:3890546/44‑1 (MQ=255)
                                                                     aTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCAcc   <  1:4670438/70‑1 (MQ=255)
                                                                     aTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCg  >  1:6016703/1‑71 (MQ=255)
                                                                     |                                                                      
GGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCG  >  AP009048/609407‑609546

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