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

TTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGC  >  NC_000913/1870896‑1871088
                                                                                                  |                                                                                               
tttGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCttt                                                                                               >  2:718667/1‑100 (MQ=255)
 ttGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttt                                                                                              <  2:459687/100‑1 (MQ=255)
  tGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCttttt                                                                                             >  1:442479/1‑100 (MQ=255)
  tGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCttttt                                                                                             >  1:216217/1‑100 (MQ=255)
   gattAACCTCGCTGGCGCTCTTTTTTTATGGAAAAGACAACATCCTTTACAACAATAAG‑AAAATCACGGGAATCCTGTTGCTGGCGCTGATCCCttttt                                                                                             <  2:617864/97‑1 (MQ=255)
    tttAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACGACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttttag                                                                                           >  2:699016/1‑98 (MQ=255)
    tttAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttttag                                                                                           >  1:68236/1‑98 (MQ=255)
         ccTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGttttt                                                                                     >  1:451550/1‑100 (MQ=255)
            cGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTccc                                                                                  <  1:1017583/100‑1 (MQ=255)
             gCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTcctc                                                                                 >  2:793428/1‑98 (MQ=255)
              cTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTccccn                                                                                <  1:1098271/100‑2 (MQ=255)
                   gcTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTgg                                                                           <  1:721716/100‑1 (MQ=255)
                   gcTCTTTTGTTATGGAAAAGACAACATCCTTGAAAAAAAAAAA‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTgg                                                                           <  2:364531/100‑1 (MQ=255)
                    ctctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGc                                                                          <  2:969841/100‑1 (MQ=255)
                    ctctTTTGTTATGGAAAAGACAACATCCTTAACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGc                                                                          <  1:46262/100‑1 (MQ=255)
                      ctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGca                                                                         >  1:981630/1‑100 (MQ=255)
                      ctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGccc                                                                        >  2:585989/1‑98 (MQ=255)
                      ctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGcac                                                                        >  1:830269/1‑100 (MQ=255)
                      ctTTTGTTATGGAAAAGACAACATCCTTGACAAAAAAAAA‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGc                                                                          <  2:1064391/98‑1 (MQ=255)
                         ttGTTATGGAAAAGACAACATCCTTGACAACAATATA‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGcacaca                                                                     <  2:289225/100‑1 (MQ=255)
                          tGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACaa                                                                    >  2:888156/1‑100 (MQ=255)
                           gTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAAt                                                                   >  1:322718/1‑100 (MQ=255)
                            ttATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATc                                                                  <  2:559943/100‑1 (MQ=255)
                            ttATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATc                                                                  <  2:303137/100‑1 (MQ=255)
                              aTGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCtn                                                                <  2:528238/100‑2 (MQ=255)
                              aTGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTg                                                                >  2:275913/1‑100 (MQ=255)
                                ggAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAg                                                              >  1:404105/1‑100 (MQ=255)
                                 gAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACACTCTGAGc                                                             >  2:1110106/1‑100 (MQ=255)
                                 gAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGc                                                             >  1:29694/1‑100 (MQ=255)
                                   aaaGACAACCTCCTTGACAACAATAAGAAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCa                                                            <  1:572383/100‑1 (MQ=255)
                                     aGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGtt                                                         <  1:720131/100‑1 (MQ=255)
                                      gACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGtt                                                         >  1:351377/1‑99 (MQ=255)
                                      gACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTta                                                        <  1:1050681/100‑1 (MQ=255)
                                      gACAACAACCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTta                                                        <  2:1134168/100‑1 (MQ=255)
                                       acaacaTCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTtat                                                       >  1:713440/1‑100 (MQ=255)
                                             tCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCt                                                 >  1:15851/1‑100 (MQ=255)
                                             tCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAACCTGAGCAGTTATAACGCt                                                 >  1:674233/1‑100 (MQ=255)
                                                ttGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGAc                                              <  2:163470/100‑1 (MQ=255)
                                                   acaacaATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTAt                                           >  2:415566/1‑100 (MQ=255)
                                                   aaaacaaAAAA‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTa                                            <  2:101989/97‑1 (MQ=255)
                                                     aacaaTAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTc                                         <  1:608864/100‑1 (MQ=255)
                                                      acaaTAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCt                                        >  1:1104030/1‑100 (MQ=255)
                                                         aTAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTg                                     <  2:1064957/100‑1 (MQ=255)
                                                         aTAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTg                                     <  2:1171266/100‑1 (MQ=255)
                                                           aaG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTa                                   <  1:1030626/100‑1 (MQ=255)
                                                           aaG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTa                                   >  2:1030617/1‑100 (MQ=255)
                                                            aG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTAt                                  >  2:909716/1‑100 (MQ=255)
                                                            aG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTAt                                  >  1:1086977/1‑100 (MQ=255)
                                                             g‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                 >  1:1063830/1‑100 (MQ=255)
                                                             g‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                 <  1:925262/100‑1 (MQ=255)
                                                             g‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                 >  2:587740/1‑100 (MQ=255)
                                                             g‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                 <  1:1133054/100‑1 (MQ=255)
                                                                acaaaaCGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATgtc                               >  1:882618/3‑100 (MQ=255)
                                                                aaaaaaCGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATgtc                               >  2:1187199/1‑100 (MQ=255)
                                                                aaaaaaCGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATgtc                               >  2:88945/1‑100 (MQ=255)
                                                                         gAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTAc                      >  2:172352/1‑100 (MQ=255)
                                                                         gAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTAc                      >  2:913250/1‑100 (MQ=255)
                                                                              cTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTcc                 <  2:1117019/100‑1 (MQ=255)
                                                                               tGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTcc                 <  2:981622/100‑1 (MQ=255)
                                                                                 ttGCTGGCGTTGATCCC‑TTTATCGTATCTGCTCTTCTCGCACACACTCTCAGACGTACTTACGCTGACTATTCTTTCTATTTCGACGATCACTGTCCg                <  2:760480/98‑1 (MQ=255)
                                                                                   gCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACa            >  1:562084/1‑100 (MQ=255)
                                                                                     tGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAAc          >  2:477466/1‑100 (MQ=255)
                                                                                        cgcTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCAt       >  1:368232/1‑100 (MQ=255)
                                                                                            gATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTg   <  2:351374/100‑1 (MQ=255)
                                                                                             aTCCC‑TTTTTAGTTTTTCCTCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGc  <  2:442475/100‑1 (MQ=255)
                                                                                                  |                                                                                               
TTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAG‑AAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGC  >  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: ‑