Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 657,802 T→A 8.6% F30Y (TTC→TAC)  cspE → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048657,8020TA 8.6% 255.6 / 14.0 93F30Y (TTC→TAC) cspEDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (45/40);  new base A (6/2);  total (51/42)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGA  >  AP009048/657734‑657871
                                                                    |                                                                     
cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                                     <  1:4821509/71‑1 (MQ=255)
cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                                     <  1:4435674/71‑1 (MQ=255)
cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                                     <  1:3431917/71‑1 (MQ=255)
 gTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                                     >  1:3369993/1‑70 (MQ=255)
   tAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa                                                                   >  1:3956411/1‑70 (MQ=255)
   tAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa                                                                   >  1:2956402/1‑70 (MQ=255)
    aaGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTac                                                                  <  1:4326880/70‑1 (MQ=255)
      gTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACt                                                               >  1:4092656/1‑71 (MQ=255)
         gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct                                                            <  1:7325972/71‑1 (MQ=255)
         gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct                                                            <  1:1870735/71‑1 (MQ=255)
         gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGAAGTGTACGAACACGtct                                                            <  1:6192711/71‑1 (MQ=25)
         gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGAAGTGTACGAACACGtct                                                            <  1:1673916/71‑1 (MQ=25)
            tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTacac                                                                >  1:3305081/1‑64 (MQ=255)
            tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg                                                         >  1:3708827/1‑71 (MQ=255)
            tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg                                                         >  1:2044733/1‑71 (MQ=255)
            tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg                                                         >  1:2713369/1‑71 (MQ=255)
                 aGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                    <  1:1736395/71‑1 (MQ=255)
                  gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                    >  1:7063905/1‑70 (MQ=255)
                  gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                    >  1:7035422/1‑70 (MQ=255)
                  gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                    >  1:5668175/1‑70 (MQ=255)
                   tCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                    >  1:1791814/1‑69 (MQ=255)
                   tCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa                                                  <  1:2587352/71‑1 (MQ=255)
                    ccAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                    >  1:5165310/1‑68 (MQ=255)
                     cAAAGTATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc                                                           <  1:5409624/60‑1 (MQ=255)
                     cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTc                                                                      >  1:4542572/1‑49 (MQ=255)
                     cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc                                                           <  1:2329395/60‑1 (MQ=255)
                      aaaGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGAc                                               >  1:7511701/1‑71 (MQ=255)
                        aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                                     >  1:6676690/1‑47 (MQ=255)
                        aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa                                                                   >  1:6912680/1‑49 (MQ=255)
                         ggATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTac                                                                  <  1:5762591/49‑1 (MQ=255)
                          gATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAAt                                           <  1:5627647/71‑1 (MQ=255)
                           aTTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                    <  1:6614943/61‑1 (MQ=255)
                             tCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGt                                        >  1:3941755/1‑71 (MQ=255)
                              cGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGt                                        <  1:7902918/70‑1 (MQ=255)
                               ggTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct                                                            >  1:2361193/1‑49 (MQ=255)
                               ggTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAAt                                           <  1:3558143/66‑1 (MQ=255)
                                gTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                                     <  1:2895236/39‑1 (MQ=255)
                                gTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGtttt                                     <  1:371197/71‑1 (MQ=255)
                                 tttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                                      >  1:4961038/1‑69 (MQ=255)
                                  ttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa                                                  >  1:805729/1‑56 (MQ=255)
                                  ttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa                                                  >  1:6116692/1‑56 (MQ=255)
                                   tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGtttt                                     >  1:1441973/1‑68 (MQ=255)
                                   tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                                   >  1:5787280/1‑70 (MQ=255)
                                   tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                                   >  1:5446529/1‑70 (MQ=255)
                                    cATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa                                 >  1:2539389/1‑71 (MQ=255)
                                     aTTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTaca                                                                 >  1:7744232/1‑38 (MQ=255)
                                     aTTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa                                 <  1:3504256/70‑1 (MQ=255)
                                     aTTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAc                                >  1:1308675/1‑71 (MQ=255)
                                      ttACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGAc                                               >  1:1283595/1‑55 (MQ=255)
                                      ttACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAct                               >  1:5507068/1‑71 (MQ=255)
                                        aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTaca                                                                 >  1:4893129/1‑35 (MQ=255)
                                        aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAActc                              <  1:6520011/70‑1 (MQ=255)
                                        aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAActc                              <  1:2718746/70‑1 (MQ=255)
                                         cTCCTGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAc                                >  1:1057902/1‑67 (MQ=255)
                                          tCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc                                                           <  1:1127338/39‑1 (MQ=255)
                                            cGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCa                                                       >  1:2385637/1‑41 (MQ=255)
                                             ggAAGACGGCAGCAAAGACGTGTTCTTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                                   <  1:364742/60‑1 (MQ=255)
                                             ggAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                                   <  1:2888888/60‑1 (MQ=255)
                                             ggAAGACGGAAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg                        <  1:7858808/71‑1 (MQ=255)
                                              gAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                    >  1:2234123/1‑42 (MQ=255)
                                              gAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa                       <  1:3585061/71‑1 (MQ=255)
                                               aaGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGa                                                >  1:2249522/1‑45 (MQ=255)
                                               aaGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg                        >  1:3603404/1‑69 (MQ=255)
                                                aGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg                        >  1:6733109/1‑68 (MQ=255)
                                                 gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa                       <  1:1332503/68‑1 (MQ=255)
                                                 gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa                       <  1:6065605/68‑1 (MQ=255)
                                                 gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa                       <  1:6056808/68‑1 (MQ=255)
                                                 gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAgg                    >  1:1862221/1‑71 (MQ=255)
                                                  aCGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGt                   <  1:6885882/71‑1 (MQ=255)
                                                   cGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGAc                                               >  1:47654/1‑42 (MQ=255)
                                                   cGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAg                     >  1:5650000/1‑68 (MQ=255)
                                                   cGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTc                  <  1:5249713/71‑1 (MQ=255)
                                                    ggCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGt                   >  1:4800274/1‑69 (MQ=255)
                                                    ggCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCa                 <  1:7136193/71‑1 (MQ=255)
                                                      cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTa                                    <  1:6929864/50‑1 (MQ=255)
                                                      cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgc               <  1:4612432/71‑1 (MQ=255)
                                                      cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgc               <  1:5720084/71‑1 (MQ=255)
                                                       agcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa                                 >  1:3387973/1‑52 (MQ=255)
                                                            aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                                      >  1:5670939/1‑42 (MQ=255)
                                                            aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                                      >  1:4568527/1‑42 (MQ=255)
                                                            aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                                      >  1:2666809/1‑42 (MQ=255)
                                                            aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAct                               <  1:317109/49‑1 (MQ=255)
                                                            aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg              >  1:7816748/1‑66 (MQ=25)
                                                            aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg              >  1:497230/1‑66 (MQ=25)
                                                            aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg              >  1:5940894/1‑66 (MQ=25)
                                                            aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg              >  1:487866/1‑66 (MQ=25)
                                                            aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGcgcg            >  1:1229900/1‑68 (MQ=25)
                                                            aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTCCTGAAGGTCAGcgcg            >  1:6715600/1‑68 (MQ=25)
                                                                gtgtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg                        <  1:7376272/52‑1 (MQ=255)
                                                                gtgtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg                        <  1:2039999/52‑1 (MQ=255)
                                                                  gtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCt                         <  1:2857452/49‑1 (MQ=255)
                                                                  gtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCt                         <  1:6076128/49‑1 (MQ=255)
                                                                  gtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGtt     <  1:3086961/69‑1 (MQ=255)
                                                                   ttCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCg   <  1:3209608/70‑1 (MQ=255)
                                                                    tCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaa                                  >  1:64293/1‑38 (MQ=255)
                                                                    tCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaa                                  >  1:536051/1‑38 (MQ=255)
                                                                    tCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaa                                  >  1:5052724/1‑38 (MQ=255)
                                                                    tCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGa  >  1:6842191/1‑70 (MQ=255)
                                                                    |                                                                     
CGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGA  >  AP009048/657734‑657871

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: