Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,094,198 A→T 9.2% T301S (ACA→TCA)  ycdT → predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,094,1980AT 9.2% 237.0 / 14.5 87T301S (ACA→TCA) ycdTpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (43/36);  new base T (6/2);  total (49/38)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAA  >  AP009048/1094130‑1094265
                                                                    |                                                                   
cAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGa                                                                     >  1:7866283/1‑69 (MQ=255)
 aaaGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAg                                                                  <  1:4025253/71‑1 (MQ=255)
   aGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGa                                                                     >  1:5927794/1‑66 (MQ=255)
   aGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGtt                                                                >  1:3684331/1‑71 (MQ=255)
   aGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGtt                                                                >  1:5369196/1‑71 (MQ=255)
    gAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGtt                                                                >  1:5260626/1‑70 (MQ=255)
     aaCATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCa                                                              <  1:3242929/71‑1 (MQ=255)
        aTTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAatc                                                           <  1:1265467/71‑1 (MQ=255)
        aTTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAat                                                            >  1:7800826/1‑70 (MQ=255)
              caTCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGTCAGTTCAatcagcatc                                                     <  1:4189390/71‑1 (MQ=255)
              caTCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagcatc                                                     <  1:7906710/71‑1 (MQ=255)
              caTCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagcatc                                                     <  1:2337804/71‑1 (MQ=255)
                tCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagcatcag                                                   >  1:5306470/1‑71 (MQ=255)
                tCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagcatcag                                                   <  1:612254/71‑1 (MQ=255)
                 cgcgATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagcatcagc                                                  >  1:3708425/1‑71 (MQ=255)
                 cgcgATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagcatcagc                                                  >  1:838388/1‑71 (MQ=255)
                     aTCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGtt                                                                >  1:4274133/1‑53 (MQ=255)
                     aTCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGtt                                                                >  1:5204653/1‑53 (MQ=255)
                     aTCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaa                                              <  1:990511/71‑1 (MQ=255)
                     aTCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaa                                              <  1:1250018/71‑1 (MQ=255)
                     aTCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaa                                              <  1:2403571/71‑1 (MQ=255)
                      tCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaa                                             <  1:5640047/71‑1 (MQ=255)
                      tCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaa                                             <  1:2588026/71‑1 (MQ=255)
                        ccTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagcatc                                                     <  1:4000253/61‑1 (MQ=255)
                        ccTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaaa                                           <  1:2791749/71‑1 (MQ=255)
                        ccTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaa                                            <  1:6235451/70‑1 (MQ=255)
                        ccTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaa                                            >  1:337556/1‑70 (MQ=255)
                         cTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCa                                                              <  1:5554411/51‑1 (MQ=255)
                           tAACCAATATATTTAAAAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagcatcagc                                                  >  1:2107656/1‑61 (MQ=255)
                            aaCCAATATATTTAACAGAAATTATTTTTTTAATGAACTGa                                                                     >  1:1537387/1‑41 (MQ=255)
                            aaCCAATATATTTAACAGAAATTATTTTTTTAATGAACTGa                                                                     >  1:1836672/1‑41 (MQ=255)
                              ccAATATATTTAACAGAAATTATTTTTTTATTGAACTGTCAGTTCat                                                             <  1:3166086/47‑2 (MQ=37)
                              ccAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCt                                     >  1:51759/1‑71 (MQ=255)
                                 atatatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagc                                                        >  1:5630021/1‑49 (MQ=255)
                                 atatatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaaaa                                          >  1:1048937/1‑63 (MQ=255)
                                  tatatTTAACAGAAATTATTTTTTTAATGAACTGAc                                                                    <  1:4980670/36‑1 (MQ=255)
                                  tatatTTAACAGAAATTATTTTTTTAATGAACTGAc                                                                    <  1:7753353/36‑1 (MQ=255)
                                  tatatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaa                                            <  1:3244611/60‑1 (MQ=255)
                                  tatatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTAt                                  <  1:1983620/70‑1 (MQ=255)
                                   atatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAat                                                            >  1:3636422/1‑43 (MQ=255)
                                    tatTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTg                                <  1:2868335/70‑1 (MQ=255)
                                       ttAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCaaaa                                            <  1:3708407/55‑1 (MQ=255)
                                        tAACAGAAATTATTTTTTTAATGAACTGACAGTTCAatc                                                           >  1:4813536/1‑39 (MQ=255)
                                        tAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTc                            <  1:6024227/70‑1 (MQ=255)
                                        tAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCa                           >  1:7075211/1‑71 (MQ=255)
                                        tAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCa                           >  1:7221252/1‑71 (MQ=255)
                                        tAACAGAAATTATTTTTTTAATGAACAGACAGTTCAatc                                                           >  1:4048177/1‑39 (MQ=255)
                                          aCAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagcatcagc                                                  >  1:3509921/1‑46 (MQ=255)
                                          aCAGAAATTATTTTTTTAATGAACTGACAGTTCAatcagcatcagc                                                  >  1:1117763/1‑46 (MQ=255)
                                          aCAGAAATTATTTTTTTAATGAACTGACACTTCAatcagcatcagc                                                  >  1:336916/1‑46 (MQ=255)
                                           cAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTa                                   >  1:1648284/1‑60 (MQ=255)
                                           cAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCAt                          >  1:3427392/1‑69 (MQ=255)
                                           cAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCAt                          >  1:2421295/1‑69 (MQ=255)
                                            aGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGAt                       <  1:4430023/71‑1 (MQ=255)
                                             gAAATTATTTTTTTAATGAACTGACAGTTCAATCTGCATCAGCCCAAAAAACGCCTTATTGCGTCAt                          >  1:1447152/1‑67 (MQ=255)
                                             gAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGAtt                      <  1:799077/71‑1 (MQ=255)
                                             gAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGAt                       >  1:636854/1‑70 (MQ=255)
                                              aaaTTATTTTTTTAATGAACTGACAGTTCAatcagcatc                                                     <  1:7617896/39‑1 (MQ=255)
                                              aaaTTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTa                     >  1:2825981/1‑71 (MQ=255)
                                               aaTTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGAtt                      >  1:7646978/1‑69 (MQ=255)
                                               aaTTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTAt                    <  1:2878635/71‑1 (MQ=255)
                                                attattTTTTTAATGAACTGACAGTTCAATCAGCATCAGccc                                                >  1:7516445/1‑42 (MQ=255)
                                                attattTTTTTAATGAACTGACAGTTCAATCAGCATCAGccc                                                >  1:6547890/1‑42 (MQ=255)
                                                attattTTTTTAATGAACTGACAGTTCAATCAGCATCAGccc                                                >  1:2014942/1‑42 (MQ=255)
                                                attattTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTAt                                  >  1:5501405/1‑56 (MQ=255)
                                                    tttttttAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTAtt                                 >  1:46895/1‑53 (MQ=255)
                                                    tttttttAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGt                             <  1:4560220/57‑1 (MQ=255)
                                                    tttttttAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGt                             <  1:7146458/57‑1 (MQ=255)
                                                     ttttttAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTa                     >  1:3050583/1‑64 (MQ=255)
                                                     ttttttAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTa                     >  1:6057511/1‑64 (MQ=255)
                                                     ttttttAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTa                     >  1:3138186/1‑64 (MQ=255)
                                                     ttttttAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGata               >  1:5799972/1‑70 (MQ=255)
                                                       ttttAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATc             <  1:5744661/70‑1 (MQ=255)
                                                        tttAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATc             <  1:5917833/69‑1 (MQ=255)
                                                        tttAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGa           >  1:2778738/1‑71 (MQ=255)
                                                          tAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGAc          <  1:5936622/70‑1 (MQ=255)
                                                          tAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGAc          <  1:1999700/70‑1 (MQ=255)
                                                          tAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGAc          >  1:7287901/1‑70 (MQ=255)
                                                            ttGAACTGTCAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGt                             >  1:6116108/2‑49 (MQ=255)
                                                            ttGAACTGTCAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGt                             >  1:1096298/2‑49 (MQ=255)
                                                            ttGAACTGTCAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGt                             >  1:3926085/2‑49 (MQ=255)
                                                            ttGAACTGTCAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGt                             >  1:6451035/2‑49 (MQ=255)
                                                            ttGAACTGTCAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGt                             >  1:4267788/2‑49 (MQ=255)
                                                            ttGAACTGTCAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGt                             >  1:5075095/2‑49 (MQ=255)
                                                             tGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCAt                          >  1:1081552/1‑51 (MQ=255)
                                                               aaCTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTAttc                                >  1:3169044/1‑42 (MQ=255)
                                                               aaCTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTc    <  1:2538985/71‑1 (MQ=255)
                                                                aCTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTa                     >  1:3562882/1‑53 (MQ=255)
                                                                aCTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACtt     <  1:4408803/69‑1 (MQ=255)
                                                                 cTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCaa  <  1:4121473/71‑1 (MQ=255)
                                                                  tGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCaa  <  1:3186875/70‑1 (MQ=255)
                                                                    |                                                                   
CAAAGAACATTGCACATCGCGATCCCTTAACCAATATATTTAACAGAAATTATTTTTTTAATGAACTGACAGTTCAATCAGCATCAGCCCAAAAAACGCCTTATTGCGTCATGATTATGGATATCGACCACTTCAA  >  AP009048/1094130‑1094265

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: