Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,233,032 C→A 7.6% intergenic (+80/+11) yjbB → / ← pepE predicted transporter/(alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,233,0320CA 7.6% 255.1 / 4.3 78intergenic (+80/+11)yjbB/pepEpredicted transporter/(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (44/28);  new base A (2/4);  total (46/32)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAG  >  AP009048/4232967‑4233101
                                                                  |                                                                     
tAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCaa                                                                  >  1:5579175/1‑71 (MQ=255)
  ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAAtt                                                                >  1:3184792/1‑71 (MQ=255)
  ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAAtt                                                                >  1:2920410/1‑71 (MQ=255)
   gCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTc                                                               <  1:1060183/71‑1 (MQ=255)
    ccTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCg                                                                      >  1:3779133/1‑63 (MQ=255)
    ccTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAAtt                                                                <  1:1974841/69‑1 (MQ=255)
      tGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACt                                                            >  1:616697/1‑71 (MQ=255)
        aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTc                                                               <  1:3490948/66‑1 (MQ=255)
        aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACtt                                                           >  1:1363010/1‑70 (MQ=255)
        aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACt                                                            <  1:3125710/69‑1 (MQ=255)
        aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACt                                                            >  1:7645225/1‑69 (MQ=255)
           aGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTaaaa                                                       <  1:228958/71‑1 (MQ=255)
           aGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTaaaa                                                       <  1:7801562/71‑1 (MQ=255)
           aGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTaaa                                                        <  1:5953320/70‑1 (MQ=255)
              cgcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACg                                                    <  1:5485127/71‑1 (MQ=255)
              cgcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACg                                                    <  1:1735210/71‑1 (MQ=255)
              cgcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACg                                                    <  1:7239894/71‑1 (MQ=255)
               gcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACgg                                                   <  1:7575328/71‑1 (MQ=255)
               gcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCAgg                                                                     <  1:3309026/53‑1 (MQ=255)
               gcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCAgg                                                                     <  1:4842987/53‑1 (MQ=255)
               gcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCAgg                                                                     <  1:1986592/53‑1 (MQ=255)
               gcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCAgg                                                                     <  1:7314295/53‑1 (MQ=255)
                cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCgg                                                                     <  1:248578/52‑1 (MQ=18)
                 gTCAAGCGTCGCATCAGGGACTATGAACAAGATGCCGGATCATGCATCAGGCa                                                                   >  1:5545544/1‑53 (MQ=37)
                 gTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCAGGCa                                                                   >  1:7273669/1‑53 (MQ=255)
                  tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATcc                                                                       <  1:5927415/48‑1 (MQ=255)
                   cAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGAc                                               >  1:5121009/1‑71 (MQ=255)
                   cAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGAc                                               >  1:2410217/1‑71 (MQ=255)
                    aaGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCaa                                                                  <  1:1186634/51‑1 (MQ=255)
                    aaGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACg                                                    <  1:37096/65‑1 (MQ=255)
                    aaGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGa                                                >  1:945209/1‑69 (MQ=255)
                     aGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCa                                             <  1:3758003/71‑1 (MQ=255)
                      gCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGCCCAg                                            >  1:4653460/1‑71 (MQ=255)
                      gCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAg                                            <  1:3141356/71‑1 (MQ=255)
                      gCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAg                                            >  1:7427982/1‑71 (MQ=255)
                         tCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTc                                                               <  1:5497408/49‑1 (MQ=255)
                         tCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCa                                             <  1:7346768/67‑1 (MQ=255)
                         tCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCtt                                         <  1:4871179/71‑1 (MQ=255)
                          cGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTc                                        >  1:3942660/1‑71 (MQ=255)
                          cGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTc                                        >  1:551897/1‑71 (MQ=255)
                           gCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACtt                                                           >  1:1668234/1‑51 (MQ=255)
                           gCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACtt                                                           >  1:854368/1‑51 (MQ=255)
                            cATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCa                                                              >  1:7570305/1‑47 (MQ=255)
                               cAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGc                                    >  1:6657413/1‑70 (MQ=255)
                               cAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCg                                   >  1:913962/1‑71 (MQ=255)
                                   cACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAAc                               >  1:4151023/1‑71 (MQ=255)
                                    aCTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAAcc                              >  1:5989640/1‑71 (MQ=255)
                                        tGAACAAGATGCCGGATCATGCATCCGGCAATTCAc                                                             >  1:3287961/1‑36 (MQ=255)
                                        tGAACAAGATGCCGGATCATGCATCCGGCAATTCAc                                                             >  1:2850664/1‑36 (MQ=255)
                                         gAACAAGATGCCGGATCATGCATCCGGCAATTCAc                                                             >  1:1521097/1‑35 (MQ=255)
                                         gAACAAGATGCCGGATCATGCATCCGGCAATTCAc                                                             >  1:7199933/1‑35 (MQ=255)
                                         gAACAAGATGCCGGATCATGCATCCGGCAATTCAc                                                             >  1:7191188/1‑35 (MQ=255)
                                         gAACAAGATGCCGGATCATGCATCCGGCAATTCAc                                                             >  1:6670754/1‑35 (MQ=255)
                                         gAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                         >  1:1130121/1‑71 (MQ=255)
                                          aaCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACgg                                                   >  1:4220177/1‑44 (MQ=255)
                                          aaCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCg                                   >  1:4498948/1‑60 (MQ=255)
                                           aCAAGATGCCGGATCATGCATCCGGCAATTCACTTaaaaa                                                      <  1:1995980/40‑1 (MQ=255)
                                           aCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttctt                       <  1:7482496/71‑1 (MQ=255)
                                           aCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttctt                       <  1:3161962/71‑1 (MQ=255)
                                           aCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttct                        >  1:6102334/1‑70 (MQ=255)
                                            cAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAg                                            >  1:935167/1‑49 (MQ=255)
                                               gATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttct                        >  1:1585079/1‑66 (MQ=255)
                                                 tGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAg                 <  1:6343939/71‑1 (MQ=255)
                                                  gCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGc                >  1:6916156/1‑71 (MQ=255)
                                                   ccGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGcc               >  1:3383136/1‑71 (MQ=255)
                                                   ccGGATCATCCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGcc               >  1:4148253/1‑71 (MQ=255)
                                                     ggATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAAc                               >  1:3462586/1‑53 (MQ=255)
                                                      gATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCtt             >  1:7234389/1‑70 (MQ=255)
                                                      gATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTa            <  1:4422874/71‑1 (MQ=255)
                                                           aTGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAac         <  1:6229175/70‑1 (MQ=255)
                                                              cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCa                                      >  1:86132/1‑38 (MQ=255)
                                                              cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCa                                      >  1:3987361/1‑38 (MQ=255)
                                                              cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAAcaca      >  1:4067573/1‑70 (MQ=255)
                                                              cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACAt     <  1:2865334/71‑1 (MQ=255)
                                                               aTCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGa                                 <  1:1459011/42‑1 (MQ=255)
                                                               aTCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACAt     >  1:618495/1‑70 (MQ=255)
                                                               aTCCGGCAATTAACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATa    <  1:7889711/71‑1 (MQ=255)
                                                                tCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCt              >  1:2314834/1‑60 (MQ=255)
                                                                tCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACAt     >  1:5479272/1‑69 (MQ=255)
                                                                 ccGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGc                            >  1:4766842/1‑45 (MQ=255)
                                                                  cGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg  <  1:5442079/70‑1 (MQ=255)
                                                                  cGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg  <  1:7749759/70‑1 (MQ=255)
                                                                  cGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg  <  1:2608468/70‑1 (MQ=255)
                                                                  cGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg  <  1:5742215/70‑1 (MQ=255)
                                                                  |                                                                     
TAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAG  >  AP009048/4232967‑4233101

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: