Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,688,553 A→T 7.8% V401D (GTC→GAC)  yfhK ← predicted sensory kinase in two‑component system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,688,5530AT 7.8% 141.5 / 4.4 51V401D (GTC→GAC) yfhKpredicted sensory kinase in two‑component system
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base T (2/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTG  >  AP009048/2688487‑2688623
                                                                  |                                                                      
cccTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACAt                                                                     >  1:1523368/1‑70 (MQ=255)
 ccTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATc                                                                    <  1:1179184/70‑1 (MQ=255)
  cTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCaa                                                                  <  1:1703745/71‑1 (MQ=255)
     aaaaaaGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAata                                                                >  1:354780/1‑70 (MQ=255)
       aaaaGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATaa                                                             <  1:2131916/71‑1 (MQ=255)
         aaGGGTTCGAAGATCATGTCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAatat                                                               <  1:127500/67‑1 (MQ=255)
         aaGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAAc                                                           >  1:3315279/1‑71 (MQ=255)
             gTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCcgc                                                       <  1:4048543/71‑1 (MQ=255)
                cGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAAcc                                                          >  1:1988788/1‑65 (MQ=255)
                 gAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAAtata                                                              >  1:1897172/1‑60 (MQ=255)
                    gATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATaa                                                             <  1:112099/58‑1 (MQ=255)
                     aTCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAat                                                                 <  1:909903/53‑1 (MQ=255)
                     aTCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAat                                                                 <  1:1439000/53‑1 (MQ=255)
                       cATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGt                                              >  1:4522996/1‑70 (MQ=255)
                          ggCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACAt                                                                     >  1:2850628/1‑44 (MQ=255)
                          ggCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAAc                                          >  1:2587014/1‑71 (MQ=255)
                             gcgTTCCTCTTGCGGAATGGGCGTGCCTATATTGATGTCATCAata                                                                <  1:3557183/46‑1 (MQ=37)
                             gcgTTCCTCTTGCGGAATGGGCGTGCCTATATTGATGTCATCAata                                                                <  1:3551969/46‑1 (MQ=37)
                              cgTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAActgct                                      >  1:2718241/1‑71 (MQ=255)
                                ttCCTCTTGCGGAATGGGCGTGCCTATATTGATGTCATCAatac                                                               >  1:862129/1‑43 (MQ=25)
                                  ccTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCcgcg                                                      >  1:342960/1‑51 (MQ=255)
                                  ccTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTa                                             <  1:1418824/60‑1 (MQ=255)
                                  ccTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCg                                   >  1:719259/1‑70 (MQ=255)
                                  ccTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCg                                   >  1:1535988/1‑70 (MQ=255)
                                   ctctTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATaa                                                             <  1:3215536/43‑1 (MQ=255)
                                    tctTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAAcc                                                          <  1:3097530/45‑1 (MQ=255)
                                      ttGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAgg                               >  1:363317/1‑70 (MQ=255)
                                      ttGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAg                                <  1:1564290/69‑1 (MQ=255)
                                       tGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGc                              >  1:3877049/1‑70 (MQ=255)
                                        gCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAAcc                                                          >  1:1871865/1‑41 (MQ=255)
                                        gCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCa                             <  1:73166/70‑1 (MQ=255)
                                         cGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCTCACCATGTAAACTGCTGCGAAGGCaa                            <  1:540950/70‑1 (MQ=255)
                                         cGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCaaa                           >  1:3711723/1‑71 (MQ=255)
                                         cGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCaa                            <  1:4452783/70‑1 (MQ=255)
                                          ggAATGGGCGTGCCTGTATTGATGACATCAATATAAAc                                                           >  1:954586/1‑38 (MQ=255)
                                          ggAATGGGCGTCCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACtgctgc                                    <  1:3645332/61‑1 (MQ=255)
                                             aTGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGtt                       <  1:3350244/71‑1 (MQ=255)
                                             aTGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGt                        >  1:1570752/1‑70 (MQ=255)
                                              tGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTa                      >  1:4106331/1‑71 (MQ=255)
                                                 gCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGa                                  <  1:2675859/56‑1 (MQ=255)
                                                  cGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTaaa                                           >  1:887739/1‑46 (MQ=255)
                                                   gTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAg                                >  1:73379/1‑56 (MQ=255)
                                                    tGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCg                                   <  1:3664415/52‑1 (MQ=255)
                                                    tGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGAt                >  1:1784830/1‑71 (MQ=255)
                                                      ccTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGaa                                 >  1:4059537/1‑52 (MQ=255)
                                                      ccTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGAtt               <  1:1194562/70‑1 (MQ=255)
                                                        tGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTCCTGCGAAGGCAAATGtt                       <  1:3184120/60‑1 (MQ=255)
                                                         gTATTGATGACATCAATATAAACCCGCGCACCATGTAAActgctg                                     >  1:3203697/1‑45 (MQ=255)
                                                         gTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAAt                          >  1:1383302/1‑56 (MQ=255)
                                                         atATTGATGTCATCAATACAAACCCGCGCACCATGTaaa                                           >  1:1104406/2‑39 (MQ=21)
                                                           aTTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGAtt               <  1:4447881/65‑1 (MQ=255)
                                                              gatgaCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGaa                                 >  1:3824487/1‑44 (MQ=255)
                                                               atgaCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGAt                <  1:473301/60‑1 (MQ=255)
                                                                tgaCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGa                 <  1:3108881/58‑1 (MQ=255)
                                                                  aCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTg  <  1:1240541/71‑1 (MQ=255)
                                                                  |                                                                      
CCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTG  >  AP009048/2688487‑2688623

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: