Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 875,714 T→G V127G (GTC→GGC)  yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,7140TG88.1% 115.2 / 1.3 44V127G (GTC→GGC) yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  major base G (18/19);  minor base A (0/5);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 6.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTC  >  NC_000913/875618‑875809
                                                                                                |                                                                                               
gTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCaa                                                                                              >  2:693038/1‑100 (MQ=255)
  ccTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCaaaa                                                                                            <  2:824200/100‑1 (MQ=255)
         tGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAgg                                                                                     <  2:173003/100‑1 (MQ=255)
             ttCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGttt                                                                                 <  2:498203/100‑1 (MQ=255)
              tCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTa                                                                                <  1:216015/100‑1 (MQ=255)
                cAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTAtc                                                                              >  2:936077/1‑100 (MQ=255)
                cAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTAtc                                                                              <  1:226195/100‑1 (MQ=255)
                 aGTGACGAATGTATAAAGGTAAAGGAGTTCGAGAAAAATAATTGTGAAGATATTATGCAGGTTGTGAAGTAAAAAAAAAAAAAAAATACAGGGTTTAtct                                                                             <  2:943135/100‑1 (MQ=255)
                    gACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCtt                                                                          >  1:604831/1‑100 (MQ=255)
                     aCGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCttt                                                                         <  1:705515/100‑1 (MQ=255)
                      cGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCtttt                                                                        >  1:894389/1‑100 (MQ=255)
                      cGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCtttt                                                                        <  1:521593/100‑1 (MQ=255)
                       gAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATAAAAAAAAAAAAAACACAGGGTTTATCTCTTTTg                                                                       <  2:928103/100‑1 (MQ=255)
                             aTTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTa                                                                 <  1:938184/100‑1 (MQ=255)
                             aTTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTa                                                                 <  2:24550/100‑1 (MQ=255)
                             aTTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTa                                                                 >  2:729704/1‑100 (MQ=255)
                                 aGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAAc                                                             <  1:727582/100‑1 (MQ=255)
                                    tAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCtt                                                          >  1:560622/1‑100 (MQ=255)
                                    tAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCtt                                                          >  1:1267442/1‑100 (MQ=255)
                                      aaGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTAAAAAAAAAAAAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCg                                                        <  2:886960/100‑1 (MQ=255)
                                           gTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTAt                                                   >  2:742735/1‑100 (MQ=255)
                                           gTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTAt                                                   >  2:1285781/1‑100 (MQ=255)
                                             tCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTa                                                 <  1:1290020/100‑1 (MQ=255)
                                               gagaAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTAcc                                               <  2:471635/100‑1 (MQ=255)
                                                 gaAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACctc                                             >  2:1075761/1‑100 (MQ=255)
                                                     aaTTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATc                                         >  1:1178577/1‑100 (MQ=255)
                                                        tattGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCgg                                      <  1:374472/100‑1 (MQ=255)
                                                          ttGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTa                                    <  1:396212/100‑1 (MQ=255)
                                                           tgtgCAGATATTATGCAGGTTGTCAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTaa                                   >  2:290196/1‑100 (MQ=255)
                                                            gtgCAGATATTATGCAGGTTGTGAAGTATCGACATAGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAAc                                  <  2:440210/100‑1 (MQ=255)
                                                               cAGATATTATGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTca                               <  2:894382/100‑1 (MQ=255)
                                                                aGATATTATGCAGGTTGTGAAGTAAAAAAAAAAAAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTcac                              <  2:149284/100‑1 (MQ=255)
                                                                     ttatGCAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCg                         >  1:189500/1‑100 (MQ=255)
                                                                      tatGCAGGTTGTGAAGTAAAAAAAAAAAAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGt                        <  2:1110900/100‑1 (MQ=255)
                                                                          cAGGTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTc                    >  2:1139614/1‑100 (MQ=255)
                                                                            ggTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTg                  >  1:817908/1‑100 (MQ=255)
                                                                            ggTTGTGAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTg                  <  2:44986/100‑1 (MQ=255)
                                                                                gtgAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGatat              >  1:147137/1‑100 (MQ=255)
                                                                                 tgAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATAtt             >  1:617859/1‑100 (MQ=255)
                                                                                  gAAGTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATAttt            >  1:701659/1‑100 (MQ=255)
                                                                                    aGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGatgcgaat          >  1:425847/1‑94 (MQ=255)
                                                                                     gTATCGACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCACTATCCGGTAACTCACAATCGTAGTCTGATATTTTTg         <  2:700524/100‑1 (MQ=255)
                                                                                          gACATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGAGATTTTGGCTTgg    >  2:106979/1‑100 (MQ=255)
                                                                                            cATGGCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTc  <  1:1169250/100‑1 (MQ=255)
                                                                                                |                                                                                               
GTCCTTATCTGATTTCCAGTGACGAATGTATTCAGGTAAAGGAGTTCGAGAAAAATTATTGTGCAGATATTATGCAGGTTGTGAAGTATCGACATGTCAAAAATACAGGGTTTATCTCTTTTGACGGTAAAACCTTCGTCTATTACCTCTATCCGGTAACTCACAATCGTAGTCTGATATTTTTGCTTGGTC  >  NC_000913/875618‑875809

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