Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 404,123 (C)5→4 coding (421/1116 nt) yaiC → diguanylate cyclase, cellulose regualtor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913404,1190C.97.8% 193.9 / ‑1.8 46coding (417/1116 nt)yaiCdiguanylate cyclase, cellulose regualtor
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (27/18);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

GCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCAT  >  NC_000913/403984‑404250
                                                                                                                                       |                                                                                                                                   
gCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGcccc                                                                                                                                  >  2:661298/1‑139 (MQ=255)
                 aCAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGt                                                                                                                <  1:447161/139‑1 (MQ=255)
                 aCAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTACTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGt                                                                                                                <  1:46221/139‑1 (MQ=255)
                    aCTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTg                                                                                                             <  2:366261/139‑1 (MQ=255)
                                      tATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGttt                                                                                             <  2:212610/137‑1 (MQ=255)
                                      tATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGttt                                                                                             >  1:212610/1‑137 (MQ=255)
                                             gggAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTg                                                                                    <  2:88569/139‑1 (MQ=255)
                                                           gCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTg                                                                      >  2:334691/1‑139 (MQ=255)
                                                              tAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGAc                                                                    <  1:8437/139‑1 (MQ=255)
                                                              tAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTcgc                                                                                                                    <  1:510974/90‑1 (MQ=255)
                                                              tAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTcgc                                                                                                                    >  2:510974/1‑90 (MQ=255)
                                                                aTGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACggg                                                                 >  1:152471/1‑139 (MQ=255)
                                                                 tGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATCGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGc                                                                >  2:506894/1‑139 (MQ=255)
                                                                     cGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTa                                                            >  1:198101/1‑139 (MQ=255)
                                                                     cGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTa                                                            >  2:724440/1‑139 (MQ=255)
                                                                     cGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTa                                                            >  1:747325/1‑139 (MQ=255)
                                                                               ctgcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGtt                                                  <  1:724440/139‑1 (MQ=255)
                                                                                            cgATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGt                                                                                                                >  2:789516/1‑64 (MQ=255)
                                                                                            cgATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGt                                                                                                                <  1:789516/64‑1 (MQ=255)
                                                                                            cgATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGcgccg                                     >  2:350107/1‑139 (MQ=255)
                                                                                            cgATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGcgccg                                     >  1:598270/1‑139 (MQ=255)
                                                                                              aTGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCt                                   >  1:620723/1‑139 (MQ=255)
                                                                                                       atTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAAtggtg                          <  2:231925/139‑1 (MQ=255)
                                                                                                         tATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGc                        <  2:598270/139‑1 (MQ=255)
                                                                                                           tgtgTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctc                      >  2:356458/1‑139 (MQ=255)
                                                                                                           tgtgTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctc                      >  2:115192/1‑139 (MQ=255)
                                                                                                           tgtgTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctc                      >  1:402216/1‑139 (MQ=255)
                                                                                                           tgtgTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctc                      >  2:308335/1‑139 (MQ=255)
                                                                                                           tgtgTCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGccc                      >  1:200198/1‑137 (MQ=255)
                                                                                                            gtgtCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGctct                     <  1:115192/139‑1 (MQ=255)
                                                                                                             tgtCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGAc                                                                    >  2:139652/1‑91 (MQ=255)
                                                                                                             tgtCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGAc                                                                    <  1:139652/91‑1 (MQ=255)
                                                                                                              gtCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGATTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTcc                   >  1:518035/1‑139 (MQ=255)
                                                                                                               tCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCgtg                                                       >  1:54136/1‑102 (MQ=255)
                                                                                                               tCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCgtg                                                       <  2:54136/102‑1 (MQ=255)
                                                                                                               tCTGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTccc                  >  1:484632/1‑139 (MQ=255)
                                                                                                                 tGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGa                                                                     <  2:104272/86‑1 (MQ=255)
                                                                                                                 tGAATTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGa                                                                     >  1:104272/1‑86 (MQ=255)
                                                                                                                   aaTTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTcg                                                                                                                     <  1:627317/36‑1 (MQ=255)
                                                                                                                   aaTTTGATGGGGGCAGGCGGCCCCGTCTGTTTGTcg                                                                                                                     >  2:627317/1‑36 (MQ=255)
                                                                                                                             gggCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGttt                                                                                             <  2:290338/50‑1 (MQ=255)
                                                                                                                             gggCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGttt                                                                                             >  1:290338/1‑50 (MQ=255)
                                                                                                                             gggCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGGTGa                                                                     <  1:517149/74‑1 (MQ=255)
                                                                                                                             gggCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGa                                                                     >  2:517149/1‑74 (MQ=255)
                                                                                                                             gggCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTc    >  2:287877/1‑139 (MQ=255)
                                                                                                                              ggCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCa   >  2:178406/1‑139 (MQ=255)
                                                                                                                               gCAGGCGGCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCAt  <  1:350107/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
GCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCAT  >  NC_000913/403984‑404250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: