Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,870,696 T→G I104M (ATT→ATG 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,6960TG100.0% 203.8 / NA 52I104M (ATT→ATGyeaIputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (29/23);  total (29/23)

TCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATTTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAG  >  NC_000913/1870549‑1870843
                                                                                                                                                   |                                                                                                                                                   
tCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCAtgttt                                                                                                                                                  >  1:642329/1‑151 (MQ=255)
tCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCAtgttt                                                                                                                                                  >  2:788710/1‑151 (MQ=255)
tCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCAtgttt                                                                                                                                                  >  1:224765/1‑151 (MQ=255)
  aaaGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATgtttct                                                                                                                                                >  1:714507/1‑151 (MQ=255)
   aaGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATgtttctt                                                                                                                                               >  1:876628/1‑151 (MQ=255)
   aaGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATgtttctt                                                                                                                                               >  2:165058/1‑151 (MQ=255)
           ttATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGt                                                                                                                                        <  2:900683/150‑1 (MQ=255)
            tATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGtt                                                                                                                                       <  2:882577/150‑1 (MQ=255)
              ttttCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAAt                                                                                                                                    >  2:412355/1‑151 (MQ=255)
                ttCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATg                                                                                                                                   <  2:708533/150‑1 (MQ=255)
                ttCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATg                                                                                                                                   >  1:108727/1‑150 (MQ=255)
                 tCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGa                                                                                                                                  >  1:200044/1‑150 (MQ=255)
                       tATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATAtt                                                                                                                            <  1:279579/150‑1 (MQ=255)
                           gCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCAc                                                                                                                        <  2:594378/150‑1 (MQ=255)
                              gACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTaac                                                                                                                    >  1:884523/1‑151 (MQ=255)
                              gACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTa                                                                                                                      <  2:442313/149‑1 (MQ=255)
                               aCGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTaac                                                                                                                    >  2:614944/1‑150 (MQ=255)
                                     attttattttGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCATCATGTTTCTTTTCATGTTAATGAAATATTTCACTAAcaaacaaa                                                                                                             <  2:8445/151‑1 (MQ=255)
                                        ttattttGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAAttt                                                                                                          >  2:303827/1‑151 (MQ=255)
                                          attttGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAAtttt                                                                                                         <  2:186171/150‑1 (MQ=255)
                                            tttGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTAc                                                                                                       <  1:393405/150‑1 (MQ=255)
                                               gCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAac                                                                                                   >  2:803798/1‑151 (MQ=255)
                                                   tactacCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTc                                                                                               >  1:2667/1‑151 (MQ=255)
                                                      tacCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCAtt                                                                                            >  1:320038/1‑151 (MQ=255)
                                                       acCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCAttt                                                                                           <  1:901937/151‑1 (MQ=255)
                                                        ccTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCAtttt                                                                                          >  2:847122/1‑151 (MQ=255)
                                                         cTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCAtttt                                                                                          >  1:670944/1‑150 (MQ=255)
                                                           tatataACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTa                                                                                         >  2:912284/1‑149 (MQ=255)
                                                           tatataACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTa                                                                                         >  2:348495/1‑149 (MQ=255)
                                                            atataACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGc                                                                                      >  2:624917/1‑151 (MQ=255)
                                                               taACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCtt                                                                                   >  2:163522/1‑151 (MQ=255)
                                                                  cGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCa                                                                                <  1:797540/151‑1 (MQ=255)
                                                                    aaaaCACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCAtt                                                                              <  2:768118/151‑1 (MQ=255)
                                                                      aaCACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCAtttt                                                                            >  1:232243/1‑151 (MQ=255)
                                                                       acacAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCAttttt                                                                           <  1:36761/151‑1 (MQ=255)
                                                                                tAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCt                                                                  >  2:354479/1‑151 (MQ=255)
                                                                                 aaaaCTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCtt                                                                 <  1:165058/151‑1 (MQ=255)
                                                                                      tCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCt                                                             >  1:274144/1‑150 (MQ=255)
                                                                                       cATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTa                                                            <  2:916471/150‑1 (MQ=255)
                                                                                        aTCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTACTGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTAtt                                                          >  1:266896/1‑151 (MQ=255)
                                                                                                 aTGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGaa                                                  >  2:701466/1‑150 (MQ=255)
                                                                                                       tACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCAtt                                            <  1:617700/150‑1 (MQ=255)
                                                                                                         ccTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGt                                          >  2:203240/1‑150 (MQ=255)
                                                                                                            gATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAAtt                                      >  1:158791/1‑151 (MQ=255)
                                                                                                               ggTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTAt                                    <  2:391540/150‑1 (MQ=255)
                                                                                                                         gCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCaa                         <  2:836625/151‑1 (MQ=255)
                                                                                                                         gCTTTGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCa                          >  1:216196/1‑150 (MQ=255)
                                                                                                                           tttGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAAtt                       >  2:858934/1‑151 (MQ=255)
                                                                                                                            ttGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAAtt                       >  2:26237/1‑150 (MQ=255)
                                                                                                                             tGTTTTTCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTa                      <  2:25591/150‑1 (MQ=255)
                                                                                                                                  ttCCTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGc                 <  2:383996/150‑1 (MQ=255)
                                                                                                                                     cTTAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGt              >  1:518666/1‑150 (MQ=255)
                                                                                                                                      ttAACGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTac            >  1:788695/1‑151 (MQ=255)
                                                                                                                                          cGCCGTCATGTTTCTTTTCATGTTAATGAAATATTTAACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTg        <  1:223207/151‑1 (MQ=255)
                                                                                                                                           gCCGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTg        >  2:218382/1‑150 (MQ=255)
                                                                                                                                             cGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATc     <  2:624828/151‑1 (MQ=255)
                                                                                                                                             cGTCATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATc     >  2:920107/1‑151 (MQ=255)
                                                                                                                                                cATGTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAg  <  2:884503/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATTTTTCTTTTCATGTTAATGAAATATTTCACTAACAAACAAATTTTACCAACACTCATTTTAAGCCTTGCATTTTTAAGTGGCCTTATCTATTTAGTTGAAACCATTGTAATTATCCATAAACCAATTAACGGCAGTACACTGATCCAG  >  NC_000913/1870549‑1870843

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