Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,065,341 G→T 10.1% I163I (ATC→ATA ygfI ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,065,3410GT10.1% 220.2 / 6.8 69I163I (ATC→ATAygfIpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/36);  new base T (3/4);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GCTCTTGCTCTTGTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATAAAATCTAA  >  AP009048/3065271‑3065406
                                                                      |                                                                 
gctcttgctcttGTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTg                                                                   <  1:1377457/71‑1 (MQ=255)
gctcttgctcttGTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTGCACTg                                                                   <  1:1352528/71‑1 (MQ=255)
 ctcttgctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGa                                                                  <  1:2279980/71‑1 (MQ=255)
  tcttgctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGa                                                                  >  1:3529125/1‑70 (MQ=255)
   cttgctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGa                                                                  >  1:809566/1‑69 (MQ=255)
   cttgctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGa                                                                  >  1:1913061/1‑69 (MQ=255)
      gctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGt                                                              <  1:3932575/70‑1 (MQ=255)
      gctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGt                                                              >  1:3184130/1‑70 (MQ=255)
      gctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGt                                                              <  1:332175/70‑1 (MQ=255)
      gctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGt                                                              >  1:1203765/1‑70 (MQ=255)
      gctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGt                                                              <  1:1014182/70‑1 (MQ=255)
       ctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTc                                                             <  1:3876210/70‑1 (MQ=255)
       ctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTc                                                             <  1:1681244/70‑1 (MQ=255)
       ctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTc                                                             <  1:3572862/70‑1 (MQ=255)
       ctcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTc                                                             <  1:424258/70‑1 (MQ=255)
        tcttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAg                                                           >  1:4611175/1‑71 (MQ=255)
          ttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGt                                                         >  1:3838615/1‑71 (MQ=255)
          ttgTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGt                                                         >  1:1986024/1‑71 (MQ=255)
            gTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAgg                                                          <  1:1427651/68‑1 (MQ=255)
            gTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAgg                                                          >  1:1104134/1‑68 (MQ=255)
            gTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGttt                                                       >  1:924124/1‑71 (MQ=255)
                ctctTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTc                                                     >  1:1675762/1‑69 (MQ=255)
                  ctTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGt                                                 >  1:2904928/1‑71 (MQ=255)
                   ttCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGt                                                              <  1:3943521/57‑1 (MQ=255)
                    tCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTc                                               >  1:22082/1‑71 (MQ=255)
                     cATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGt                                                         >  1:294553/1‑60 (MQ=255)
                     cATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTcc                                              >  1:1694210/1‑71 (MQ=255)
                      aTAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGa                                                                  <  1:3917386/50‑1 (MQ=255)
                       tAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTg                                            <  1:4295985/71‑1 (MQ=255)
                       tAACGGCGAGAAACGAGAAAGACTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTg                                            <  1:2720846/71‑1 (MQ=255)
                        aaCGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTcc                                                    <  1:3146317/62‑1 (MQ=255)
                        aaCGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTg                                            >  1:1563466/1‑70 (MQ=255)
                          cGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAAc                                         <  1:3601514/71‑1 (MQ=255)
                           ggCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGt                                                              <  1:2383971/49‑1 (MQ=255)
                            gCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTc                                       >  1:3929438/1‑71 (MQ=255)
                              gagaaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGt                                     <  1:1215688/71‑1 (MQ=255)
                              gagaaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGt                                     <  1:3301759/71‑1 (MQ=255)
                               agaaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTg                                    <  1:4144912/71‑1 (MQ=255)
                               agaaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTg                                    >  1:2424210/1‑71 (MQ=255)
                                gaaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTc                                                     >  1:3062958/1‑53 (MQ=255)
                                gaaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTg                                    <  1:1993708/70‑1 (MQ=255)
                                gaaacgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTg                                    >  1:476674/1‑70 (MQ=255)
                                   acgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGa                                                                  >  1:1238111/1‑37 (MQ=255)
                                   acgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTc                                               <  1:2693570/56‑1 (MQ=255)
                                    cgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGa                                >  1:2775108/1‑70 (MQ=255)
                                    cgagaaaGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGt                                                         >  1:4109007/1‑45 (MQ=255)
                                        aaaGCCTAATCGCAGCGGTGGTAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACAt                           >  1:4525425/1‑71 (MQ=255)
                                          aGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGc                         <  1:705165/71‑1 (MQ=255)
                                           gCCTAATCGCAGCGGTTGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCt                        <  1:4586514/71‑1 (MQ=255)
                                           gCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCt                        <  1:2748749/71‑1 (MQ=255)
                                            ccTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTc                                                     <  1:2705631/41‑1 (MQ=255)
                                             cTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAg                                 <  1:3256789/60‑1 (MQ=255)
                                              tAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCtgtg                     >  1:2873045/1‑71 (MQ=255)
                                               aaTCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCGGATACATGCTGTgg                    >  1:3211094/1‑71 (MQ=255)
                                                 tCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTc                                                     <  1:2778509/36‑1 (MQ=255)
                                                 tCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCg                                      <  1:3504381/51‑1 (MQ=255)
                                                          tggtAGTTCACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGataa            <  1:636245/68‑1 (MQ=255)
                                                          tggtAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCtgt                      <  1:3117979/58‑1 (MQ=255)
                                                             tAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGACCTCGTGCAGATACATGCTGTggg                   <  1:1629056/58‑1 (MQ=255)
                                                             tAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCtg                       <  1:1820927/54‑1 (MQ=255)
                                                              aGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATAc                             >  1:1003155/1‑47 (MQ=255)
                                                               gTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTc                                       <  1:4493143/36‑1 (MQ=255)
                                                               gTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATaaa        >  1:2074279/1‑67 (MQ=255)
                                                               gTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATAAAATCt    <  1:3972314/71‑1 (MQ=255)
                                                                ttCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCa                                  <  1:2639324/40‑1 (MQ=255)
                                                                ttCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATa                              <  1:3022154/44‑1 (MQ=255)
                                                                ttCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATAAAATc     >  1:1942308/1‑69 (MQ=255)
                                                                 tCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATaaaa       <  1:4480405/66‑1 (MQ=255)
                                                                  cACTTATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATAAAATCTaa  <  1:1169324/70‑1 (MQ=255)
                                                                    cTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATg                          <  1:3352663/44‑1 (MQ=255)
                                                                    cTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCtgt                      <  1:4372834/48‑1 (MQ=255)
                                                                      |                                                                 
GCTCTTGCTCTTGTCGCTCTTCATAACGGCGAGAAACGAGAAAGCCTAATCGCAGCGGTGGTAGTTCACTGATAGTCAGGTTTTCCAGTTCCTGAACTCGTGCAGATACATGCTGTGGGGAGATAATAAAATCTAA  >  AP009048/3065271‑3065406

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: