Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 609,474 G→A 5.4% intergenic (‑163/+3) 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,4740GA 5.4% 486.7 / 10.8 147intergenic (‑163/+3)entD/fepAphosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/iron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (68/71);  new base A (5/3);  total (73/74)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 148 total aligned reads displayed.

GGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGG  >  AP009048/609407‑609541
                                                                   |                                                                   
ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCAt                                                                  <  1:4467367/71‑1 (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)
     gATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGc                                                                    >  1:7656646/1‑64 (MQ=255)
     gATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATcgg                                                               >  1:4416089/1‑67 (MQ=255)
     gATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGaa                                                             <  1:586060/71‑1 (MQ=255)
       tAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGt                                                           <  1:5034899/71‑1 (MQ=255)
       tAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGACAACAATCTGACGTTAGCATCAGAAGt                                                           <  1:1454242/71‑1 (MQ=255)
        aaGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTg                                                          >  1:1323080/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)
              ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGa                                                              <  1:4637372/61‑1 (MQ=255)
               tCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                    >  1:1908663/1‑70 (MQ=255)
                catgccggATGCGGCGTGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTa                                                  >  1:3800433/4‑71 (MQ=37)
                cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGttt                                                   >  1:2336040/1‑70 (MQ=255)
                cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                    >  1:2026877/1‑69 (MQ=255)
                cACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTa                                                  >  1:2946457/1‑71 (MQ=255)
                  cgccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACg                                                >  1:5025603/1‑71 (MQ=255)
                  cgccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACg                                                >  1:1718169/1‑71 (MQ=255)
                   gccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGc                                               <  1:7654689/71‑1 (MQ=255)
                    ccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAg                                                               <  1:1963804/54‑1 (MQ=255)
                    ccgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAg                                                               <  1:4237727/54‑1 (MQ=255)
                     cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCa                                                                <  1:4488259/52‑1 (MQ=255)
                     cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCa                                                                <  1:5541681/52‑1 (MQ=255)
                     cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCa                                                                <  1:5272035/52‑1 (MQ=255)
                     cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCt                                              <  1:4021075/70‑1 (MQ=255)
                       cATCCGGCATGAACGACGCTCACTTTGTCAACAATCTGACGATAGCATCAGAAGTGGGTGTTTACGCTCat                                           <  1:5640832/71‑1 (MQ=38)
                       cATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGTTGTTTACGCTCat                                           <  1:4239660/71‑1 (MQ=255)
                       cATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCat                                           <  1:369017/71‑1 (MQ=255)
                          ccGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTa                                                  >  1:7610779/1‑61 (MQ=255)
                           cGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTggg                                                        >  1:2516240/1‑54 (MQ=255)
                           cGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTa                                                  >  1:3253834/1‑60 (MQ=255)
                           cGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCatat                                         >  1:4351272/1‑69 (MQ=255)
                              cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                    <  1:3423054/71‑1 (MQ=255)
                              cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                    <  1:528128/71‑1 (MQ=255)
                                tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGa                                                              <  1:3238108/43‑1 (MQ=255)
                                tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGa                                                              <  1:329285/43‑1 (MQ=255)
                                tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGa                                                              <  1:5321793/43‑1 (MQ=255)
                                tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGc                                               >  1:7731972/1‑58 (MQ=255)
                                tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacg                                   >  1:6348942/1‑70 (MQ=255)
                                tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacg                                   >  1:416687/1‑70 (MQ=255)
                                 gAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGtt                                                    <  1:5664479/52‑1 (MQ=255)
                                  aaCGACGCGCACTTTGTCAACAATCTGACGTTAGCATCa                                                                >  1:1603819/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:3206132/1‑45 (MQ=255)
                                   acgacgCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTgg                                                         >  1:6992482/1‑45 (MQ=255)
                                    cgacgCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgtacgt                              >  1:110934/1‑71 (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:7366537/1‑71 (MQ=255)
                                             cTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCg                     <  1:1309926/71‑1 (MQ=255)
                                               ttGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTAc                                                 >  1:5555082/1‑41 (MQ=255)
                                                tGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacg                                   >  1:4598344/1‑54 (MQ=255)
                                                tGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacg                                   >  1:7589959/1‑54 (MQ=255)
                                                tGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGtt                   <  1:7372040/70‑1 (MQ=255)
                                                tGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTta                  >  1:6634635/1‑71 (MQ=255)
                                                 gTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTta                  <  1:6970738/70‑1 (MQ=255)
                                                  tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGc                                               <  1:4844260/40‑1 (MQ=255)
                                                  tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCg                          >  1:890839/1‑61 (MQ=255)
                                                  tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtat                 <  1:4302317/70‑1 (MQ=255)
                                                  tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtat                 <  1:5388897/70‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCg                          <  1:3039766/59‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCg                          <  1:4836489/59‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCg                          <  1:3989879/59‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg              <  1:2904530/71‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg              <  1:2022847/71‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg               <  1:1472404/70‑1 (MQ=255)
                                                    aacaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg               <  1:4172284/70‑1 (MQ=255)
                                                     acaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgta                                 >  1:3543146/1‑51 (MQ=255)
                                                     acaaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCacgta                                 >  1:165879/1‑51 (MQ=255)
                                                       aaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                    >  1:2278659/1‑46 (MQ=255)
                                                       aaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                    >  1:5460437/1‑46 (MQ=255)
                                                       aaTCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                    >  1:150906/1‑46 (MQ=255)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCacgta                                 >  1:7492396/2‑48 (MQ=21)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCacgta                                 >  1:6642135/2‑48 (MQ=21)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCacgta                                 >  1:7257973/2‑48 (MQ=21)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTc                             >  1:7937092/2‑52 (MQ=25)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCg                     <  1:1262546/59‑1 (MQ=37)
                                                        ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCg                     <  1:7248046/59‑1 (MQ=37)
                                                         tCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg               >  1:2593560/1‑65 (MQ=255)
                                                         tCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg               >  1:775177/1‑65 (MQ=255)
                                                         tCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg               >  1:5424951/1‑65 (MQ=255)
                                                         tCTGACGTTAGCATAAGAAGTGGGTGTTTACGCTCATATAc                                       >  1:7903829/1‑41 (MQ=255)
                                                           tGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                         <  1:7089220/53‑1 (MQ=255)
                                                              cGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGcccc    >  1:1416213/1‑71 (MQ=255)
                                                               gTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAgg              >  1:4107978/1‑60 (MQ=255)
                                                               gTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCg   >  1:6270301/1‑71 (MQ=255)
                                                                 tAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCTGCTCGTTATAGGTATACGCCCCgg  >  1:3315960/1‑70 (MQ=255)
                                                                   |                                                                   
GGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGG  >  AP009048/609407‑609541

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: