Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 281,469 A→T 7.3% intergenic (‑262/‑12) yagA ← / → yagE predicted DNA‑binding transcriptional regulator/predicted lyase/synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048281,4690AT 7.3% 268.2 / 11.2 96intergenic (‑262/‑12)yagA/yagEpredicted DNA‑binding transcriptional regulator/predicted lyase/synthase
Reads supporting (aligned to +/- strand):  ref base A (45/44);  new base T (3/4);  total (48/48)
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.60e-01

CCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCC  >  AP009048/281405‑281538
                                                                |                                                                     
ccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAg                                                                  >  1:3757659/1‑70 (MQ=255)
ccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAg                                                                  >  1:490412/1‑70 (MQ=255)
ccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAg                                                                  >  1:6126176/1‑70 (MQ=255)
ccGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAg                                                                  >  1:6255713/1‑70 (MQ=255)
 cGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATATTATAGaa                                                                <  1:6575941/71‑1 (MQ=255)
  gATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGaa                                                                >  1:6917864/1‑70 (MQ=255)
    tttACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACa                                                              >  1:7639474/1‑70 (MQ=255)
    tttACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTACATTATAGaa                                                                >  1:6536296/1‑68 (MQ=255)
    tttACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTACATTATAGaa                                                                >  1:2424635/1‑68 (MQ=255)
    tttACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTACATTATAGaa                                                                >  1:2430532/1‑68 (MQ=255)
     ttACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATaa                                                                       >  1:730736/1‑60 (MQ=255)
       aCCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACaaaa                                                           <  1:5751288/70‑1 (MQ=255)
         ctctCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGt                                                                    >  1:6079022/1‑59 (MQ=255)
         ctctCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAAt                                                          <  1:4607699/69‑1 (MQ=255)
         ctctCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGa                                                        <  1:5333116/71‑1 (MQ=255)
           ctcCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTa                                                                   >  1:432553/1‑58 (MQ=255)
              cTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGAtt                                                      <  1:4509334/68‑1 (MQ=255)
                 gTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGc                                                  <  1:5005513/69‑1 (MQ=255)
                       cGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCATCAAGgag                                            <  1:1087739/69‑1 (MQ=255)
                       cGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGgaga                                           <  1:7039784/70‑1 (MQ=255)
                       cGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGgag                                            <  1:4181812/69‑1 (MQ=255)
                        gTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAtc                                         <  1:4137231/71‑1 (MQ=255)
                         tCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATg                                                                     <  1:6118683/42‑1 (MQ=255)
                         tCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAg                                                   <  1:4265104/60‑1 (MQ=255)
                         tCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAg                                                   <  1:4119151/60‑1 (MQ=255)
                         tCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGgag                                            <  1:2692638/67‑1 (MQ=255)
                         tCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGgag                                            <  1:4179382/67‑1 (MQ=255)
                         tCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAtct                                        >  1:6153583/1‑71 (MQ=255)
                         tCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAtct                                        <  1:1064177/71‑1 (MQ=255)
                         tCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAtc                                         >  1:6273007/1‑70 (MQ=255)
                           aTAGTATGATCGTTAAATAAACGAACTCTGTTCTATAATGTAGAACAAAATGATTCAGCa                                                 <  1:4873981/60‑1 (MQ=255)
                           aTAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCa                                                 <  1:453463/60‑1 (MQ=255)
                           aTAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCa                                                 >  1:672493/1‑60 (MQ=255)
                           aTAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCa                                                 <  1:609944/60‑1 (MQ=255)
                           aTAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCa                                      >  1:543205/1‑71 (MQ=255)
                           aTAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCa                                      >  1:6965878/1‑71 (MQ=255)
                           aTAGTATGATCGATAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAtctc                                       >  1:2609032/1‑70 (MQ=255)
                           aTAGTAGGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCa                                      >  1:5341/1‑71 (MQ=255)
                            tAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGt                                                                    <  1:2848051/40‑1 (MQ=255)
                            tAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTa                                                                   <  1:5558990/41‑1 (MQ=255)
                             aGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGa                                                        <  1:4856174/51‑1 (MQ=255)
                              gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTc                                                     >  1:5283960/1‑53 (MQ=255)
                              gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGc                                   >  1:3071101/1‑71 (MQ=255)
                              gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGc                                   >  1:4381938/1‑71 (MQ=255)
                               tATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGcc                                  >  1:2538096/1‑71 (MQ=255)
                               tATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGc                                   <  1:6776212/70‑1 (MQ=255)
                                aTGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAAt                                                          <  1:3023817/46‑1 (MQ=255)
                                 tGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCg                                 >  1:5475826/1‑70 (MQ=255)
                                   aTCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGAtt                                                      >  1:5987454/1‑47 (MQ=255)
                                   aTCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCa                               <  1:5862311/70‑1 (MQ=255)
                                     cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAAt                                                          >  1:5366599/1‑41 (MQ=255)
                                     cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTc                            >  1:1950095/1‑71 (MQ=255)
                                     cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTc                            >  1:3251992/1‑71 (MQ=255)
                                     cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTc                            >  1:3372525/1‑71 (MQ=255)
                                      gTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAgg                                              <  1:5223138/52‑1 (MQ=255)
                                       ttAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATg                                                         <  1:6492225/40‑1 (MQ=255)
                                        taaataaaCGAACGCTGTTCTATAATGTAGAACAAAAt                                                          <  1:6641321/38‑1 (MQ=255)
                                         aaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCcgcg                        >  1:2345891/1‑71 (MQ=255)
                                         aaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCcgc                         <  1:3693098/70‑1 (MQ=255)
                                          aataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt                       >  1:5974070/1‑71 (MQ=255)
                                           ataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCa                                                    <  1:445809/41‑1 (MQ=255)
                                           ataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCa                                                    <  1:7811676/41‑1 (MQ=255)
                                           ataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAtctc                                       <  1:6313564/54‑1 (MQ=255)
                                           ataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAtctc                                       <  1:4152298/54‑1 (MQ=255)
                                           ataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAtctc                                       <  1:1868853/54‑1 (MQ=255)
                                           ataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAtctc                                       <  1:7882778/54‑1 (MQ=255)
                                           ataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAtctc                                       <  1:4483953/54‑1 (MQ=255)
                                             aaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCaa                                                >  1:2179177/1‑43 (MQ=255)
                                             aaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAg                                               >  1:452925/1‑44 (MQ=255)
                                             aaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGgag                                            >  1:5851716/1‑47 (MQ=255)
                                             aaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGgag                                            >  1:4163657/1‑47 (MQ=255)
                                              aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCaa                                                <  1:2981770/42‑1 (MQ=255)
                                              aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCaa                                                <  1:1401791/42‑1 (MQ=255)
                                              aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgttgtt                   <  1:2661578/71‑1 (MQ=255)
                                                 gaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAc                >  1:2511477/1‑71 (MQ=255)
                                                   acgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgttgt                    >  1:501378/1‑65 (MQ=255)
                                                    ttttGTTCTACATTATAGAACAGAATGATTCAGCAAGGAGATCTCATGCCGCAGTCcgcg                        <  1:6542623/57‑1 (MQ=25)
                                                    ttttGTTCTACATTATAGAACAGAATGATTCAGCAAGGAGATCTCATGCCGCAGTCcgcg                        <  1:1620353/57‑1 (MQ=25)
                                                    cgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACgg              >  1:7927700/1‑70 (MQ=255)
                                                     tttGTTCTACATTATAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTcc                           <  1:4026281/54‑1 (MQ=37)
                                                     gCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAt                                          >  1:1861016/1‑41 (MQ=255)
                                                     gCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAt                                          >  1:945408/1‑41 (MQ=255)
                                                     gCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAt                                          >  1:7684649/1‑41 (MQ=255)
                                                     gCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAt                                          >  1:4977244/1‑41 (MQ=255)
                                                     gCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAt                                          >  1:6179362/1‑41 (MQ=255)
                                                     gCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGa            <  1:2400836/71‑1 (MQ=255)
                                                       tGTTCTATAATGTAGAACAAAATGATTCAGCAAGga                                             <  1:11211/36‑1 (MQ=255)
                                                       tGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgttgtt                   <  1:4116765/62‑1 (MQ=255)
                                                       tGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAAt          >  1:727059/1‑71 (MQ=255)
                                                         ttCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGc                                <  1:3221803/47‑1 (MQ=255)
                                                         ttCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATc         <  1:7899816/70‑1 (MQ=255)
                                                          tCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCCCGGGAATCAt       >  1:1087706/1‑71 (MQ=255)
                                                           cTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCAtt      >  1:7548570/1‑71 (MQ=255)
                                                            tataATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCcgcg                        <  1:6784368/52‑1 (MQ=255)
                                                             ataATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAAt          >  1:7749231/1‑65 (MQ=255)
                                                              taATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGaa           >  1:4460215/1‑63 (MQ=255)
                                                              taATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTccc   >  1:1393273/1‑71 (MQ=255)
                                                               aaTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACgg              >  1:6108386/1‑59 (MQ=255)
                                                               aaTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACgg              >  1:1013926/1‑59 (MQ=255)
                                                               aaTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTcccc  >  1:6474076/1‑71 (MQ=255)
                                                                |                                                                     
CCGATTTACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCC  >  AP009048/281405‑281538

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: