Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,870,998 (T)6→5 coding (614/1476 nt) yeaI → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,870,9930T.96.5% 248.4 / 2.0 58coding (609/1476 nt)yeaIputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (29/27);  total (30/28)
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.58e-01

TTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGC  >  NC_000913/1870896‑1871088
                                                                                                  |                                                                                               
tttGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCttt                                                                                               >  2:718667/1‑100 (MQ=255)
 ttGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttt                                                                                              <  2:459687/100‑1 (MQ=255)
  tGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCttttt                                                                                             >  1:442479/1‑100 (MQ=255)
  tGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCttttt                                                                                             >  1:216217/1‑100 (MQ=255)
   gattAACCTCGCTGGCGCTCTTTTTTTATGGAAAAGACAACATCCTTTACAACAATAAGAAAATCACGGGAATCCTGTTGCTGGCGCTGATCCCttttt                                                                                             <  2:617864/97‑1 (MQ=255)
    tttAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACGACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttttag                                                                                           >  2:699016/1‑98 (MQ=255)
    tttAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttttag                                                                                           >  1:68236/1‑98 (MQ=255)
         ccTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGttttt                                                                                     >  1:451550/1‑100 (MQ=255)
            cGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTccc                                                                                  <  1:1017583/100‑1 (MQ=255)
             gCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTcctc                                                                                 >  2:793428/1‑98 (MQ=255)
              cTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTccccn                                                                                <  1:1098271/100‑2 (MQ=255)
                   gcTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTgg                                                                           <  1:721716/100‑1 (MQ=255)
                   gcTCTTTTGTTATGGAAAAGACAACATCCTTGAAAAAAAAAAAAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTgg                                                                           <  2:364531/100‑1 (MQ=255)
                    ctctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGc                                                                          <  2:969841/100‑1 (MQ=255)
                    ctctTTTGTTATGGAAAAGACAACATCCTTAACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGc                                                                          <  1:46262/100‑1 (MQ=255)
                      ctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGca                                                                         >  1:981630/1‑100 (MQ=255)
                      ctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGccc                                                                        >  2:585989/1‑98 (MQ=255)
                      ctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGcac                                                                        >  1:830269/1‑100 (MQ=255)
                      ctTTTGTTATGGAAAAGACAACATCCTTGACAAAAAAAAAAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGc                                                                          <  2:1064391/98‑1 (MQ=255)
                         ttGTTATGGAAAAGACAACATCCTTGACAACAATATAAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGcacaca                                                                     <  2:289225/100‑1 (MQ=255)
                          tGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACaa                                                                    >  2:888156/1‑100 (MQ=255)
                           gTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAAt                                                                   >  1:322718/1‑100 (MQ=255)
                            ttATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATc                                                                  <  2:559943/100‑1 (MQ=255)
                            ttATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATc                                                                  <  2:303137/100‑1 (MQ=255)
                              aTGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCtn                                                                <  2:528238/100‑2 (MQ=255)
                              aTGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTg                                                                >  2:275913/1‑100 (MQ=255)
                                ggAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAg                                                              >  1:404105/1‑100 (MQ=255)
                                 gAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACACTCTGAGc                                                             >  2:1110106/1‑100 (MQ=255)
                                 gAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGc                                                             >  1:29694/1‑100 (MQ=255)
                                   aaaGACAACCTCCTTGACAACAATAAGAAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCa                                                            <  1:572383/100‑1 (MQ=255)
                                     aGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGtt                                                         <  1:720131/100‑1 (MQ=255)
                                      gACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGtt                                                         >  1:351377/1‑99 (MQ=255)
                                      gACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTta                                                        <  1:1050681/100‑1 (MQ=255)
                                      gACAACAACCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTta                                                        <  2:1134168/100‑1 (MQ=255)
                                       acaacaTCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTtat                                                       >  1:713440/1‑100 (MQ=255)
                                             tCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCt                                                 >  1:15851/1‑100 (MQ=255)
                                             tCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAACCTGAGCAGTTATAACGCt                                                 >  1:674233/1‑100 (MQ=255)
                                                ttGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGAc                                              <  2:163470/100‑1 (MQ=255)
                                                   acaacaATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTAt                                           >  2:415566/1‑100 (MQ=255)
                                                   aaaacaaAAAAAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTa                                            <  2:101989/97‑1 (MQ=255)
                                                     aacaaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTc                                         <  1:608864/100‑1 (MQ=255)
                                                      acaaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCt                                        >  1:1104030/1‑100 (MQ=255)
                                                         aTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTg                                     <  2:1064957/100‑1 (MQ=255)
                                                         aTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTg                                     <  2:1171266/100‑1 (MQ=255)
                                                           aaGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTa                                   <  1:1030626/100‑1 (MQ=255)
                                                           aaGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTa                                   >  2:1030617/1‑100 (MQ=255)
                                                            aGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTAt                                  >  2:909716/1‑100 (MQ=255)
                                                            aGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTAt                                  >  1:1086977/1‑100 (MQ=255)
                                                             gAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                 >  1:1063830/1‑100 (MQ=255)
                                                             gAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                 <  1:925262/100‑1 (MQ=255)
                                                             gAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                 >  2:587740/1‑100 (MQ=255)
                                                             gAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                 <  1:1133054/100‑1 (MQ=255)
                                                                acaaaaCGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATgtc                               >  1:882618/3‑100 (MQ=255)
                                                                aaaaaaCGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATgtc                               >  2:1187199/1‑100 (MQ=255)
                                                                aaaaaaCGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATgtc                               >  2:88945/1‑100 (MQ=255)
                                                                         gAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTAc                      >  2:172352/1‑100 (MQ=255)
                                                                         gAATCCTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTAc                      >  2:913250/1‑100 (MQ=255)
                                                                              cTGTTGCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTcc                 <  2:1117019/100‑1 (MQ=255)
                                                                               tGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTcc                 <  2:981622/100‑1 (MQ=255)
                                                                                 ttGCTGGCGTTGATCCCTTTATCGTATCTGCTCTTCTCGCACACACTCTCAGACGTACTTACGCTGACTATTCTTTCTATTTCGACGATCACTGTCCg                <  2:760480/98‑1 (MQ=255)
                                                                                   gCTGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACa            >  1:562084/1‑100 (MQ=255)
                                                                                     tGGCGCTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAAc          >  2:477466/1‑100 (MQ=255)
                                                                                        cgcTGATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCAt       >  1:368232/1‑100 (MQ=255)
                                                                                            gATCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTg   <  2:351374/100‑1 (MQ=255)
                                                                                             aTCCCTTTTTAGTTTTTCCTCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGc  <  2:442475/100‑1 (MQ=255)
                                                                                                  |                                                                                               
TTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGC  >  NC_000913/1870896‑1871088

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