Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 875,787 T→C 6.9% R151R (CGT→CGC yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,7870TC 6.9% 162.5 / 3.7 58R151R (CGT→CGCyliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (29/25);  new base C (2/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGTGATGAATATAATGCTAAAAACGCCA  >  NC_000913/875653‑875926
                                                                                                                                      |                                                                                                                                           
gTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCt                                                                                                                                        >  2:73377/1‑140 (MQ=255)
gTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCt                                                                                                                                        >  1:1321150/1‑140 (MQ=255)
   aaGGAGTTCGAGAAAAAGTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTAACCGTTAACTCACAATCGTATTCTGat                                                                                                                                     >  1:496866/1‑140 (MQ=255)
         ttCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATAttttt                                                                                                                               >  1:1658920/1‑140 (MQ=255)
          tCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTg                                                                                                                              >  2:1604717/1‑140 (MQ=255)
              gtaaaaTTATTGTTCATATATTATGCAGGTTGTGTTGTATCTAAATGTCAATAATTCTGGGTTTTTCTCTTTTGTCGGTAAAACCTTCGTCTATTACCTCTATCCGGTTACTCACAATCGTAGTCTGATATTTTTGCTTg                                                                                                                          <  1:1418217/138‑1 (MQ=255)
              gaAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTg                                                                                                                          <  1:1418816/140‑1 (MQ=255)
                   attattGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTg                                                                                                                     <  1:862715/140‑1 (MQ=255)
                            cAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAAAACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTATGCTTGGTCTGGAGCGattt                                                                                                            >  1:291352/1‑140 (MQ=255)
                                         ggTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAATCCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCaa                                                                                               >  1:1664182/1‑140 (MQ=255)
                                           ttGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCaaaa                                                                                             >  1:797790/1‑140 (MQ=255)
                                               gAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGc                                                                                         <  2:932511/140‑1 (MQ=255)
                                               gAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGc                                                                                         <  2:586459/140‑1 (MQ=255)
                                               gAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGc                                                                                         <  2:287510/140‑1 (MQ=255)
                                               gAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGc                                                                                         <  2:1649405/140‑1 (MQ=255)
                                                     tCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGa                                                                                   >  1:295862/1‑140 (MQ=255)
                                                          aTGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGGTGGCGATGGAc                                                                              >  2:1240783/1‑140 (MQ=255)
                                                             tCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGc                                                                           >  2:1307435/1‑140 (MQ=255)
                                                             tCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGc                                                                           >  1:1452264/1‑140 (MQ=255)
                                                              cAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCg                                                                          <  2:1151560/140‑1 (MQ=255)
                                                              cAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCg                                                                          <  1:650419/140‑1 (MQ=255)
                                                                       aGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGa                                                                 >  2:848557/1‑140 (MQ=255)
                                                                       aGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGAGAGAACCTGa                                                                 >  2:712278/1‑140 (MQ=255)
                                                                             tATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTTACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTtct                                                           <  1:529331/140‑1 (MQ=255)
                                                                             tATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTtct                                                           <  2:1658920/140‑1 (MQ=255)
                                                                              aTCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTtctc                                                          >  1:640810/1‑140 (MQ=255)
                                                                                ctctTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGCtgtctc                                                        <  1:1612704/140‑1 (MQ=255)
                                                                                     ttGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTAttt                                                   >  1:1379364/1‑140 (MQ=255)
                                                                                     ttGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTAttt                                                   >  1:1468465/1‑140 (MQ=255)
                                                                                     ttGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTAttt                                                   >  1:1379605/1‑140 (MQ=255)
                                                                                     ttGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTAttt                                                   >  1:1234894/1‑140 (MQ=255)
                                                                                          ggTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGaa                                              <  1:848557/140‑1 (MQ=255)
                                                                                             aaaaCCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTtctct                                                         <  1:805160/126‑1 (MQ=255)
                                                                                             aaaaCCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTtctct                                                         >  2:805160/1‑126 (MQ=255)
                                                                                                    tCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGGTCTCTCTATTTAAGAACGGTAAaccg                                    >  1:615739/1‑140 (MQ=255)
                                                                                                           ttACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAaccggtgaccg                             <  2:1417816/140‑1 (MQ=255)
                                                                                                              ccTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAaccggtgaccggtg                          <  2:1067339/140‑1 (MQ=255)
                                                                                                              ccTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAaccggtgaccggtg                          <  2:442320/140‑1 (MQ=255)
                                                                                                              ccTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAaccggtgaccggtg                          <  2:1004128/140‑1 (MQ=255)
                                                                                                              ccTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAaccggtgaccggtg                          <  2:1028724/140‑1 (MQ=255)
                                                                                                              ccTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGATTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAaccggtgaccggtg                          <  2:442479/140‑1 (MQ=255)
                                                                                                                tctATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGtgat                        <  1:1252063/140‑1 (MQ=255)
                                                                                                                tctATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGtgat                        <  1:1561698/140‑1 (MQ=255)
                                                                                                                 ctATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGtgatg                       >  2:1409477/1‑140 (MQ=255)
                                                                                                                  tATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGtgatga                      >  2:551667/1‑140 (MQ=255)
                                                                                                                   aTCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGTGATGaa                     >  2:1582564/1‑140 (MQ=255)
                                                                                                                   aTCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGTGATGaa                     >  2:306274/1‑140 (MQ=255)
                                                                                                                    tCCGGTAACTCACAATCGCAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAaccggtgaccg                             >  2:309791/1‑131 (MQ=255)
                                                                                                                    tCCGGTAACTCACAATCGCAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAaccggtgaccg                             >  2:309643/1‑131 (MQ=255)
                                                                                                                    tCCGGTAACTCACAATCGCAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAaccggtgaccg                             <  1:309643/131‑1 (MQ=255)
                                                                                                                    tCCGGTAACTCACAATCGCAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAaccggtgaccg                             <  1:309791/131‑1 (MQ=255)
                                                                                                                     ccGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGTGATGAata                   >  1:1517194/1‑140 (MQ=255)
                                                                                                                       tgTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGTGATGAAtata                 <  2:555036/139‑1 (MQ=255)
                                                                                                                         tAACTCACAAACGTAGTCTGATATTTAAGCTTCGTCTGGAGAGTTTATCTTTACTGTCAAAATCGCTGGCGATGGAAAGCGAGAACCAGATCATCTCTCGAATAAAGAACGGTAAACCGGTGACCGGTGATGAATATAAt               >  1:648810/1‑140 (MQ=255)
                                                                                                                            cTCACAATCGTAGTCTGAGATTTTTGCTTGGTCTGGAGCGTTTTTCGTTACTGTCAAAATCGCTGGCGAGGGACAGCGAGAACCTGATGTTCTCTCTCTTTAAGAACGGTAAACCGGTGACCGGTGATGAAGATAATGCt            >  1:680321/1‑140 (MQ=255)
                                                                                                                              cacaATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGTGATGAATATAATGCTaa          >  1:411330/1‑140 (MQ=255)
                                                                                                                                    cGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGTGATGAATATAATGCTAAAAACGc    <  2:1577123/140‑1 (MQ=255)
                                                                                                                                    cGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTATTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGTGATGAATATAATGCTAAAAACGc    <  2:674086/140‑1 (MQ=255)
                                                                                                                                      tAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGTGATGAATATAATGCTAAAAACGCCa  >  2:470440/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                           
GTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTCTGGAGCGTTTTTCTTTACTGTCAAAATCGCTGGCGATGGACAGCGAGAACCTGATGTTCTCTCTATTTAAGAACGGTAAACCGGTGACCGGTGATGAATATAATGCTAAAAACGCCA  >  NC_000913/875653‑875926

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: