Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 609,469 C→T 5.5% intergenic (‑158/+8) 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,4690CT 5.5% 475.9 / 10.7 146intergenic (‑158/+8)entD/fepAphosphopantetheinyltransferase component of enterobactin synthase multienzyme complex/iron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (62/76);  new base T (5/3);  total (67/79)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 149 total aligned reads displayed.

TAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCC  >  AP009048/609405‑609539
                                                                |                                                                      
  ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCAt                                                                <  1:4467367/71‑1 (MQ=255)
   gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCAt                                                                >  1:2214797/1‑70 (MQ=255)
   gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATc                                                               <  1:217419/71‑1 (MQ=255)
   gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATc                                                               <  1:6263926/71‑1 (MQ=255)
    ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCAt                                                                >  1:4119157/1‑69 (MQ=255)
       gATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGc                                                                  >  1:7656646/1‑64 (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)
               gTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGt                                                                      >  1:6252568/1‑52 (MQ=255)
               gTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAg                                                             <  1:2178768/61‑1 (MQ=255)
               gTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgtgt                                                   >  1:7474475/1‑71 (MQ=255)
                ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGtt                                                                     <  1:5919367/52‑1 (MQ=255)
                ttCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGtt                                                                     <  1:6380892/52‑1 (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:6969354/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)
                       cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCa                                                              <  1:4488259/52‑1 (MQ=255)
                       cgcATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCt                                            <  1:4021075/70‑1 (MQ=255)
                         cATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGTTGTTTACGCTCat                                         <  1:4239660/71‑1 (MQ=255)
                         cATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCat                                         <  1:369017/71‑1 (MQ=255)
                          aTCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCata                                        >  1:6872651/1‑71 (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)
                             cGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTggg                                                      >  1:2516240/1‑54 (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:3423054/71‑1 (MQ=255)
                                cATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCac                                  <  1:528128/71‑1 (MQ=255)
                                  tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGa                                                            <  1:329285/43‑1 (MQ=255)
                                  tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGa                                                            <  1:3238108/43‑1 (MQ=255)
                                  tGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGc                                             >  1:7731972/1‑58 (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)
                                     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)
                                        acgCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTcc                          <  1:7034687/71‑1 (MQ=255)
                                          gcgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCg                        <  1:5412467/71‑1 (MQ=255)
                                           cgcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                       <  1:5534423/71‑1 (MQ=255)
                                            gcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                       >  1:4400363/1‑70 (MQ=255)
                                            gcACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGc                      <  1:3011917/71‑1 (MQ=255)
                                             cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTAc                                               <  1:925556/45‑1 (MQ=255)
                                             cACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTAc                                               <  1:6953943/45‑1 (MQ=255)
                                              aCTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGgtgt                                                   >  1:6388943/1‑40 (MQ=255)
                                               cTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCt                                            <  1:680149/46‑1 (MQ=255)
                                               cTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCg                   <  1:1309926/71‑1 (MQ=255)
                                               cTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCg                   >  1:7366537/1‑71 (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:7589959/1‑54 (MQ=255)
                                                  tGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGtt                 <  1:7372040/70‑1 (MQ=255)
                                                  tGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTta                >  1:6634635/1‑71 (MQ=255)
                                                    tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGc                                             <  1:4844260/40‑1 (MQ=255)
                                                    tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCg                        >  1:890839/1‑61 (MQ=255)
                                                    tCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTtat               <  1:5388897/70‑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:150906/1‑46 (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:7257973/2‑48 (MQ=21)
                                                          ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCacgta                               >  1:7492396/2‑48 (MQ=21)
                                                          ttCTGATGCTAACGTCAGAAGTGGGTGTTTACGCTCATATACCacgta                               >  1:5666961/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:1262546/59‑1 (MQ=37)
                                                           tCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg             >  1:2593560/1‑65 (MQ=255)
                                                           tCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAg             >  1:775177/1‑65 (MQ=255)
                                                           tCTGACGTTAGCATAAGAAGTGGGTGTTTACGCTCATATAc                                     >  1:7903829/1‑41 (MQ=255)
                                                             tGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                       <  1:7089220/53‑1 (MQ=255)
                                                             tGACGTCAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCgg                       <  1:7104968/53‑1 (MQ=255)
                                                                cGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGcccc  >  1:1416213/1‑71 (MQ=255)
                                                                cGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGcccc  >  1:4299405/1‑71 (MQ=255)
                                                                |                                                                      
TAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTTGTCAACAATCTGACGTTAGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCC  >  AP009048/609405‑609539

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: