Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 19,348 G→A 16.2% G335D (GGT→GAT)  USA300HOU_0014 → DHH subfamily phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073019,3480GA16.2% 103.7 / 13.2 49G335D (GGT→GAT) USA300HOU_0014DHH subfamily phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (26/15);  new base A (3/5);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCAACATTAC  >  CP000730/19209‑19488
                                                                                                                                           |                                                                                                                                            
cgcgGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGt                                                                                                                                             <  1:351565/141‑1 (MQ=255)
  cgGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGa                                                                                                                                           >  1:27365/1‑141 (MQ=255)
        cGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATCTCCTTGCAGCGGGTGACAAAGt                                                                                                                                     >  2:114064/1‑141 (MQ=255)
          aCCAAGTTTCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCa                                                                                                                                   <  1:257643/141‑1 (MQ=255)
            cAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCAtt                                                                                                                                 <  1:121629/141‑1 (MQ=255)
             aaGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCTCATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTa                                                                                                                                <  1:19193/141‑1 (MQ=255)
               gTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATACTTGGAGAGGGAGACAAAGTCCTTATc                                                                                                                              >  2:554758/1‑141 (MQ=255)
                     aTTAAAAGTATTAATGGTAAAGTGCGTTTCTATGGCGGTAAGACGGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAAAGGGTGACAAGGTCATTATCATGGGa                                                                                                                        >  1:347044/1‑141 (MQ=255)
                               ttAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTc                                                                                                              >  1:254316/1‑141 (MQ=255)
                                  aTGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAAAGCACGTGTGATCTCACATGCATCAAAAGAGATCCTCGCAAATGGGGACAAAGTTCTTAGCATCGGACATCAACGTCCTg                                                                                                           >  2:457812/1‑141 (MQ=255)
                                     gTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACt                                                                                                        >  1:251525/1‑141 (MQ=255)
                                             cGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTCATCGTGTAAGAGCACGTGTGACCACACATGCGGTAAAAGATAACCTTGCAGAGGGTGAAAAAGACATTATCATGGGCCAAAAACGTCCTAACTTataagcc                                                                                                >  2:220362/1‑135 (MQ=255)
                                                   tATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGCCATAAACGTCCTGACTTAGCTGCAATTGGt                                                                                          >  1:36534/1‑141 (MQ=255)
                                                    aTGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTg                                                                                         >  1:100561/1‑141 (MQ=255)
                                                     tGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTgc                                                                                        >  1:451604/1‑141 (MQ=255)
                                                       gCGGTAAGACTGACCCGATGGATAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAACAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTgcag                                                                                      <  1:22606/141‑1 (MQ=255)
                                                       gCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACAGGCGTTAAAAGATATCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTgcag                                                                                      >  1:242229/1‑141 (MQ=255)
                                                       gCGGCAAGACGGTCCCGATTGCGAAACGCTCTCGTGTAAGAGTATGTGTGACCTCACATTCGTTAAATTATTTCCTTGGAGAGGGTGACAATGTCAGTATCATGGGACATAAACGTCCTGACTTAGATGCGATTGGTgtag                                                                                      <  2:142902/141‑1 (MQ=255)
                                                        gtgTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTgcagc                                                                                     <  1:375797/139‑1 (MQ=255)
                                                          gTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCaa                                                                                   <  2:100561/141‑1 (MQ=255)
                                                           tAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATCGGTGCAGCAAt                                                                                  >  2:336781/1‑141 (MQ=255)
                                                           tAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAAt                                                                                  <  1:440463/141‑1 (MQ=255)
                                                             aGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCg                                                                                >  2:320887/1‑141 (MQ=255)
                                                              gACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCgg                                                                               >  2:395005/1‑141 (MQ=255)
                                                                cTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGgtg                                                                             <  2:287558/141‑1 (MQ=255)
                                                                  gACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACCTAGATGCAATTGGTGCAGCAATCGgtgtg                                                                           <  2:308510/141‑1 (MQ=255)
                                                                             gaAACGGACTCGTGTAAGAGCACGTGTGATCTCACATGCATTAAAAGATATCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGc                                                                <  1:46683/141‑1 (MQ=255)
                                                                                  gTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAatga                                                           >  1:328736/1‑141 (MQ=255)
                                                                                  gTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAatga                                                           <  1:213983/141‑1 (MQ=255)
                                                                                        gtgtAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGaata                                                     >  1:513228/1‑141 (MQ=255)
                                                                                          gtAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGaataat                                                   >  1:523213/1‑141 (MQ=255)
                                                                                                 cACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGAGGAATAATTTAGAAg                                            <  1:299852/141‑1 (MQ=255)
                                                                                                   cGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCa                                          >  1:330190/1‑141 (MQ=255)
                                                                                                    gtgtGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCAt                                         >  2:393642/1‑141 (MQ=255)
                                                                                                       tgATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACa                                      >  1:52673/1‑141 (MQ=255)
                                                                                                       tgATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACa                                      >  2:479172/1‑141 (MQ=255)
                                                                                                       tgATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACa                                      >  1:324420/1‑141 (MQ=255)
                                                                                                             cacaTGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTCCCATCGGACCTAAACGGCCATACTTAGATGCAAGTGGGGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTAt                                <  1:164615/141‑1 (MQ=255)
                                                                                                               caTGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGCTGCAATGGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTa                              >  2:487196/1‑141 (MQ=255)
                                                                                                                     ttAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAATCATACATCGTATTAAATgag                        >  1:17277/1‑141 (MQ=255)
                                                                                                                       aaaaGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATGTAGAAGCATACATCGTATTAAATGAGAc                      >  1:303990/1‑141 (MQ=255)
                                                                                                                            atatCCTTGCAGAGTATCACAAATTCATTATCATGGTACATAAACGTCCTGACTTAGATGCAATTGCTGCAACAATCGGTCTGTCTACATTTACAACGATGCATAAATTAAAACCATACATCGTATTAAATgcaaccaaca                 <  2:564855/141‑10 (MQ=255)
                                                                                                                            atatCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACa                 <  1:471021/141‑1 (MQ=255)
                                                                                                                             tatCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACAt                >  1:249482/1‑141 (MQ=255)
                                                                                                                              atCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAACCGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACAtt               <  2:293477/141‑1 (MQ=255)
                                                                                                                                 cTTGCAGAGGGTGACAAAGTCGTTATGATGGGACATAAACGTCCTGACTTAGATGCAAGTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGAt            <  2:435867/141‑1 (MQ=255)
                                                                                                                                  ttGCAGAGGGTTACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATc           <  1:263576/141‑1 (MQ=255)
                                                                                                                                   tGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGGCTGCCATTGATcc          >  1:206088/1‑141 (MQ=255)
                                                                                                                                     cAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCaa        >  1:325296/1‑141 (MQ=255)
                                                                                                                                         gggTGACAAGGTCATTATCATTGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCAACAtt    <  2:514071/141‑1 (MQ=255)
                                                                                                                                          ggTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCAACATTa   <  1:320887/141‑1 (MQ=255)
                                                                                                                                           gTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCAACATTac  >  2:230191/1‑141 (MQ=255)
                                                                                                                                           gTGACAAAGTCATTATCATGGGACATAAACGTCCGGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAAGAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCAACATTac  <  2:50747/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
CGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCAACATTAC  >  CP000730/19209‑19488

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: